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

snapshot...

parent a1ee7c2a
Branches
No related tags found
No related merge requests found
...@@ -92,11 +92,11 @@ def process_cld_opd(grd_k): ...@@ -92,11 +92,11 @@ 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_thin = np.where(keep, np.logical_and(grd_k > 0.1, grd_k < 4.0), False) keep_thin = np.where(keep, np.logical_and(grd_k > 0.1, grd_k < 4.0), False)
keep_thik = np.where(keep, np.logical_and(grd_k > 20.0, grd_k < 158.00), False) keep_thik = np.where(keep, np.logical_and(grd_k > 4.0, grd_k < 40.0), False)
frac_thin = np.sum(keep_thin) / num_keep frac_thin = np.sum(keep_thin) / num_keep
frac_thik = np.sum(keep_thik) / num_keep frac_thik = np.sum(keep_thik) / num_keep
# if not (frac_thin >= 0.10 and frac_thik >= 0.10): if not (frac_thin >= 0.20 and frac_thik >= 0.20):
if not (frac_thin >= 0.20): # if not (frac_thin >= 0.20):
return None return None
grd_k = np.where(np.invert(keep), 0, grd_k) grd_k = np.where(np.invert(keep), 0, grd_k)
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