From c795468b7311e0a916bfb9d3e958b2dd824d5eb0 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 20 Dec 2021 14:49:38 -0600
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 14 ++++++++++----
1 file changed, 10 insertions(+), 4 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 4baa96e4..0c157f20 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -2051,7 +2051,13 @@ def run_icing_predict(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', output_dir=h
data_dct, ll, cc = make_for_full_domain_predict(h5f, name_list=train_params, satellite=satellite, domain=domain)
# ancil_data_dct, _, _ = make_for_full_domain_predict(h5f, name_list=['cld_height_acha', 'cld_geo_thick'])
+
solzen, satzen = make_for_full_domain_predict2(h5f, satellite=satellite, domain=domain)
+ keep = np.logical_or(lats_2d > -63.0, lats_2d < 63.0)
+ keep = np.where(keep, satzen < 70, False)
+ if day_night == 'DAY':
+ keep = np.where(keep, solzen < 80, False)
+
if fidx == 0:
num_elems = len(cc)
num_lines = len(ll)
@@ -2065,10 +2071,10 @@ def run_icing_predict(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', output_dir=h
for flvl in flt_lvls:
probs = probs_2d_dct[flvl]
preds = preds_2d_dct[flvl]
- keep = np.logical_or(lats_2d > -63.0, lats_2d < 63.0)
- keep = np.where(keep, satzen < 70, False)
- if day_night == 'DAY':
- keep = np.where(keep, solzen < 80, False)
+ # keep = np.logical_or(lats_2d > -63.0, lats_2d < 63.0)
+ # keep = np.where(keep, satzen < 70, False)
+ # if day_night == 'DAY':
+ # keep = np.where(keep, solzen < 80, False)
np.where(keep, preds, -1)
np.where(keep, probs, -1.0)
--
GitLab