diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 8d963234cd8fd1b6a98a2005c2b850cc617bf242..06d0841d455f6c970bad70807f3b0917aa316ad5 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -1210,8 +1210,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd opd_thick_threshold = 2 opd_thin_threshold = 0.1 - closeness_top = 200.0 # meters - max_cld_depth = 1500.0 + closeness_top = 100.0 # meters + max_cld_depth = 500.0 max_altitude = 4000.0 max_cld_altitude = 4000.0 @@ -1252,12 +1252,13 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd 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,] - max_cld_depth) < 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, 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, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 274.0), False) # keep = np.where(keep, bt_11um[i,] < 275.0, False) - keep = np.where(keep, cld_top_temp[i,] < 274.2, 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]) @@ -1288,7 +1289,7 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_ opd_thick_threshold = 2 opd_thin_threshold = 0.1 - closeness_top = 200.0 # meters + closeness_top = 100.0 # meters max_cld_depth = 1500.0 max_altitude = 5000.0 max_cld_altitude = 5000.0 @@ -1330,12 +1331,13 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_ 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,] - max_cld_depth) < 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, 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, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 274.0), False) # keep = np.where(keep, bt_11um[i,] < 275.0, False) - keep = np.where(keep, cld_top_temp[i,] < 274.2, 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])