diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 1239d07510b6135cb3af0e96914d10b5222ca18f..4fa25560b5389e469e570923cccebc06b5a590cc 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -940,6 +940,15 @@ def fov_extract(outfile='/home/rink/fovs_out.h5', train_params=train_params_day)
     # icing_lons = icing_lons[ds_indexes]
     # icing_lats = icing_lats[ds_indexes]
 
+    # do sort
+    ds_indexes = np.argsort(icing_times)
+    for ds_name in train_params:
+        data_dct[ds_name] = data_dct[ds_name][ds_indexes]
+    icing_intensity = icing_intensity[ds_indexes]
+    icing_times = icing_times[ds_indexes]
+    icing_lons = icing_lons[ds_indexes]
+    icing_lats = icing_lats[ds_indexes]
+
     h5f_expl = h5py.File(a_clvr_file, 'r')
     h5f_out = h5py.File(outfile, 'w')
 
@@ -1084,6 +1093,15 @@ def tile_extract(outfile='/home/rink/tiles_out.h5', train_params=train_params_da
     # icing_lons = icing_lons[ds_indexes]
     # icing_lats = icing_lats[ds_indexes]
 
+    # do sort
+    ds_indexes = np.argsort(icing_times)
+    for ds_name in train_params:
+        data_dct[ds_name] = data_dct[ds_name][ds_indexes]
+    icing_intensity = icing_intensity[ds_indexes]
+    icing_times = icing_times[ds_indexes]
+    icing_lons = icing_lons[ds_indexes]
+    icing_lats = icing_lats[ds_indexes]
+
     h5f_expl = h5py.File(a_clvr_file, 'r')
     h5f_out = h5py.File(outfile, 'w')