snapshot...

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)