snapshot...

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