diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index f39a86866eebae92b7d0444ec36ea6d7d9bde18d..e507ad1c24230cdc5554d276c86e4da5b01502c7 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1144,17 +1144,32 @@ def run_mean_std():
 
     for dname in ds_list:
         data = []
+        data_i = []
+        data_ni = []
         for idx, ice_h5f in enumerate(ice_h5f_lst):
             no_ice_h5f = no_ice_h5f_lst[idx]
             data.append(ice_h5f[dname][:,].flatten())
             data.append(no_ice_h5f[dname][:,].flatten())
 
+            data_i.append(ice_h5f[dname][:,].flatten())
+            data_ni.append(no_ice_h5f[dname][:,].flatten())
+
         data = np.concatenate(data)
         mean = np.nanmean(data)
         data -= mean
         std = np.nanstd(data)
 
-        print(dname,': ', mean, std)
+        data_i = np.concatenate(data_i)
+        mean_i = np.nanmean(data_i)
+        data_i -= mean_i
+        std_i = np.nanstd(data_i)
+
+        data_ni = np.concatenate(data_ni)
+        mean_ni = np.nanmean(data_ni)
+        data_ni -= mean_ni
+        std_ni = np.nanstd(data_ni)
+
+        print(dname,': ', mean, mean_i, mean_ni, std, std_i, std_ni)
 
         mean_dct[dname] = mean
         std_dct[dname] = std