diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index c7fbe65f43a51535bb4f143ee0e382aa13e800b5..8d81aca98cea183588e1112319ea6166f56c46b7 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -811,17 +811,20 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, np.logical_and(cld_top_hgt[i,]+closeness > icing_alt[i], cld_top_hgt[i,]-closeness < icing_alt[i])), False) # Test3 - 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, (cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4) & (cld_top_hgt[i,] > 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) - # Test5, Should this go before Test2? (Would it flip some keepers from #2) + # Test5 keep = np.where(keep, np.invert(bt_11um[i,] > 270.0), False) # Test6 keep = np.where(keep, np.invert(bt_11um[i,] < 228.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) + mask.append(keep) idxs.append(i) num_tested.append(num_keep)