From ca133b2f87aeee9872e69f1d5b4869830ed070f5 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Thu, 23 Mar 2023 09:28:52 -0500
Subject: [PATCH] snapshot...
---
modules/util/viirs_surfrad.py | 24 ++++++++++++++++++------
1 file changed, 18 insertions(+), 6 deletions(-)
diff --git a/modules/util/viirs_surfrad.py b/modules/util/viirs_surfrad.py
index 6a83650f..ad03cd7a 100644
--- a/modules/util/viirs_surfrad.py
+++ b/modules/util/viirs_surfrad.py
@@ -29,8 +29,8 @@ keep_out_opd = ['/ships19/cloud/scratch/cphillips/clavrx/run_viirs_superres/site
keep_out = keep_out_opd
-target_param = 'cloud_probability'
-# target_param = 'cld_opd_dcomp'
+# target_param = 'cloud_probability'
+target_param = 'cld_opd_dcomp'
group_name_i = 'super/'
group_name_m = 'orig/'
@@ -78,14 +78,26 @@ 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)
- 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:
+ keep_vthin = np.where(keep, grd_k < 0.01, False)
+ keep_vthik = np.where(keep, grd_k < 150.0, False)
+ frac_vthin = np.sum(keep_vthin) / num_keep
+ frac_vthik = np.sum(keep_vthik) / num_keep
+ if not (frac_vthin >= 0.10 and frac_vthik >= 0.10):
return None
+ grd_k = np.where(np.invert(keep), 0, grd_k)
return grd_k
--
GitLab