diff --git a/modules/icing/util.py b/modules/icing/util.py
index 31d8014227115ffb394a94d654ac7335989ef563..20ce4232549c3095d5848df12315e781fe60cbc3 100644
--- a/modules/icing/util.py
+++ b/modules/icing/util.py
@@ -1266,7 +1266,7 @@ def run_icing_predict_image(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', output
             keep = np.where(keep, obs_intensity >= icing_intensity_plot_min, False)
             keep_lons = obs_lons[keep]
             keep_lats = obs_lats[keep]
-            obs_intensity = obs_intensity[keep]
+            keep_obs_intensity = obs_intensity[keep]
         else:
             keep_lons = None
             keep_lats = None
@@ -1281,7 +1281,7 @@ def run_icing_predict_image(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', output
         max_prob = np.where(max_prob < prob_thresh, np.nan, max_prob)
 
         make_icing_image(h5f, max_prob, None, None, clvrx_str_time, satellite, domain,
-                         ice_lons_vld=keep_lons, ice_lats_vld=keep_lats, ice_intensity=obs_intensity, extent=extent)
+                         ice_lons_vld=keep_lons, ice_lats_vld=keep_lats, ice_intensity=keep_obs_intensity, extent=extent)
 
         print('Done: ', clvrx_str_time)
         h5f.close()