diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 177e8bd2fd12466f4db2e6f9049da28229cb826e..052a3323aaba002142a8c69febdb7f572099a2c0 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -869,23 +869,23 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mas # Test1 keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False) - # Test2 - keep = np.where(keep, - ((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, ((cld_phase[i,] == 4) & (cld_opd[i,] < opd_thin_threshold) & (cld_top_hgt[i,] > icing_alt[i])), False) - - # Test5 - keep = np.where(keep, (bt_11um[i,] > 270.0), False) - - # Test6 - keep = np.where(keep, (bt_11um[i,] < 228.0), False) - - # Test3 - keep = np.where(keep, np.invert((cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4) & (cld_top_hgt[i,] > icing_alt[i])), False) + # # Test2 + # keep = np.where(keep, + # ((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, ((cld_phase[i,] == 4) & (cld_opd[i,] < opd_thin_threshold) & (cld_top_hgt[i,] > icing_alt[i])), False) + # + # # Test5 + # keep = np.where(keep, (bt_11um[i,] > 270.0), False) + # + # # Test6 + # keep = np.where(keep, (bt_11um[i,] < 228.0), False) + # + # # Test3 + # keep = np.where(keep, np.invert((cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4) & (cld_top_hgt[i,] > icing_alt[i])), False) mask.append(keep) idxs.append(i)