From 8c39ebbe8b9f77459615dba79498b4a0ce863b4f Mon Sep 17 00:00:00 2001 From: tomrink <rink@ssec.wisc.edu> Date: Thu, 17 Mar 2022 14:28:05 -0500 Subject: [PATCH] snapshot... --- modules/icing/pirep_goes.py | 8 +++----- 1 file changed, 3 insertions(+), 5 deletions(-) diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index eee6d50f..b56cba31 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -961,11 +961,11 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, # keep = np.where(keep, np.invert((cld_phase[i,] == 4) & (cld_opd[i,] < opd_thin_threshold) & (cld_top_hgt[i,] > icing_alt[i])), False) # Test5 - keep = np.where(keep, (bt_11um[i,] < 270.0), False) + #keep = np.where(keep, (bt_11um[i,] < 270.0), False) # 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) @@ -1024,8 +1024,7 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mas if (np.sum(keep) / num_keep) < 0.60: continue - # keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 273.0), False) - keep = np.where(keep, 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) @@ -1726,7 +1725,6 @@ def split_data(times): [get_timestamp('2022-01-01_00:00'), get_timestamp('2022-01-07_23:59')], [get_timestamp('2022-03-01_00:00'), get_timestamp('2022-03-07_23:59')]] - test_time_idxs = [] for t_rng in time_ranges: tidxs = np.searchsorted(times, t_rng) -- GitLab