diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index b9536fa55b6df596060ee8d8ed125a8287d1bcc6..07562157b4f712d6ddaf996bf61efa89a957d6e7 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -950,7 +950,7 @@ def tile_extract(trnfile='/home/rink/tiles_l1b_train.h5', tstfile='/home/rink/ti
     trn_data_dct = {}
     for ds_name in train_params:
         trn_data_dct[ds_name] = data_dct[ds_name][trn_idxs,]
-    trn_icing_intesity = icing_intensity[trn_idxs,]
+    trn_icing_intensity = icing_intensity[trn_idxs,]
     trn_icing_times = icing_times[trn_idxs,]
     trn_icing_lons = icing_lons[trn_idxs,]
     trn_icing_lats = icing_lats[trn_idxs,]
@@ -958,13 +958,13 @@ def tile_extract(trnfile='/home/rink/tiles_l1b_train.h5', tstfile='/home/rink/ti
     #  Data augmentation -------------------------------------------------------------
     if augment:
         trn_data_dct_aug = {ds_name: [] for ds_name in train_params}
-        trn_icing_intesity_aug = []
+        trn_icing_intensity_aug = []
         trn_icing_times_aug = []
         trn_icing_lons_aug = []
         trn_icing_lats_aug = []
 
-        for k in range(trn_icing_intesity.shape[0]):
-            iceint = trn_icing_intesity[k]
+        for k in range(trn_icing_intensity.shape[0]):
+            iceint = trn_icing_intensity[k]
             icetime = trn_icing_times[k]
             icelon = trn_icing_lons[k]
             icelat = trn_icing_lats[k]
@@ -975,9 +975,9 @@ def tile_extract(trnfile='/home/rink/tiles_l1b_train.h5', tstfile='/home/rink/ti
                     trn_data_dct_aug[ds_name].append(np.flipud(dat[k,]))
                     trn_data_dct_aug[ds_name].append(np.rot90(dat[k,]))
 
-                trn_icing_intesity_aug.append(iceint)
-                trn_icing_intesity_aug.append(iceint)
-                trn_icing_intesity_aug.append(iceint)
+                trn_icing_intensity_aug.append(iceint)
+                trn_icing_intensity_aug.append(iceint)
+                trn_icing_intensity_aug.append(iceint)
 
                 trn_icing_times_aug.append(icetime)
                 trn_icing_times_aug.append(icetime)
@@ -993,14 +993,14 @@ def tile_extract(trnfile='/home/rink/tiles_l1b_train.h5', tstfile='/home/rink/ti
 
         for ds_name in train_params:
             trn_data_dct_aug[ds_name] = np.stack(trn_data_dct_aug[ds_name])
-        trn_icing_intesity_aug = np.stack(trn_icing_intesity_aug)
-        trn_icing_times_aug = np.stack(trn_icing_intesity_aug)
+        trn_icing_intensity_aug = np.stack(trn_icing_intensity_aug)
+        trn_icing_times_aug = np.stack(trn_icing_intensity_aug)
         trn_icing_lons_aug = np.stack(trn_icing_lons_aug)
         trn_icing_lats_aug = np.stack(trn_icing_lats_aug)
 
         for ds_name in train_params:
             trn_data_dct[ds_name] = np.concatenate([trn_data_dct[ds_name], trn_data_dct_aug[ds_name]])
-        trn_icing_intensity = np.concatenate([trn_icing_intesity, trn_icing_intesity_aug])
+        trn_icing_intensity = np.concatenate([trn_icing_intensity, trn_icing_intensity_aug])
         trn_icing_times = np.concatenate([trn_icing_times, trn_icing_times_aug])
         trn_icing_lons = np.concatenate([trn_icing_lons, trn_icing_lons_aug])
         trn_icing_lats = np.concatenate([trn_icing_lats, trn_icing_lats_aug])