snapshot...

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