From 3c1f00304197a308739e82dcf9bc848bce5cf1d8 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Thu, 22 Dec 2022 14:55:01 -0600
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 20 ++++++++++----------
1 file changed, 10 insertions(+), 10 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 8602cf37..52f320f2 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1229,8 +1229,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
# if not (is_day(solzen[i,]) or is_night(solzen[i,])):
# continue
- if not (icing_alt[i] < max_altitude):
- continue
+ # if not (icing_alt[i] < max_altitude):
+ # continue
keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy
keep_1 = np.invert(np.isnan(cld_top_hgt[i,]))
@@ -1243,13 +1243,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_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,] < max_cld_altitude, 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, 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,] < 274.2, False)
# cld_hgt = cld_top_hgt[i, ].flatten()
# med_cld_hgt = np.median(cld_hgt[keep])
@@ -1307,8 +1307,8 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
# if not (is_day(solzen[i,]) or is_night(solzen[i,])):
# continue
- if not (icing_alt[i] < max_altitude):
- continue
+ # if not (icing_alt[i] < max_altitude):
+ # continue
keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy
keep_1 = np.invert(np.isnan(cld_top_hgt[i,]))
@@ -1321,13 +1321,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,] < max_cld_altitude, 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, 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,] < 274.2, False)
# cld_hgt = cld_top_hgt[i, ].flatten()
# med_cld_hgt = np.median(cld_hgt[keep])
--
GitLab