snapshot...

modules/icing/pirep_goes.py

@@ -629,6 +629,8 @@ def run_qc(filename, filename_l1b, day_night='ANY'):
             keep_idxs.append(idxs[i])
         else:
             bts_reject.append((bt_11um[idxs[i],].flatten())[mask[i]])
+            phs_reject.append((cld_phase[idxs[i],].flatten())[mask[i]])
+            opd_reject.append((cld_opd[idxs[i],].flatten())[mask[i]])
 
     print('num valid pireps: ', len(keep_idxs))
     keep_idxs = np.array(keep_idxs)
@@ -671,6 +673,14 @@ 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)
+
+    print(np.histogram(bts, bins=10))
+    print(np.histogram(opd, bins=10))
+    print(np.histogram(phs, bins=6))
+
 
 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