diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 7ac608c64e4c7057b76bd19fc13b580f2c92c497..9361e7ba1817f2df2abfb7d5dc0899f0db5367c8 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -313,31 +313,31 @@ def run_qc(filename, filename_l1b):
     f_l1b = h5py.File(filename_l1b, 'r')
     bt_11um = f_l1b['temp_11_0um_nom'][:, 10:30, 10:30]
 
+    print('num pireps: ', len(icing_alt))
+
     mask = apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, bt_11um)
 
     f.close()
     f_l1b.close()
 
-    cnt = 0
     bts = []
     phs = []
     opd = []
     for i in range(len(mask)):
-        if (np.sum(mask[i]) / 400) > 0.25:
-            cnt += 1
+        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]])
         #else:
         #    bts.append((bt_11um[i,].flatten())[mask[i]])
-    print(cnt)
+    print('num valid pireps: ', len(bts))
 
     bts = np.concatenate(bts)
     phs = np.concatenate(phs)
     opd = np.concatenate(opd)
-    print(bts.shape)
-    print(np.histogram(bts, bins=20))
-    print(np.histogram(opd, bins=20))
+
+    #print(np.histogram(bts, bins=20))
+    #print(np.histogram(opd, bins=20))
 
     return mask