From e7c875b0954f62f3cf162ed2f81430058bb0e5df Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Wed, 24 Mar 2021 20:40:30 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 8 ++++++--
1 file changed, 6 insertions(+), 2 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 686daf80..ba2b1d87 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -525,8 +525,9 @@ def run_qc(filename, filename_l1b):
phs = np.concatenate(phs)
opd = np.concatenate(opd)
- #print(np.histogram(bts, bins=20))
- #print(np.histogram(opd, bins=20))
+ print(np.histogram(bts, bins=20))
+ print(np.histogram(opd, bins=20))
+ print(np.histogram(phs, bins=8))
return mask
@@ -541,6 +542,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask,
cld_opd = cld_opd.reshape((num_obs, -1))
bt_11um = bt_11um.reshape((num_obs, -1))
+ skip = True
mask = []
for i in range(num_obs):
keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy
@@ -548,6 +550,8 @@ 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
+ if skip:
+ continue
keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False)
--
GitLab