diff --git a/modules/icing/util.py b/modules/icing/util.py
index 618a1ca9955a8d2d2768ec79ad112f69c28a97fd..e9523e8be863d5bedd9c73a01180c8614c4d1241 100644
--- a/modules/icing/util.py
+++ b/modules/icing/util.py
@@ -98,7 +98,7 @@ def run_icing_predict(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', output_dir=h
         bt_10_4 = ancil_data_dct['temp_10_4um_nom']
         day_idxs = []
         nght_idxs = []
-        avg_bts = []
+        avg_bt = []
         for j in range(num_lines):
             for i in range(num_elems):
                 k = i + j*num_elems
@@ -108,7 +108,7 @@ def run_icing_predict(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', output_dir=h
                     day_idxs.append(k)
                 else:
                     nght_idxs.append(k)
-                avg_bts.append(get_average(bt_10_4[k]))
+                avg_bt.append(get_average(bt_10_4[k]))
 
         num_tiles = num_lines * num_elems
         num_day_tiles = len(day_idxs)
@@ -179,7 +179,8 @@ def run_icing_predict(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', output_dir=h
             preds_2d_dct[flvl] = fd_preds.reshape((num_lines, num_elems))
             probs_2d_dct[flvl] = fd_probs.reshape((num_lines, num_elems))
 
-        bt_10_4_2d = bt_10_4.reshape((num_lines, num_elems))
+        avg_bt = np.array(avg_bt)
+        bt_10_4_2d = avg_bt.reshape((num_lines, num_elems))
 
         write_icing_file_nc4(clvrx_str_time, output_dir, preds_2d_dct, probs_2d_dct,
                              x_rad, y_rad, lons_2d, lats_2d, cc, ll,