diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index a773a6f41307eeac7f7ab21d90071e29cf709e87..5bc168ead1af051232e2be932bc01eaf707e2c92 100644 --- a/modules/icing/pirep_goes.py +++ b/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,],