snapshot...

modules/icing/pirep_goes.py

@@ -717,7 +717,7 @@ def run_qc(filename, filename_l1b, day_night='ANY'):
     keep_idxs = []
 
     for i in range(len(mask)):
-        if (np.sum(mask[i]) / num_tested[i]) > 0.20:
+        if (np.sum(mask[i]) / 400) > 0.20:
             keep_idxs.append(idxs[i])
 
     print('num valid pireps: ', len(keep_idxs))
@@ -782,7 +782,9 @@ 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
-        num_tested.append(np.sum(keep))
+        num_keep = np.sum(keep)
+        if num_keep == 0:
+            continue
 
         keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False)
 
@@ -801,6 +803,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask,
 
         mask.append(keep)
         idxs.append(i)
+        num_tested.append(num_keep)
 
     return mask, idxs, num_tested