snapshot...

modules/util/viirs_l1b_l2.py

@@ -64,8 +64,10 @@ 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, np.logical_and(0.10 < grd_k, grd_k < 0.90), False)
-    if np.sum(keep)/num_keep < 0.30:
+    # keep = np.where(keep, np.logical_and(0.10 < grd_k, grd_k < 0.90), False)
+    keep = np.where(keep, grd_k < 0.50, False)
+    frac_keep = np.sum(keep)/num_keep
+    if 0.30 < frac_keep < 0.70:
         return None
     grd_k = np.where(np.invert(keep), 0, grd_k)
     return grd_k
@@ -89,7 +91,8 @@ def process_cld_opd_(grd_k, dum):
         return None
     grd_k = np.where(np.invert(keep), 0, grd_k)
     keep = np.where(keep, np.logical_and(0.1 < grd_k, grd_k < 158.0), False)
-    if np.sum(keep)/num_keep < 0.50:
+    frac_keep = np.sum(keep)/num_keep
+    if frac_keep < 0.50:
         return None
     return grd_k