diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index c464a4aae98dd6c98a90fa9ee1e1cad6c919e0b2..21f24ff1e15aed44e55637b1b85f2f1ebd34abf4 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1084,7 +1084,7 @@ def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20, cloud
         mask, idxs, num_tested = apply_qc_icing_pireps_v3(icing_alt, cld_top_hgt, cld_geo_dz, cld_phase, cld_opd, cld_mask, bt_11um, sol_zen, sat_zen, cld_temp, day_night=day_night, cloudy_frac=cloudy_frac)
 
     keep_idxs = []
-    mask_keep_idxs = []
+    keep_mask = []
 
     for i in range(len(mask)):
         # keep_idxs.append(idxs[i])
@@ -1092,10 +1092,10 @@ def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20, cloud
         if icing:
             if frac > pass_thresh_frac:
                 keep_idxs.append(idxs[i])
-                mask_keep_idxs.append(i)
+                keep_mask.append(mask)
         elif frac > pass_thresh_frac:
             keep_idxs.append(idxs[i])
-            mask_keep_idxs.append(i)
+            keep_mask.append(mask)
 
     if len(keep_idxs) == 0:
         f.close()
@@ -1105,9 +1105,7 @@ def run_qc(filename, filename_l1b, day_night='ANY', pass_thresh_frac=0.20, cloud
     print('day_night, icing, all, valid, pass: ', day_night, icing, len(icing_alt), len(mask), len(keep_idxs))
     print('-----------------------')
     keep_idxs = np.array(keep_idxs)
-    mask_keep_idxs = np.array(mask_keep_idxs)
     mask = np.concatenate(mask).reshape((-1, ny, nx))
-    mask = mask[mask_keep_idxs,]
 
     data_dct = {}
     for didx, ds_name in enumerate(ds_list):