diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 5d1eb8095bd2ce7ef91c4d48c832b22edeee0e1f..9d966cd5c4d7853484b056c87f372b75ac3667c8 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -1237,8 +1237,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd # continue # if not (500.0 < icing_alt[i] < 3000.0): - # if not (icing_alt[i] < 3000.0): - # continue + if not (icing_alt[i] < 4000.0): + continue keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy keep_1 = np.invert(np.isnan(cld_top_hgt[i,])) @@ -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, bt_11um[i,] < 274.0, False) + keep = np.where(keep, bt_11um[i,] < 275.0, False) mask.append(keep) idxs.append(i) @@ -1297,8 +1297,8 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_ # continue # if not (500.0 < icing_alt[i] < 3000.0): - # if not (icing_alt[i] < 3000.0): - # continue + if not (icing_alt[i] < 4000.0): + continue keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy keep_1 = np.invert(np.isnan(cld_top_hgt[i,])) @@ -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, bt_11um[i,] < 274.0, False) + keep = np.where(keep, bt_11um[i,] < 275.0, False) mask.append(keep) idxs.append(i)