diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index e637b0e418279e9e4d792844862382081787595a..9fd0ea3ecc7f9ceefd0e9c99b36cd192e8a3a393 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -1041,7 +1041,6 @@ def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20, icing keep_idxs = [] for i in range(len(mask)): - keep_idxs.append(idxs[i]) # frac = np.sum(mask[i]) / nx_x_ny frac = np.sum(mask[i]) / num_tested[i] if icing: @@ -1049,6 +1048,7 @@ def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20, icing keep_idxs.append(idxs[i]) elif frac > pass_thresh_frac: keep_idxs.append(idxs[i]) + # keep_idxs.append(idxs[i]) print('day_night, icing, all, valid, pass: ', day_night, icing, len(icing_alt), len(mask), len(keep_idxs)) print('-----------------------') @@ -1277,8 +1277,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd keep = np.where(keep, (cld_top_hgt[i,] - max_cld_depth) < icing_alt[i], False) # keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 270.0), False) - keep = np.where(keep, bt_11um[i,] < 275.0, False) - # keep = np.where(keep, cld_top_temp[i,] < 275.0, False) + # keep = np.where(keep, bt_11um[i,] < 275.0, False) + keep = np.where(keep, cld_top_temp[i,] < 275.0, False) # cld_hgt = cld_top_hgt[i, ].flatten() # med_cld_hgt = np.median(cld_hgt[keep]) @@ -1355,8 +1355,8 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_ keep = np.where(keep, (cld_top_hgt[i,] - max_cld_depth) < icing_alt[i], False) # keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 270.0), False) - keep = np.where(keep, bt_11um[i,] < 275.0, False) - # keep = np.where(keep, cld_top_temp[i,] < 275.0, False) + # keep = np.where(keep, bt_11um[i,] < 275.0, False) + keep = np.where(keep, cld_top_temp[i,] < 275.0, False) # cld_hgt = cld_top_hgt[i, ].flatten() # med_cld_hgt = np.median(cld_hgt[keep])