diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 886e7d08af287f8e23f865e636cda22cb98f8ce0..9a61ffb041ca6cf63c1ed43ddf9673a2dcec7d37 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1750,7 +1750,7 @@ def time_filter_2(times, dt_str_0=None, dt_str_1=None, format_code='%Y-%m-%d_%H:
     return keep_times, keep_idxs
 
 
-def time_filter_3(icing_dct, ts_0, ts_1):
+def time_filter_3(icing_dct, ts_0, ts_1, alt_lo=None, alt_hi=None):
     keep_reports = []
     keep_times = []
     keep_lons = []
@@ -1760,10 +1760,12 @@ def time_filter_3(icing_dct, ts_0, ts_1):
         if ts_0 <= ts < ts_1:
             rpts = icing_dct[ts]
             for idx, tup in enumerate(rpts):
-                keep_reports.append(tup)
-                keep_times.append(ts)
-                keep_lats.append(tup[0])
-                keep_lons.append(tup[1])
+                falt = tup[2]
+                if alt_lo is not None and (alt_lo < falt <= alt_hi):
+                    keep_reports.append(tup)
+                    keep_times.append(ts)
+                    keep_lats.append(tup[0])
+                    keep_lons.append(tup[1])
 
     return keep_times, keep_lons, keep_lats, keep_reports
 
@@ -1850,6 +1852,14 @@ def tiles_info(filename):
     print('Icing 6:  ', np.sum(iint == 6))
 
 
+flt_level_ranges = {k:None for k in range(5)}
+flt_level_ranges[0] = [0.0, 2000.0]
+flt_level_ranges[1] = [2000.0, 4000.0]
+flt_level_ranges[2] = [4000.0, 6000.0]
+flt_level_ranges[3] = [6000.0, 8000.0]
+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',
@@ -1863,6 +1873,8 @@ def run_make_images(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', ckpt_dir_s_pat
         clvrx_ds = CLAVRx(clvrx_dir)
     clvrx_files = clvrx_ds.flist
 
+    alt_lo, alt_hi = flt_level_ranges[flight_level]
+
     for fname in clvrx_files:
         h5f = h5py.File(fname, 'r')
         dto = clvrx_ds.get_datetime(fname)
@@ -1875,7 +1887,7 @@ def run_make_images(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', ckpt_dir_s_pat
         ts_0 = dto_0.timestamp()
         ts_1 = dto_1.timestamp()
         if pirep_file is not None:
-            _, keep_lons, keep_lats, _ = time_filter_3(ice_dict, ts_0, ts_1)
+            _, 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_lons = obs_lons[keep]