From 9e90b054ce9390b68596835adec80381e32efbaf Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 29 Mar 2021 08:44:39 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 24 ++++++++++++++++++++++--
1 file changed, 22 insertions(+), 2 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 880e636f..e9516e6a 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -593,7 +593,7 @@ 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):
+def run_qc(filename, filename_l1b, outfile, outfile_l1b):
f = h5py.File(filename, 'r')
icing_alt = f['icing_altitude'][:]
cld_top_hgt = f['cld_height_acha'][:, 10:30, 10:30]
@@ -638,6 +638,26 @@ def run_qc(filename, filename_l1b):
#print(np.histogram(opd_dc, bins=10))
print(np.histogram(phs, bins=6))
+ keep_idxs = np.array(keep_idxs)
+
+ data_dct = {}
+ for didx, ds_name in enumerate(ds_list):
+ data_dct[ds_name] = f[ds_name][keep_idxs,]
+ lon_c = f['longitude'][keep_idxs]
+ lat_c = f['latitude'][keep_idxs]
+ time_s = f['time'][keep_idxs]
+ fl_alt_s = f['icing_altidue'][keep_idxs]
+ ice_int_s = f['icing_intensity'][keep_idxs]
+ unq_ids = f['unique_ids'][keep_idxs]
+
+ create_file(outfile, data_dct, ds_list, ds_types, lon_c, lat_c, time_s, fl_alt_s, ice_int_s, unq_ids)
+
+ data_dct = {}
+ for didx, ds_name in enumerate(l1b_ds_list):
+ data_dct[ds_name] = f[ds_name][keep_idxs]
+
+ create_file(outfile_l1b, data_dct, l1b_ds_list, l1b_ds_types, lon_c, lat_c, time_s, fl_alt_s, ice_int_s, unq_ids)
+
return mask
@@ -651,7 +671,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
+ skip = False
mask = []
for i in range(num_obs):
keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3) # cloudy
--
GitLab