From 41a3af9423b753c01bc6bdd404ca694e592df77a Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Tue, 30 Mar 2021 16:29:32 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 12 +++++-------
1 file changed, 5 insertions(+), 7 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 8864bd4f..0ac29ffd 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -594,13 +594,15 @@ def process_1(ice_dct, no_ice_dct, neg_ice_dct):
return new_ice_dct, new_no_ice_dct, new_neg_ice_dct
-def run_qc(filename, filename_l1b, outfile, outfile_l1b, day_night='ANY'):
+def run_qc(filename, filename_l1b, day_night='ANY'):
f = h5py.File(filename, 'r')
icing_alt = f['icing_altitude'][:]
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]
+ if day_night == 'DAY':
+ cld_opd = f['cld_opd_dcomp'][:, 10:30, 10:30]
+ else:
+ cld_opd = f['cld_opd_acha'][:, 10:30, 10:30]
cld_mask = f['cloud_mask'][:, 10:30, 10:30]
sol_zen = f['solar_zenith_angle'][:, 10:30, 10:30]
@@ -614,7 +616,6 @@ def run_qc(filename, filename_l1b, outfile, outfile_l1b, day_night='ANY'):
bts_pass = []
phs_pass = []
opd_pass = []
- opd_dc = []
keep_idxs = []
bts_reject = []
phs_reject = []
@@ -625,7 +626,6 @@ def run_qc(filename, filename_l1b, outfile, outfile_l1b, day_night='ANY'):
bts_pass.append((bt_11um[idxs[i],].flatten())[mask[i]])
phs_pass.append((cld_phase[idxs[i],].flatten())[mask[i]])
opd_pass.append((cld_opd[idxs[i],].flatten())[mask[i]])
- #opd_dc.append(cld_opd_dc[i,].flatten())[mask[i]]
keep_idxs.append(idxs[i])
else:
bts_reject.append((bt_11um[idxs[i],].flatten())[mask[i]])
@@ -666,11 +666,9 @@ def run_qc(filename, filename_l1b, outfile, outfile_l1b, day_night='ANY'):
bts = np.concatenate(bts_pass)
phs = np.concatenate(phs_pass)
opd = np.concatenate(opd_pass)
- #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))
--
GitLab