diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index afc286e68080d6176eae30f257901762a9aad298..7e55945892fad3eb9515e399f737f785886b5b9a 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -1045,7 +1045,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 = 50.0 # meters + closeness_top = 100.0 # meters + closeness_bot = 200.0 # meters num_obs = len(icing_alt) cld_mask = cld_mask.reshape((num_obs, -1)) cld_top_hgt = cld_top_hgt.reshape((num_obs, -1)) @@ -1077,8 +1078,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd continue # Test 1 - 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) # # Test2 # keep = np.where(keep,