diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 8602cf378bd8206607cd1efd5ce1d7f8e2c91afe..52f320f2b4c5b3054a43672443ebccc9434e6af4 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -1229,8 +1229,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd # if not (is_day(solzen[i,]) or is_night(solzen[i,])): # continue - if not (icing_alt[i] < max_altitude): - continue + # if not (icing_alt[i] < max_altitude): + # continue keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy keep_1 = np.invert(np.isnan(cld_top_hgt[i,])) @@ -1243,13 +1243,13 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd if num_keep / nx_x_ny < cloud_frac: # At least this fraction cloudy continue - ##keep = np.where(keep, (cld_top_hgt[i,] + closeness_top) > icing_alt[i], False) - keep = np.where(keep, cld_top_hgt[i,] < max_cld_altitude, False) + keep = np.where(keep, (cld_top_hgt[i,] + closeness_top) > icing_alt[i], False) + #keep = np.where(keep, cld_top_hgt[i,] < max_cld_altitude, False) ##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, cld_top_temp[i,] < 274.2, False) # cld_hgt = cld_top_hgt[i, ].flatten() # med_cld_hgt = np.median(cld_hgt[keep]) @@ -1307,8 +1307,8 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_ # if not (is_day(solzen[i,]) or is_night(solzen[i,])): # continue - if not (icing_alt[i] < max_altitude): - continue + # if not (icing_alt[i] < max_altitude): + # continue keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy keep_1 = np.invert(np.isnan(cld_top_hgt[i,])) @@ -1321,13 +1321,13 @@ def apply_qc_no_icing_pireps(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 - ##keep = np.where(keep, (cld_top_hgt[i,] + closeness_top) > icing_alt[i], False) - keep = np.where(keep, cld_top_hgt[i,] < max_cld_altitude, False) + keep = np.where(keep, (cld_top_hgt[i,] + closeness_top) > icing_alt[i], False) + #keep = np.where(keep, cld_top_hgt[i,] < max_cld_altitude, False) ##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, cld_top_temp[i,] < 274.2, False) # cld_hgt = cld_top_hgt[i, ].flatten() # med_cld_hgt = np.median(cld_hgt[keep])