snapshot...

modules/icing/pirep_goes.py

@@ -1397,8 +1397,8 @@ 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 = 200.0  # meters
-    max_cld_depth = 1000.0
-    max_altitude = 4000.0
+    max_cld_depth = 1200.0
+    max_altitude = 4500.0
     max_cld_altitude = 4000.0
 
     num_obs = len(icing_alt)
@@ -1433,8 +1433,15 @@ def apply_qc_icing_pireps_v3(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
         if num_keep / nx_x_ny < cloudy_frac:  # At least this fraction cloudy
             continue
 
+        # ctt_mean = np.nanmean(cld_top_temp[i,])
+        # if ctt_mean > 274.0:
+        #     continue
+
+        # if icing_alt[i,] > max_altitude:
+        #     continue
+
         cth_mean = np.nanmean(cld_top_hgt[i,])
-        if not (closeness_top+cth_mean > icing_alt[i,] > cth_mean - 1200.0):
+        if not (closeness_top+cth_mean > icing_alt[i,] > cth_mean - max_cld_depth):
             continue
 
         # keep = np.where(keep, np.logical_and(cld_top_hgt[i,] > icing_alt[i,],