snapshot...

modules/icing/pirep_goes.py

@@ -1277,8 +1277,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
         keep = np.where(keep, (cld_top_hgt[i,] - max_cld_depth) < icing_alt[i], 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,] < 275.0, False)
-        keep = np.where(keep, cld_top_temp[i,] < 275.0, False)
+        keep = np.where(keep, bt_11um[i,] < 275.0, False)
+        # keep = np.where(keep, cld_top_temp[i,] < 275.0, False)
 
         # cld_hgt = cld_top_hgt[i, ].flatten()
         # med_cld_hgt = np.median(cld_hgt[keep])
@@ -1310,7 +1310,7 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
         opd_thin_threshold = 0.1
 
     closeness_top = 100.0  # meters
-    max_cld_depth = 1500.0
+    max_cld_depth = 2000.0
     max_altitude = 5000.0
     max_cld_altitude = 5000.0
 
@@ -1355,8 +1355,8 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
         keep = np.where(keep, (cld_top_hgt[i,] - max_cld_depth) < icing_alt[i], 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,] < 275.0, False)
-        keep = np.where(keep, cld_top_temp[i,] < 275.0, False)
+        keep = np.where(keep, bt_11um[i,] < 275.0, False)
+        # keep = np.where(keep, cld_top_temp[i,] < 275.0, False)
 
         # cld_hgt = cld_top_hgt[i, ].flatten()
         # med_cld_hgt = np.median(cld_hgt[keep])