snapshot...

modules/icing/pirep_goes.py

@@ -1466,7 +1466,7 @@ def apply_qc_icing_pireps_exp2(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cl
 def apply_qc_icing_pireps_v3(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd, cld_mask, bt_11um, solzen, satzen, cld_top_temp, day_night='ANY', cloudy_frac=0.5):
 
     closeness_top = 50.0  # meters
-    max_cld_depth = 800.0
+    max_cld_depth = 500.0
     max_altitude = 4500.0
     max_cld_altitude = 4000.0
 
@@ -1493,9 +1493,10 @@ def apply_qc_icing_pireps_v3(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
         # keep_0 = cld_mask[i, ] == 3  # Confident cloudy
         keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3)  # cloudy
         keep_1 = np.invert(np.isnan(cld_top_hgt[i,]))
-        keep_2 = np.invert(np.isnan(bt_11um[i,]))
-        keep_3 = np.invert(np.isnan(cld_geo_dz[i,]))
-        keep = keep_0 & keep_1 & keep_2 & keep_3
+        # keep_2 = np.invert(np.isnan(bt_11um[i,]))
+        # keep_3 = np.invert(np.isnan(cld_geo_dz[i,]))
+        # keep = keep_0 & keep_1 & keep_2 & keep_3
+        keep = keep_0 & keep_1
         num_keep = np.sum(keep)
         if num_keep == 0:
             continue