From 0829cc8385beec8ce207ad06dc02970064a1ecfb Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 29 Apr 2024 13:06:06 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 12 ++++++------
1 file changed, 6 insertions(+), 6 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 78a53c6b..e1b3f497 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1464,9 +1464,9 @@ def apply_qc_icing_pireps_exp2(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cl
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 = 100.0 # meters
- max_cld_depth = 800.0
- max_altitude = 4800.0
- max_cld_altitude = 4800.0
+ max_cld_depth = 1500.0
+ max_altitude = 5000.0
+ max_cld_altitude = 5000.0
num_obs = len(icing_alt)
cld_mask = cld_mask.reshape((num_obs, -1))
@@ -1491,10 +1491,10 @@ def apply_qc_icing_pireps_v3(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
# 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(cld_geo_dz[i,]))
+ # keep_2 = np.invert(np.isnan(cld_geo_dz[i,]))
# keep_3 = np.invert(np.isnan(bt_11um[i,]))
# keep = keep_0 & keep_1 & keep_2 & keep_3
- keep = keep_0 & keep_1 & keep_2
+ keep = keep_0 & keep_1
num_keep = np.sum(keep)
if num_keep == 0:
continue
@@ -1515,7 +1515,7 @@ def apply_qc_icing_pireps_v3(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
keep = np.where(keep, cld_top_hgt[i,] < max_cld_altitude, False)
keep = np.where(keep, np.logical_and((cld_top_hgt[i,] + closeness_top) > icing_alt[i,],
- (cld_top_hgt[i,] - cld_geo_dz[i,]) < icing_alt[i,]), False)
+ (cld_top_hgt[i,] - max_cld_depth) < icing_alt[i,]), False)
# keep = np.where(keep, np.logical_and(cld_top_temp[i,] > 228.0, cld_top_temp[i,] < 274.0), False)
--
GitLab