From cdad9919398c2bbdf4c2e526a3446758b40187ce Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Fri, 26 Apr 2024 15:15:40 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 9 +++++----
1 file changed, 5 insertions(+), 4 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 3e508d63..d2028a9b 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1466,7 +1466,7 @@ 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 = 50.0 # meters
- max_cld_depth = 800.0
+ max_cld_depth = 500.0
max_altitude = 4500.0
max_cld_altitude = 4000.0
@@ -1493,9 +1493,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(bt_11um[i,]))
- keep_3 = np.invert(np.isnan(cld_geo_dz[i,]))
- keep = keep_0 & keep_1 & keep_2 & keep_3
+ # keep_2 = np.invert(np.isnan(bt_11um[i,]))
+ # keep_3 = np.invert(np.isnan(cld_geo_dz[i,]))
+ # keep = keep_0 & keep_1 & keep_2 & keep_3
+ keep = keep_0 & keep_1
num_keep = np.sum(keep)
if num_keep == 0:
continue
--
GitLab