From 1f9cfce8fd0b1afcf5c15a13a7b1ab087e01ba60 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sat, 3 Feb 2024 16:16:36 -0600
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 13 ++++++++++---
1 file changed, 10 insertions(+), 3 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index a773a6f4..5bc168ea 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1397,8 +1397,8 @@ 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 = 200.0 # meters
- max_cld_depth = 1000.0
- max_altitude = 4000.0
+ max_cld_depth = 1200.0
+ max_altitude = 4500.0
max_cld_altitude = 4000.0
num_obs = len(icing_alt)
@@ -1433,8 +1433,15 @@ def apply_qc_icing_pireps_v3(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
+ # ctt_mean = np.nanmean(cld_top_temp[i,])
+ # if ctt_mean > 274.0:
+ # continue
+
+ # if icing_alt[i,] > max_altitude:
+ # continue
+
cth_mean = np.nanmean(cld_top_hgt[i,])
- if not (closeness_top+cth_mean > icing_alt[i,] > cth_mean - 1200.0):
+ if not (closeness_top+cth_mean > icing_alt[i,] > cth_mean - max_cld_depth):
continue
# keep = np.where(keep, np.logical_and(cld_top_hgt[i,] > icing_alt[i,],
--
GitLab