From 6ea461f22015dcbd73a62a862fb278a2df4c8988 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 10 Jul 2023 15:47:57 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 9 ++++-----
1 file changed, 4 insertions(+), 5 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 0d14f719..68ca3cc7 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1025,7 +1025,6 @@ def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20, cloud
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('-----------------------')
@@ -1251,8 +1250,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,] + closeness_top) > icing_alt[i], False)
# keep = np.where(keep, cld_top_hgt[i,] < max_cld_altitude, False)
- ## keep = np.where(keep, cld_top_hgt[i,] > max_cld_depth, False)
- ## keep = np.where(keep, (cld_top_hgt[i,] - max_cld_depth) < icing_alt[i], False)
+ # keep = np.where(keep, cld_top_hgt[i,] > max_cld_depth, False)
+ # 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,] < 274.0), False)
# keep = np.where(keep, bt_11um[i,] < 275.0, False)
@@ -1330,8 +1329,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,] + closeness_top) > icing_alt[i], False)
# keep = np.where(keep, cld_top_hgt[i,] < max_cld_altitude, False)
- ## keep = np.where(keep, cld_top_hgt[i,] > max_cld_depth, False)
- ## keep = np.where(keep, (cld_top_hgt[i,] - max_cld_depth) < icing_alt[i], False)
+ # keep = np.where(keep, cld_top_hgt[i,] > max_cld_depth, False)
+ # 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,] < 274.0), False)
# keep = np.where(keep, bt_11um[i,] < 275.0, False)
--
GitLab