From 709d5c46602577642748d18e336ef9d4e10c5be8 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Thu, 27 Jan 2022 15:37:19 -0600
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 8 +++-----
1 file changed, 3 insertions(+), 5 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 273a62da..13307d4b 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -2214,12 +2214,15 @@ def run_icing_predict(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', output_dir=h
['solar_zenith_angle', 'sensor_zenith_angle', 'cld_height_acha', 'cld_geo_thick'],
satellite=satellite, domain=domain)
+ satzen = ancil_data_dct['sensor_zenith_angle']
solzen = ancil_data_dct['solar_zenith_angle']
day_idxs = []
nght_idxs = []
for j in range(num_lines):
for i in range(num_elems):
k = i + j*num_elems
+ if not check_oblique(satzen[k]):
+ continue
if is_day(solzen[k]):
day_idxs.append(k)
else:
@@ -2279,11 +2282,6 @@ def run_icing_predict(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', output_dir=h
fd_preds[nght_idxs] = nght_preds[:]
fd_probs[nght_idxs] = nght_probs[:]
- # 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)
for flvl in flight_levels:
fd_preds = preds_2d_dct[flvl]
--
GitLab