snapshot...

modules/icing/pirep_goes.py

@@ -616,12 +616,14 @@ def run_qc(filename, filename_l1b):
     phs = []
     opd = []
     opd_dc = []
+    keep_idxs =[]
     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]]
+            keep_idxs.append(i)
         #else:
         #    bts.append((bt_11um[i,].flatten())[mask[i]])
     print('num valid pireps: ', len(bts))