From 4b05cd1247d985f4a3019c0beb8c230ed7bae351 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Fri, 14 Jan 2022 12:07:00 -0600
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 18 ++++++++++++++++++
1 file changed, 18 insertions(+)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 08f1e15a..4f35c749 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -873,6 +873,21 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask,
# Exp
# keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False)
+ # Test1
+ keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False)
+
+ # Test2
+ keep = np.where(keep,
+ np.invert((cld_phase[i,] == 4) &
+ np.logical_and(cld_top_hgt[i,]+closeness > icing_alt[i], cld_top_hgt[i,]-closeness < icing_alt[i])),
+ False)
+
+ # Test4
+ keep = np.where(keep, np.invert((cld_phase[i,] == 4) & (cld_opd[i,] < opd_thin_threshold) & (cld_top_hgt[i,] > icing_alt[i])), False)
+
+ # Test5
+ keep = np.where(keep, np.invert(bt_11um[i,] > 270.0), False)
+
# Test6
keep = np.where(keep, np.invert(bt_11um[i,] < 228.0), False)
@@ -929,6 +944,9 @@ def apply_qc_no_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mas
# Exp
# keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False)
+ # Test1
+ keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False)
+
# Test3
keep = np.where(keep, np.invert((cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4) & (cld_top_hgt[i,] > icing_alt[i])), False)
--
GitLab