diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index db7b1864106b4693be173c990cac95cd78cdd39d..c84b0ed500584fb203914767339211d1282aa37e 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -944,23 +944,24 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, # Test 1 keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False) - if (np.sum(keep) / num_keep) < 0.80: - continue + # 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) + # 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)), False) + # keep = np.where(keep, np.invert((cld_phase[i,] == 4) & (cld_opd[i,] < opd_thin_threshold)), False) # Test5 and Test6 - keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 270.0), False) + # keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 270.0), False) + keep = np.where(keep, bt_11um[i,] < 270.0, False) # Test3 - keep = np.where(keep, (cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4), False) + # keep = np.where(keep, (cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4), False) mask.append(keep) idxs.append(i)