Skip to content
Snippets Groups Projects
Commit fc2cd71f authored by tomrink's avatar tomrink
Browse files

snapshot...

parent 8e8a119f
Branches
No related tags found
No related merge requests found
...@@ -79,9 +79,11 @@ def process_cld_prob(grd_k): ...@@ -79,9 +79,11 @@ def process_cld_prob(grd_k):
def process_cld_opd(grd_k): def process_cld_opd(grd_k):
keep = np.invert(np.isnan(grd_k)) keep = np.invert(np.isnan(grd_k))
num_keep = np.sum(keep) num_keep = np.sum(keep)
keep_cld = np.where(keep, np.logical_and(0.1 < grd_k, grd_k < 158.0), False) # keep_cld = np.where(keep, np.logical_and(0.1 < grd_k, grd_k < 158.0), False)
keep_cld = np.where(keep, (0.1 < grd_k), False)
frac_cld = np.sum(keep_cld)/num_keep frac_cld = np.sum(keep_cld)/num_keep
if not (0.10 < frac_cld < 0.90): # if not (0.50 < frac_cld < 0.85):
if not (0.50 < frac_cld < 1.0):
return None return None
grd_k = np.where(np.invert(keep), 0, grd_k) # Convert NaN to 0 grd_k = np.where(np.invert(keep), 0, grd_k) # Convert NaN to 0
return grd_k return grd_k
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment