From d21a604cbc5e86cd7263590c57962282594f9fb5 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 24 Apr 2023 15:56:26 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 24 +++++++++++++-----------
1 file changed, 13 insertions(+), 11 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 8d963234..06d0841d 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1210,8 +1210,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
opd_thick_threshold = 2
opd_thin_threshold = 0.1
- closeness_top = 200.0 # meters
- max_cld_depth = 1500.0
+ closeness_top = 100.0 # meters
+ max_cld_depth = 500.0
max_altitude = 4000.0
max_cld_altitude = 4000.0
@@ -1252,12 +1252,13 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
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,] - max_cld_depth) < 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, 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, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 274.0), False)
# keep = np.where(keep, bt_11um[i,] < 275.0, False)
- keep = np.where(keep, cld_top_temp[i,] < 274.2, 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])
@@ -1288,7 +1289,7 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
opd_thick_threshold = 2
opd_thin_threshold = 0.1
- closeness_top = 200.0 # meters
+ closeness_top = 100.0 # meters
max_cld_depth = 1500.0
max_altitude = 5000.0
max_cld_altitude = 5000.0
@@ -1330,12 +1331,13 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_
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,] - max_cld_depth) < 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, 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, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 274.0), False)
# keep = np.where(keep, bt_11um[i,] < 275.0, False)
- keep = np.where(keep, cld_top_temp[i,] < 274.2, 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