diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index a4e057470b09d04bb4144945d88dc94279f4568e..d08ac74c1f554ff0c95bf11c683747e4a807dfbc 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 = 100.0 # meters - max_cld_depth = 500.0 + closeness_top = 150.0 # meters + max_cld_depth = 700.0 max_altitude = 4000.0 max_cld_altitude = 4000.0 @@ -1253,8 +1253,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,] + 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, 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_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,] < 274.0), False) # keep = np.where(keep, bt_11um[i,] < 275.0, False) @@ -1289,8 +1289,8 @@ 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 = 100.0 # meters - max_cld_depth = 1500.0 + closeness_top = 150.0 # meters + max_cld_depth = 700.0 max_altitude = 5000.0 max_cld_altitude = 5000.0 @@ -1332,8 +1332,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,] + 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, 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_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,] < 274.0), False) # keep = np.where(keep, bt_11um[i,] < 275.0, False)