From 03a3cc13e8920877af558fcb0ff5349bb293c40a Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Thu, 13 Apr 2023 13:20:24 -0500
Subject: [PATCH] snapshot...
---
modules/util/viirs_surfrad.py | 9 ++++++++-
1 file changed, 8 insertions(+), 1 deletion(-)
diff --git a/modules/util/viirs_surfrad.py b/modules/util/viirs_surfrad.py
index eed4fc92..5e6aa1c8 100644
--- a/modules/util/viirs_surfrad.py
+++ b/modules/util/viirs_surfrad.py
@@ -81,7 +81,7 @@ def process_cld_opd(grd_k):
keep_cld = np.where(keep, (0.1 < grd_k), False)
frac_cld = np.sum(keep_cld)/num_keep
# if not (0.50 < frac_cld < 0.85):
- if not (0.60 < frac_cld < 1.0):
+ if not (0.40 < frac_cld < 1.0):
return None
grd_k = np.where(np.invert(keep), 0, grd_k) # Convert NaN to 0
return grd_k
@@ -151,6 +151,7 @@ def run_all(directory, out_directory, day_night='ANY', pattern='clavrx_*.nc', st
h, b = np.histogram(valid_i.flatten(), bins=20, range=[0.0, 160.0])
hist_accum_valid_i += h
+
h, b = np.histogram(valid_m.flatten(), bins=20, range=[0.0, 160.0])
hist_accum_valid_m += h
@@ -164,6 +165,7 @@ def run_all(directory, out_directory, day_night='ANY', pattern='clavrx_*.nc', st
h, b = np.histogram(train_i.flatten(), bins=20, range=[0.0, 160.0])
hist_accum_train_i += h
+
h, b = np.histogram(train_m.flatten(), bins=20, range=[0.0, 160.0])
hist_accum_train_m += h
@@ -192,6 +194,7 @@ def run_all(directory, out_directory, day_night='ANY', pattern='clavrx_*.nc', st
h, b = np.histogram(valid_i.flatten(), bins=20, range=[0.0, 160.0])
hist_accum_valid_i += h
+
h, b = np.histogram(valid_m.flatten(), bins=20, range=[0.0, 160.0])
hist_accum_valid_m += h
@@ -205,6 +208,7 @@ def run_all(directory, out_directory, day_night='ANY', pattern='clavrx_*.nc', st
h, b = np.histogram(train_i.flatten(), bins=20, range=[0.0, 160.0])
hist_accum_train_i += h
+
h, b = np.histogram(train_m.flatten(), bins=20, range=[0.0, 160.0])
hist_accum_train_m += h
@@ -219,8 +223,11 @@ def run_all(directory, out_directory, day_night='ANY', pattern='clavrx_*.nc', st
print('** total_num_train_samples, total_num_valid_samples: ', total_num_train_samples, total_num_valid_samples)
print('--------------------------------------------------')
print(hist_accum_train_i)
+ print('----------')
print(hist_accum_train_m)
+ print('----------')
print(hist_accum_valid_i)
+ print('----------')
print(hist_accum_valid_m)
--
GitLab