diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 0a9629f66b395ee4ac06fe7c4230aa3b2f8cb54c..7ac608c64e4c7057b76bd19fc13b580f2c92c497 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -321,20 +321,23 @@ def run_qc(filename, filename_l1b):
     cnt = 0
     bts = []
     phs = []
+    opd = []
     for i in range(len(mask)):
-        if (np.sum(mask[i]) / 400) > 0.15:
+        if (np.sum(mask[i]) / 400) > 0.25:
             cnt += 1
             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)
 
     bts = np.concatenate(bts)
     phs = np.concatenate(phs)
+    opd = np.concatenate(opd)
     print(bts.shape)
     print(np.histogram(bts, bins=20))
-    print(np.histogram(phs, bins=6))
+    print(np.histogram(opd, bins=20))
 
     return mask