diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 08f1e15a2b8d3bac5998ce55057fb3bef8c6ab14..4f35c749e06e7f4c4932276594ac2e7515c844f4 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -873,6 +873,21 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, # Exp # keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False) + # Test1 + keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False) + + # 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) + + # 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) @@ -929,6 +944,9 @@ 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) + # Test1 + keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], 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)