From 788cad9c68c534667f319234dabf4373a55cb04f Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 19 Dec 2022 13:44:56 -0600
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 18 +++++++++---------
1 file changed, 9 insertions(+), 9 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 876207d6..02611741 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1012,7 +1012,7 @@ def run_daynight(filename, filename_l1b, day_night='ANY'):
f_l1b.close()
-def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20, icing=True):
+def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20, cloudy_fraction=0.5, icing=True):
f = h5py.File(filename, 'r')
icing_alt = f['icing_altitude'][:]
@@ -1034,9 +1034,9 @@ def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20, icing
bt_11um = f_l1b['temp_11_0um_nom'][:, y_a:y_b, x_a:x_b]
if icing:
- mask, idxs, num_tested = apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd, cld_mask, bt_11um, sol_zen, sat_zen, cld_temp, day_night=day_night)
+ mask, idxs, num_tested = apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd, cld_mask, bt_11um, sol_zen, sat_zen, cld_temp, day_night=day_night, cloud_fraction=cloudy_fraction)
else:
- mask, idxs, num_tested = apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd, cld_mask, bt_11um, sol_zen, sat_zen, cld_temp, day_night=day_night)
+ mask, idxs, num_tested = apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd, cld_mask, bt_11um, sol_zen, sat_zen, cld_temp, day_night=day_night, cloudy_frac=cloudy_fraction)
keep_idxs = []
@@ -1272,13 +1272,13 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
if num_keep / nx_x_ny < cloud_fraction: # At least this fraction cloudy
continue
- keep = np.where(keep, (cld_top_hgt[i,] + closeness_top) > icing_alt[i], False)
+ ##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) < icing_alt[i], 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,] < 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, cld_top_temp[i,] < 275.0, False)
# cld_hgt = cld_top_hgt[i, ].flatten()
# med_cld_hgt = np.median(cld_hgt[keep])
@@ -1350,13 +1350,13 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
if num_keep / nx_x_ny < cloudy_frac: # At least this fraction cloudy
continue
- keep = np.where(keep, (cld_top_hgt[i,] + closeness_top) > icing_alt[i], False)
+ ##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) < icing_alt[i], 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,] < 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, cld_top_temp[i,] < 275.0, False)
# cld_hgt = cld_top_hgt[i, ].flatten()
# med_cld_hgt = np.median(cld_hgt[keep])
--
GitLab