diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 248105c5f3e34a37368326cf070b0a286e2323ba..e637b0e418279e9e4d792844862382081787595a 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -1277,8 +1277,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd 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, bt_11um[i,] < 275.0, False) + # keep = np.where(keep, cld_top_temp[i,] < 275.0, False) # cld_hgt = cld_top_hgt[i, ].flatten() # med_cld_hgt = np.median(cld_hgt[keep]) @@ -1310,7 +1310,7 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_ opd_thin_threshold = 0.1 closeness_top = 100.0 # meters - max_cld_depth = 1500.0 + max_cld_depth = 2000.0 max_altitude = 5000.0 max_cld_altitude = 5000.0 @@ -1355,8 +1355,8 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_ 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, bt_11um[i,] < 275.0, False) + # keep = np.where(keep, cld_top_temp[i,] < 275.0, False) # cld_hgt = cld_top_hgt[i, ].flatten() # med_cld_hgt = np.median(cld_hgt[keep])