snapshot...

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))