From 78c72ca8e9a32428debc5065f1510ac65486c360 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Fri, 23 Apr 2021 13:01:27 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 7 +++++--
1 file changed, 5 insertions(+), 2 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 3335b159..ceed7a65 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -717,7 +717,7 @@ def run_qc(filename, filename_l1b, day_night='ANY'):
keep_idxs = []
for i in range(len(mask)):
- if (np.sum(mask[i]) / num_tested[i]) > 0.20:
+ if (np.sum(mask[i]) / 400) > 0.20:
keep_idxs.append(idxs[i])
print('num valid pireps: ', len(keep_idxs))
@@ -782,7 +782,9 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask,
keep_2 = np.invert(np.isnan(bt_11um[i,]))
keep_3 = np.invert(np.isnan(cld_opd[i,]))
keep = keep_0 & keep_1 & keep_2 & keep_3
- num_tested.append(np.sum(keep))
+ num_keep = np.sum(keep)
+ if num_keep == 0:
+ continue
keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False)
@@ -801,6 +803,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask,
mask.append(keep)
idxs.append(i)
+ num_tested.append(num_keep)
return mask, idxs, num_tested
--
GitLab