diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 273a62dae2ad36b9c5e10147ddce99e4f3ecc0f9..13307d4ba567e0bc00f57314515a54efa549f07e 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]