From 623bfc3410839137ff58983a4fe41aa213cd2100 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Tue, 16 Mar 2021 09:52:05 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 14 +++++++-------
1 file changed, 7 insertions(+), 7 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 7ac608c6..9361e7ba 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -313,31 +313,31 @@ def run_qc(filename, filename_l1b):
f_l1b = h5py.File(filename_l1b, 'r')
bt_11um = f_l1b['temp_11_0um_nom'][:, 10:30, 10:30]
+ print('num pireps: ', len(icing_alt))
+
mask = apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, bt_11um)
f.close()
f_l1b.close()
- cnt = 0
bts = []
phs = []
opd = []
for i in range(len(mask)):
- if (np.sum(mask[i]) / 400) > 0.25:
- cnt += 1
+ if (np.sum(mask[i]) / 400) > 0.20:
bts.append((bt_11um[i,].flatten())[mask[i]])
phs.append((cld_phase[i,].flatten())[mask[i]])
opd.append((cld_opd[i,].flatten())[mask[i]])
#else:
# bts.append((bt_11um[i,].flatten())[mask[i]])
- print(cnt)
+ print('num valid pireps: ', len(bts))
bts = np.concatenate(bts)
phs = np.concatenate(phs)
opd = np.concatenate(opd)
- print(bts.shape)
- print(np.histogram(bts, bins=20))
- print(np.histogram(opd, bins=20))
+
+ #print(np.histogram(bts, bins=20))
+ #print(np.histogram(opd, bins=20))
return mask
--
GitLab