From 3c0c00c527dec9fd51e43ca4ac6d089727b41dc4 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sun, 14 Mar 2021 15:13:15 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 13 ++++++-------
1 file changed, 6 insertions(+), 7 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index e159744b..c99e444c 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -294,15 +294,14 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, bt_11um, c
keep = keep_0 & keep_1 & keep_2 & keep_3
keep = np.where(keep, cld_top_hgt[:, ] > icing_alt[:], False)
- keep = np.invert(
- np.where(keep,
- np.logical_and(cld_phase[:, ] == 3,
+ keep = np.where(keep,
+ np.invert((cld_phase[:, ] == 3) &
np.logical_and(cld_top_hgt[:, ]+closeness > icing_alt[:], cld_top_hgt[:, ]-closeness < icing_alt[:])),
- False))
+ False)
keep = np.where(keep, (cld_opd[:, ] > opd_threshold) & (cld_phase[:, ] == 3) & (cld_top_hgt[:, ] > icing_alt[:]), False)
keep = np.where(keep, np.invert((cld_phase[:, ] == 3) & (cld_opd[:,] < 0.1) & (cld_top_hgt[:, ] > icing_alt[:]), False))
- keep = np.where(keep, np.invert(bt_11um[:,] > 270.0), False)
- keep = np.where(keep, np.invert(bt_11um[:,] < 228.0), False)
+ keep = np.where(keep, np.invert(bt_11um[:, ] > 270.0), False)
+ keep = np.where(keep, np.invert(bt_11um[:, ] < 228.0), False)
- return keep
\ No newline at end of file
+ return keep
--
GitLab