diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 058e5110ff4f1b55799ef089a7cda1ee322fd4e0..371ffc307dc4148e22cbbfc4956e6a2c87f7e61e 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -444,6 +444,10 @@ def process(ice_dct, no_ice_dct, neg_ice_dct):
     ice_keys = np.concatenate([ice_keys_5_6, ice_keys_1, ice_keys_2, ice_keys_3, ice_keys_4])
     ice_tidx = np.concatenate([ice_tidx_5_6, ice_tidx_1, ice_tidx_2, ice_tidx_3, ice_tidx_4])
 
+    sidxs = np.argsort(ice_keys)
+    ice_keys = ice_keys[sidxs]
+    ice_tidx = ice_tidx[sidxs]
+
     for idx, key in enumerate(ice_keys):
         rpts = ice_dct[key]
         tup = rpts[ice_tidx[idx]]
@@ -466,6 +470,10 @@ def process(ice_dct, no_ice_dct, neg_ice_dct):
     no_ice_tidx = no_ice_tidx[::20]
     print('no ice   : ', no_ice_keys.shape[0])
 
+    sidxs = np.argsort(no_ice_keys)
+    no_ice_keys = no_ice_keys[sidxs]
+    no_ice_tidx = no_ice_tidx[sidxs]
+
     for idx, key in enumerate(no_ice_keys):
         rpts = no_ice_dct[key]
         tup = rpts[no_ice_tidx[idx]]