From f86f0631509b2c8b25b3978d3eb0be2f684baaf4 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sat, 17 Dec 2022 09:17:37 -0600
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 10 +++++-----
1 file changed, 5 insertions(+), 5 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index e637b0e4..9fd0ea3e 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1041,7 +1041,6 @@ def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20, icing
keep_idxs = []
for i in range(len(mask)):
- keep_idxs.append(idxs[i])
# frac = np.sum(mask[i]) / nx_x_ny
frac = np.sum(mask[i]) / num_tested[i]
if icing:
@@ -1049,6 +1048,7 @@ def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20, icing
keep_idxs.append(idxs[i])
elif frac > pass_thresh_frac:
keep_idxs.append(idxs[i])
+ # keep_idxs.append(idxs[i])
print('day_night, icing, all, valid, pass: ', day_night, icing, len(icing_alt), len(mask), len(keep_idxs))
print('-----------------------')
@@ -1277,8 +1277,8 @@ 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_cld_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,] < 275.0, False)
- # keep = np.where(keep, cld_top_temp[i,] < 275.0, False)
+ # keep = np.where(keep, bt_11um[i,] < 275.0, False)
+ keep = np.where(keep, cld_top_temp[i,] < 275.0, False)
# cld_hgt = cld_top_hgt[i, ].flatten()
# med_cld_hgt = np.median(cld_hgt[keep])
@@ -1355,8 +1355,8 @@ 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_cld_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,] < 275.0, False)
- # keep = np.where(keep, cld_top_temp[i,] < 275.0, False)
+ # keep = np.where(keep, bt_11um[i,] < 275.0, False)
+ keep = np.where(keep, cld_top_temp[i,] < 275.0, False)
# cld_hgt = cld_top_hgt[i, ].flatten()
# med_cld_hgt = np.median(cld_hgt[keep])
--
GitLab