From 0eb6f6bc8b964b250aa4a2136b3dbd6f3aac9173 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Fri, 2 Dec 2022 11:37:36 -0600
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 12 +++++++-----
1 file changed, 7 insertions(+), 5 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index eb2ef35a..d02979ec 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1212,7 +1212,7 @@ 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
+ closeness_top = 100.0 # meters
max_depth = 3000.0
num_obs = len(icing_alt)
@@ -1251,9 +1251,10 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd
# Test 1
keep = np.where(keep, (cld_top_hgt[i,] + closeness_top) > icing_alt[i], False)
- keep = np.where(keep, (cld_top_hgt[i,] - max_depth) < icing_alt[i], False)
+ # keep = np.where(keep, (cld_top_hgt[i,] - max_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,] < 270.0), False)
+ keep = np.where(keep, np.logical_and(bt_11um[i,] < 274.0), False)
mask.append(keep)
idxs.append(i)
@@ -1309,9 +1310,10 @@ 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_depth) < icing_alt[i], False)
+ # keep = np.where(keep, (cld_top_hgt[i,] - max_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,] < 270.0), False)
+ keep = np.where(keep, np.logical_and(bt_11um[i,] < 274.0), False)
mask.append(keep)
idxs.append(i)
--
GitLab