diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index f552edb2d6b5d57ab08e2e346977263545dad8f8..880e636f7564f72a1002d9a001b56ed5cf8336d1 100644
--- a/modules/icing/pirep_goes.py
+++ b/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))