diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 9c32175eaf61313edccfc953298630ad6a21adae..0f364e246ab42fc7ff840f04047464f79442f322 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1966,7 +1966,7 @@ flt_level_ranges[4] = [8000.0, 15000.0]
 def run_make_images(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', ckpt_dir_s_path='/Users/tomrink/tf_model/', prob_thresh=0.5, satellite='GOES16', domain='CONUS',
                     extent=[-105, -70, 15, 50],
                     pirep_file='/Users/tomrink/data/pirep/pireps_202109200000_202109232359.csv',
-                    obs_lons=None, obs_lats=None, obs_times=None, flight_level=None,
+                    obs_lons=None, obs_lats=None, obs_times=None, obs_alt=None, flight_level=None,
                     day_night='DAY', l1b_andor_l2='l2'):
 
     if pirep_file is not None:
@@ -2000,6 +2000,7 @@ def run_make_images(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', ckpt_dir_s_pat
             _, keep_lons, keep_lats, _ = time_filter_3(ice_dict, ts_0, ts_1, alt_lo, alt_hi)
         elif obs_times is not None:
             keep = np.logical_and(obs_times >= ts_0, obs_times < ts_1)
+            keep = np.where(keep, np.logical_and(obs_alt >= alt_lo, obs_alt < alt_hi), False)
             keep_lons = obs_lons[keep]
             keep_lats = obs_lats[keep]
         else: