From 8c978a180dc87a8397e96be25a375cf1e45274c3 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Thu, 25 Mar 2021 10:30:15 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 5 +++++
1 file changed, 5 insertions(+)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 1cfdf1ed..ff842153 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -497,6 +497,7 @@ def run_qc(filename, filename_l1b):
cld_top_hgt = f['cld_height_acha'][:, 10:30, 10:30]
cld_phase = f['cloud_phase'][:, 10:30, 10:30]
cld_opd = f['cld_opd_acha'][:, 10:30, 10:30]
+ cld_opd_dc = f['cld_opd_dcomp'][:, 10:30, 10:30]
cld_mask = f['cloud_mask'][:, 10:30, 10:30]
f_l1b = h5py.File(filename_l1b, 'r')
@@ -512,11 +513,13 @@ def run_qc(filename, filename_l1b):
bts = []
phs = []
opd = []
+ opd_dc = []
for i in range(len(mask)):
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]])
+ opd_dc.append(cld_opd_dc[i,].flatten())[mask[i]]
#else:
# bts.append((bt_11um[i,].flatten())[mask[i]])
print('num valid pireps: ', len(bts))
@@ -524,9 +527,11 @@ def run_qc(filename, filename_l1b):
bts = np.concatenate(bts)
phs = np.concatenate(phs)
opd = np.concatenate(opd)
+ opd_dc = np.concatenate(opd_dc)
print(np.histogram(bts, bins=10))
print(np.histogram(opd, bins=10))
+ print(np.histogram(opd_dc, bins=10))
print(np.histogram(phs, bins=6))
return mask
--
GitLab