diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 4baa96e40ee4d1ca3fbf4a3ca7aff908655242b0..0c157f205c8e1fcaacd0dbe48b18a48ff2065935 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)