diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 17261d03d17871bda33bb4ca9fadc67aa2c6934a..84c3c826876fd27a94dddbd73f1e90adcfefe706 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -1111,7 +1111,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 = 50.0 # meters + closeness_top = 100.0 # meters + closeness_bot = 200.0 num_obs = len(icing_alt) cld_mask = cld_mask.reshape((num_obs, -1)) cld_top_hgt = cld_top_hgt.reshape((num_obs, -1)) @@ -1142,8 +1143,8 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_ if num_keep == 0: continue - keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False) - keep = np.where(keep, (cld_top_hgt[i,] - cld_geo_dz[i,] - closeness) < icing_alt[i], False) + keep = np.where(keep, (cld_top_hgt[i,] + closeness_top) > icing_alt[i], False) + keep = np.where(keep, (cld_top_hgt[i,] - cld_geo_dz[i,] - closeness_bot) < icing_alt[i], False) # keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 273.0), False)