diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 84d0c872bf5deb73bc7d235968cd6616672d4e53..3268dea404fea54d53823022dc95d8cb1b30f7c3 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -768,8 +768,10 @@ def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20, icing def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, bt_11um, solzen, satzen, day_night='ANY'): opd_thick_threshold = 2 + opd_thin_threshold = 0.1 if day_night == 'DAY': opd_thick_threshold = 20 + opd_thin_threshold = 1 closeness = 100.0 # meters num_obs = len(icing_alt) cld_mask = cld_mask.reshape((num_obs, -1)) @@ -810,7 +812,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, keep = np.where(keep, (cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4) & (cld_top_hgt[i,] > icing_alt[i]), False) - keep = np.where(keep, np.invert((cld_phase[i,] == 4) & (cld_opd[i,] < 0.1) & (cld_top_hgt[i,] > icing_alt[i])), False) + keep = np.where(keep, np.invert((cld_phase[i,] == 4) & (cld_opd[i,] < opd_thin_threshold) & (cld_top_hgt[i,] > icing_alt[i])), False) keep = np.where(keep, np.invert(bt_11um[i,] > 270.0), False) @@ -870,9 +872,9 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mas # # keep = np.where(keep, np.invert((cld_phase[i,] == 4) & (cld_opd[i,] < 0.1) & (cld_top_hgt[i,] > icing_alt[i])), False) # - # keep = np.where(keep, np.invert(bt_11um[i,] > 270.0), False) - # - # keep = np.where(keep, np.invert(bt_11um[i,] < 228.0), False) + keep = np.where(keep, bt_11um[i,] > 270.0, False) + + keep = np.where(keep, bt_11um[i,] < 228.0, False) mask.append(keep) idxs.append(i)