snapshot...

modules/util/viirs_l1b_l2.py

@@ -64,8 +64,9 @@ def process_cld_prob_(grd_k, dum):
     if num_keep / grd_k.size < 0.98:
         return None
     # hist_10 += np.histogram(grd_k.flatten(), range=[0.0, 1.0], bins=10)[0]
-    keep = np.where(keep, grd_k < 0.20, False)
-    frac_keep = np.sum(keep)/num_keep
+    keep_clr = np.where(keep, grd_k < 0.20, False)
+    # frac_keep = np.sum(keep)/num_keep
+    frac_keep = np.sum(keep_clr)/num_keep
     if not (0.40 < frac_keep < 0.60):
         return None
     grd_k = np.where(np.invert(keep), 0, grd_k)  # Convert NaN to 0