diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index bde3c47a66baec7c67585985957e050e6e5077a2..5fada8e08283789660940fdc6013e6cfc0483fca 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -990,7 +990,21 @@ def tile_extract(trnfile='/home/rink/tiles_l1b_train.h5', tstfile='/home/rink/ti
     for ds_name in train_params:
         data_dct[ds_name] = np.concatenate([icing_data_dct[ds_name], no_icing_data_dct[ds_name]])
 
-    trn_idxs, tst_idxs = split_data(icing_intensity.shape[0], shuffle=False, perc=split)
+    # 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]
+
+    #trn_idxs, tst_idxs = split_data(icing_intensity.shape[0], shuffle=False, perc=split)
+    all_idxs = np.arange(icing_intensity.shape[0])
+    splt_idx = int(icing_intensity.shape[0] * 0.8)
+    trn_idxs = all_idxs[0:splt_idx]
+    tst_idxs = all_idxs[splt_idx:]
+    # ---------------------------------------------
 
     trn_data_dct = {}
     for ds_name in train_params: