diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index ddf4d8cad6cfe05029e8285fcbd544ef2d761582..e0cdd8aee5d1b9b35300c023e59431d194ce62d6 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -943,29 +943,22 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, if num_keep == 0: continue - # Exp + # Test 1 keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False) - #if (np.sum(keep) / num_keep) < 0.50: - # continue - - # Test1 - # keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False) + if (np.sum(keep) / num_keep) < 0.80: + continue # Test2 - # keep = np.where(keep, - # np.invert((cld_phase[i,] == 4) & - # np.logical_and(cld_top_hgt[i,]+closeness > icing_alt[i], cld_top_hgt[i,]-closeness < icing_alt[i])), - # False) - # - # # Test4 - # 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((cld_phase[i,] == 4) & + np.logical_and(cld_top_hgt[i,]+closeness > icing_alt[i], cld_top_hgt[i,]-closeness < icing_alt[i])), + False) - # Test5 - #keep = np.where(keep, (bt_11um[i,] < 270.0), False) + # Test4 + # keep = np.where(keep, np.invert((cld_phase[i,] == 4) & (cld_opd[i,] < opd_thin_threshold) & (cld_top_hgt[i,] > icing_alt[i])), False) - # Test6 - # keep = np.where(keep, (bt_11um[i,] > 228.0), False) - # keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 270.0), False) + # Test5 and Test6 + keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 270.0), False) # Test3 # keep = np.where(keep, (cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4) & (cld_top_hgt[i,] > icing_alt[i]), False) @@ -1021,7 +1014,7 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mas # Exp keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False) - #if (np.sum(keep) / num_keep) < 0.50: + #if (np.sum(keep) / num_keep) < 0.80: # continue # keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 273.0), False)