diff --git a/modules/util/util.py b/modules/util/util.py
index 47df416bd4a94efd24f02418597aa5a1e4c94d4a..be0a25f6c37b5091856b4a1c832036cf40d4352d 100644
--- a/modules/util/util.py
+++ b/modules/util/util.py
@@ -812,7 +812,7 @@ def write_icing_file_nc4(clvrx_str_time, output_dir, preds_dct, probs_dct,
         icing_pred_ds.setncattr('grid_mapping', 'Projection')
         icing_pred_ds.setncattr('missing', -1)
         if has_time:
-            preds = preds.reshap((1, y.shape[0], x.shape[0]))
+            preds = preds.reshape((1, y.shape[0], x.shape[0]))
         icing_pred_ds[:,] = preds
 
     for flvl in flt_lvls:
@@ -824,7 +824,7 @@ def write_icing_file_nc4(clvrx_str_time, output_dir, preds_dct, probs_dct,
         icing_prob_ds.setncattr('grid_mapping', 'Projection')
         icing_prob_ds.setncattr('missing', -1.0)
         if has_time:
-            probs = probs.reshap((1, y.shape[0], x.shape[0]))
+            probs = probs.reshape((1, y.shape[0], x.shape[0]))
         icing_prob_ds[:,] = probs
 
     prob_s = np.stack(prob_s, axis=-1)
@@ -840,7 +840,7 @@ def write_icing_file_nc4(clvrx_str_time, output_dir, preds_dct, probs_dct,
 
     max_lvl = np.argmax(prob_s, axis=2)
     if has_time:
-        max_lvl = max_lvl.reshap((1, y.shape[0], x.shape[0]))
+        max_lvl = max_lvl.reshape((1, y.shape[0], x.shape[0]))
 
     icing_pred_ds = rootgrp.createVariable('max_icing_probability_level', 'i2', var_dim_list)
     icing_pred_ds.setncattr('coordinates', geo_coords)