snapshot...

modules/icing/pirep_goes.py

@@ -525,8 +525,9 @@ def run_qc(filename, filename_l1b):
     phs = np.concatenate(phs)
     opd = np.concatenate(opd)
 
-    #print(np.histogram(bts, bins=20))
-    #print(np.histogram(opd, bins=20))
+    print(np.histogram(bts, bins=20))
+    print(np.histogram(opd, bins=20))
+    print(np.histogram(phs, bins=8))
 
     return mask
 
@@ -541,6 +542,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask,
     cld_opd = cld_opd.reshape((num_obs, -1))
     bt_11um = bt_11um.reshape((num_obs, -1))
 
+    skip = True
     mask = []
     for i in range(num_obs):
         keep_0 = np.logical_or(cld_mask[i,] == 2, cld_mask[i,] == 3)  # cloudy
@@ -548,6 +550,8 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask,
         keep_2 = np.invert(np.isnan(bt_11um[i,]))
         keep_3 = np.invert(np.isnan(cld_opd[i,]))
         keep = keep_0 & keep_1 & keep_2 & keep_3
+        if skip:
+            continue
 
         keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False)