diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index d02979ecd44eb754aa0c573ba5b9cd9e9dad355f..5d1eb8095bd2ce7ef91c4d48c832b22edeee0e1f 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -1254,7 +1254,7 @@ 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_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,] < 274.0), False) + keep = np.where(keep, bt_11um[i,] < 274.0, False) mask.append(keep) idxs.append(i) @@ -1272,7 +1272,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_depth = 3000.0 num_obs = len(icing_alt) @@ -1313,7 +1313,7 @@ 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_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,] < 274.0), False) + keep = np.where(keep, bt_11um[i,] < 274.0, False) mask.append(keep) idxs.append(i)