snapshot...

modules/icing/pirep_goes.py

@@ -710,14 +710,14 @@ def run_qc(filename, filename_l1b, day_night='ANY'):
 
     print('num pireps all: ', len(icing_alt))
 
-    mask, idxs = apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, bt_11um, sol_zen, day_night=day_night)
+    mask, idxs, num_tested = apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask, bt_11um, sol_zen, day_night=day_night)
 
-    print('num pireps day_night: ', len(mask), day_night)
+    print('num pireps, day_night: ', len(mask), day_night)
 
     keep_idxs = []
 
     for i in range(len(mask)):
-        if (np.sum(mask[i]) / 400) > 0.20:
+        if (np.sum(mask[i]) / num_tested[i]) > 0.20:
             keep_idxs.append(idxs[i])
 
     print('num valid pireps: ', len(keep_idxs))
@@ -768,6 +768,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask,
 
     mask = []
     idxs = []
+    num_tested = []
     for i in range(num_obs):
         if day_night == 'NIGHT':
             if is_day(solzen[i,]):
@@ -781,6 +782,7 @@ 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))
 
         keep = np.where(keep, cld_top_hgt[i,] > icing_alt[i], False)
 
@@ -800,7 +802,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask,
         mask.append(keep)
         idxs.append(i)
 
-    return mask, idxs
+    return mask, idxs, num_tested
 
 
 def fov_extract(outfile='/home/rink/fovs_l1b_out.h5', train_params=l1b_ds_list):