diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 25a0873489208dc58ad0d8302aa965d56f5f527b..0c82802129ba2fd3e38b07ba37c9ee3c05f4caa5 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -318,6 +318,7 @@ def analyze(ice_dct, no_ice_dct):
     ice_indexes = np.arange(len(ice_times))
 
     ucomm2 = np.setxor1d(comm2, ice_indexes)
+    np.random.seed(42)
     np.random.shuffle(ucomm2)
     ucomm2 = ucomm2[0:8000]
 
@@ -446,7 +447,9 @@ def process_1(ice_dct, no_ice_dct, neg_ice_dct):
     ice_times_2 = np.array(ice_times_2)
     ice_keys_2 = np.array(ice_keys_2)
     print('2: ', ice_times_2.shape)
+    np.random.seed(42)
     np.random.shuffle(ice_times_2)
+    np.random.seed(42)
     np.random.shuffle(ice_keys_2)
     ice_keys_2 = ice_keys_2[0:30000]
 
@@ -470,6 +473,7 @@ def process_1(ice_dct, no_ice_dct, neg_ice_dct):
     print('no ice: ', no_ice_times.shape)
     print('neg ice: ', neg_ice_times.shape)
     no_ice_keys = np.array(no_ice_keys)
+    np.random.seed(42)
     np.random.shuffle(no_ice_keys)
     no_ice_keys = no_ice_keys[0:50000]
     uniq_sorted_no_ice = np.unique(no_ice_keys)
@@ -480,6 +484,7 @@ def process_1(ice_dct, no_ice_dct, neg_ice_dct):
         new_no_ice_dct[key] = no_ice_dct[key]
 
     neg_ice_keys = np.array(neg_ice_keys)
+    np.random.seed(42)
     np.random.shuffle(neg_ice_keys)
     neg_ice_keys = neg_ice_keys[0:5000]
     uniq_sorted_neg_ice = np.unique(neg_ice_keys)