diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index e159744b154236c0268bbdcb3b3f4df6844666d5..c99e444c60a885626edd5b581773cf5996f3620b 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -294,15 +294,14 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, bt_11um, c keep = keep_0 & keep_1 & keep_2 & keep_3 keep = np.where(keep, cld_top_hgt[:, ] > icing_alt[:], False) - keep = np.invert( - np.where(keep, - np.logical_and(cld_phase[:, ] == 3, + keep = np.where(keep, + np.invert((cld_phase[:, ] == 3) & np.logical_and(cld_top_hgt[:, ]+closeness > icing_alt[:], cld_top_hgt[:, ]-closeness < icing_alt[:])), - False)) + False) keep = np.where(keep, (cld_opd[:, ] > opd_threshold) & (cld_phase[:, ] == 3) & (cld_top_hgt[:, ] > icing_alt[:]), False) keep = np.where(keep, np.invert((cld_phase[:, ] == 3) & (cld_opd[:,] < 0.1) & (cld_top_hgt[:, ] > icing_alt[:]), False)) - keep = np.where(keep, np.invert(bt_11um[:,] > 270.0), False) - keep = np.where(keep, np.invert(bt_11um[:,] < 228.0), False) + keep = np.where(keep, np.invert(bt_11um[:, ] > 270.0), False) + keep = np.where(keep, np.invert(bt_11um[:, ] < 228.0), False) - return keep \ No newline at end of file + return keep