From f7f37676c8e580fd47b8c9e352f6952e3b2ac6da Mon Sep 17 00:00:00 2001 From: tomrink <rink@ssec.wisc.edu> Date: Thu, 17 Mar 2022 11:26:58 -0500 Subject: [PATCH] snapshot... --- modules/icing/pirep_goes.py | 14 ++++++++------ 1 file changed, 8 insertions(+), 6 deletions(-) diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 31a19598..4a5f07ac 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -945,7 +945,7 @@ 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.50: + if (np.sum(keep) / num_keep) < 0.60: continue # Test1 @@ -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) @@ -1021,10 +1021,11 @@ 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.50: + 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, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 273.0), False) + keep = np.where(keep, bt_11um[i,] < 273.0, False) # Test6 # keep = np.where(keep, (bt_11um[i,] < 228.0), False) @@ -1719,7 +1720,8 @@ def split_data(times): [get_timestamp('2019-05-01_00:00'), get_timestamp('2019-05-07_23:59')], [get_timestamp('2019-07-01_00:00'), get_timestamp('2019-07-07_23:59')], [get_timestamp('2019-09-01_00:00'), get_timestamp('2019-09-07_23:59')], - [get_timestamp('2019-11-01_00:00'), get_timestamp('2019-11-07_23:59')]] + [get_timestamp('2019-11-01_00:00'), get_timestamp('2019-11-07_23:59')], + [get_timestamp('2021-09-24_00:00'), get_timestamp('2021-10-01_23:59')]] test_time_idxs = [] -- GitLab