From 52752273936c2325ed9ab8d44f7acbf9036bf4be Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Fri, 24 Nov 2023 11:27:46 -0600
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 16 ++++++++--------
1 file changed, 8 insertions(+), 8 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index e22d28ad..ff23e752 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1399,9 +1399,9 @@ def apply_qc_icing_pireps_exp2(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cl
return mask, idxs, num_tested
-def apply_qc_icing_pireps_new(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd, cld_mask, bt_11um, solzen, satzen, cld_top_temp, day_night='ANY', cloudy_frac=0.5):
+def apply_qc_icing_pireps_v3(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd, cld_mask, bt_11um, solzen, satzen, cld_top_temp, day_night='ANY', cloudy_frac=0.5):
- closeness_top = 200.0 # meters
+ closeness_top = 100.0 # meters
max_cld_depth = 1000.0
max_altitude = 4000.0
max_cld_altitude = 4000.0
@@ -1426,8 +1426,8 @@ def apply_qc_icing_pireps_new(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld
elif day_night == 'ANY':
pass
- # keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy
- keep_0 = cld_mask[i, ] == 3 # Confident cloudy
+ # keep_0 = cld_mask[i, ] == 3 # Confident cloudy
+ keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy
keep_1 = np.invert(np.isnan(cld_top_hgt[i,]))
keep_2 = np.invert(np.isnan(bt_11um[i,]))
keep_3 = np.invert(np.isnan(cld_geo_dz[i,]))
@@ -1439,12 +1439,12 @@ def apply_qc_icing_pireps_new(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld
continue
cth_mean = np.nanmean(cld_top_hgt[i,])
- if not (closeness_top+cth_mean > icing_alt[i,] > cth_mean - 1000.0):
+ if not (closeness_top+cth_mean > icing_alt[i,] > cth_mean - 1200.0):
continue
- keep = np.where(keep, np.logical_and(cld_top_hgt[i,] > icing_alt[i,],
- (cld_top_hgt[i,] - 3000.0) < icing_alt[i,]), False)
- keep = np.where(keep, np.logical_and(cld_top_temp[i,] > 228.0, cld_top_temp[i,] < 274.0), False)
+ # keep = np.where(keep, np.logical_and(cld_top_hgt[i,] > icing_alt[i,],
+ # (cld_top_hgt[i,] - 3000.0) < icing_alt[i,]), False)
+ # keep = np.where(keep, np.logical_and(cld_top_temp[i,] > 228.0, cld_top_temp[i,] < 274.0), False)
mask.append(keep)
idxs.append(i)
--
GitLab