From 7fbaee29fef1668041fabcc5f99b5de8503f7ffc Mon Sep 17 00:00:00 2001 From: tomrink <rink@ssec.wisc.edu> Date: Fri, 18 Mar 2022 11:47:00 -0500 Subject: [PATCH] snapshot... --- modules/icing/pirep_goes.py | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index b56cba31..ddf4d8ca 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -945,8 +945,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, # Exp keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False) - if (np.sum(keep) / num_keep) < 0.60: - continue + #if (np.sum(keep) / num_keep) < 0.50: + # continue # Test1 # keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False) @@ -965,7 +965,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, # Test6 # keep = np.where(keep, (bt_11um[i,] > 228.0), 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,] > 228.0, bt_11um[i,] < 270.0), False) # Test3 # keep = np.where(keep, (cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4) & (cld_top_hgt[i,] > icing_alt[i]), False) @@ -1021,10 +1021,10 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mas # Exp keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False) - if (np.sum(keep) / num_keep) < 0.60: - continue + #if (np.sum(keep) / num_keep) < 0.50: + # continue - keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 273.0), False) + # keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 273.0), False) # Test6 # keep = np.where(keep, (bt_11um[i,] < 228.0), False) -- GitLab