From 6eef64c4bf923783019aeffb334bab668821adc9 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sat, 19 Mar 2022 13:36:47 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 31 ++++++++++++-------------------
1 file changed, 12 insertions(+), 19 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index ddf4d8ca..e0cdd8ae 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -943,29 +943,22 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask,
if num_keep == 0:
continue
- # Exp
+ # Test 1
keep = np.where(keep, (cld_top_hgt[i,] + closeness) > icing_alt[i], False)
- #if (np.sum(keep) / num_keep) < 0.50:
- # continue
-
- # Test1
- # keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False)
+ if (np.sum(keep) / num_keep) < 0.80:
+ continue
# 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)
+ 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)
- # Test5
- #keep = np.where(keep, (bt_11um[i,] < 270.0), 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)
- # Test6
- # keep = np.where(keep, (bt_11um[i,] > 228.0), False)
- # keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 270.0), False)
+ # Test5 and Test6
+ keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 270.0), False)
# Test3
# keep = np.where(keep, (cld_opd[i,] >= opd_thick_threshold) & (cld_phase[i,] == 4) & (cld_top_hgt[i,] > icing_alt[i]), False)
@@ -1021,7 +1014,7 @@ 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)
- #if (np.sum(keep) / num_keep) < 0.50:
+ #if (np.sum(keep) / num_keep) < 0.80:
# continue
# keep = np.where(keep, np.logical_and(bt_11um[i,] > 228.0, bt_11um[i,] < 273.0), False)
--
GitLab