From 637db97dd09363904946be8b92925a11bd3f6717 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 27 Mar 2023 16:47:48 -0500
Subject: [PATCH] snapshot...
---
modules/util/viirs_surfrad.py | 31 +++++++++++++++----------------
1 file changed, 15 insertions(+), 16 deletions(-)
diff --git a/modules/util/viirs_surfrad.py b/modules/util/viirs_surfrad.py
index db9aa86e..273d5d69 100644
--- a/modules/util/viirs_surfrad.py
+++ b/modules/util/viirs_surfrad.py
@@ -78,29 +78,28 @@ def process_cld_prob(grd_k):
return grd_k
-# def process_cld_opd(grd_k):
-# keep = np.invert(np.isnan(grd_k))
-# num_keep = np.sum(keep)
-# 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)
-# frac_keep = np.sum(keep)/num_keep
-# if frac_keep < 0.50:
-# return None
-# return grd_k
-
def process_cld_opd(grd_k):
keep = np.invert(np.isnan(grd_k))
num_keep = np.sum(keep)
- 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 > 4.0, grd_k < 40.0), False)
- frac_thin = np.sum(keep_thin) / num_keep
- frac_thik = np.sum(keep_thik) / num_keep
- if not (frac_thin >= 0.20 and frac_thik >= 0.20):
- # if not (frac_thin >= 0.20):
+ keep_cld = np.where(keep, np.logical_and(0.1 < grd_k, grd_k < 158.0), False)
+ frac_cld = np.sum(keep_cld)/num_keep
+ if 0.60 < frac_cld < 0.85:
return None
grd_k = np.where(np.invert(keep), 0, grd_k)
return grd_k
+# def process_cld_opd(grd_k):
+# keep = np.invert(np.isnan(grd_k))
+# num_keep = np.sum(keep)
+# 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 > 4.0, grd_k < 40.0), False)
+# frac_thin = np.sum(keep_thin) / num_keep
+# frac_thik = np.sum(keep_thik) / num_keep
+# if not (frac_thin >= 0.20 and frac_thik >= 0.20):
+# return None
+# grd_k = np.where(np.invert(keep), 0, grd_k)
+# return grd_k
+
def run_all(directory, out_directory, day_night='ANY', pattern='clavrx_*.nc', start=10):
cnt = start
--
GitLab