diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 2592ef432452a8de08e8faeebeeafd0d61fcfb5f..122fc9a1d99a3c9c309a8b742724ed53e1282116 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -673,17 +673,21 @@ def run_qc(filename, filename_l1b, day_night='ANY'):
     print(np.histogram(opd, bins=10))
     print(np.histogram(phs, bins=6))
 
-    bts = np.concatenate(bts_reject)
-    phs = np.concatenate(phs_reject)
-    opd = np.concatenate(opd_reject)
+    bts_r = np.concatenate(bts_reject)
+    phs_r = np.concatenate(phs_reject)
+    opd_r = np.concatenate(opd_reject)
 
-    print(np.histogram(bts, bins=10))
-    print(np.histogram(opd, bins=10))
-    print(np.histogram(phs, bins=6))
+    print(np.histogram(bts_r, bins=10))
+    print(np.histogram(opd_r, bins=10))
+    print(np.histogram(phs_r, bins=6))
+
+    return bts, bts_r, phs, phs_r, opd, opd_r
 
 
 def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, bt_11um, solzen, day_night='ANY'):
     opd_threshold = 2
+    if day_night == 'DAY':
+        opd_threshold = 5
     closeness = 100.0  # meters
     num_obs = len(icing_alt)
     cld_mask = cld_mask.reshape((num_obs, -1))