diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index d5dd60b80045e1e07f1e297196c8857efa0a2307..6886ebb1275dbc763b2eaf3f0450ff557182d4ba 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1415,10 +1415,10 @@ def run_mean_std(check_cloudy=False, no_icing_to_icing_ratio=5):
             keep = np.logical_or(cld_msk_i == 2, cld_msk_i == 3)
             data_i = data_i[keep]
         mean_i = np.nanmean(data_i)
-        data_i -= mean_i
-        std_i = np.nanstd(data_i)
         lo_i = np.nanmin(data_i)
         hi_i = np.nanmax(data_i)
+        data_i -= mean_i
+        std_i = np.nanstd(data_i)
         cnt_i = np.sum(np.invert(np.isnan(data_i)))
 
         data_ni = np.concatenate(data_ni)
@@ -1426,10 +1426,10 @@ def run_mean_std(check_cloudy=False, no_icing_to_icing_ratio=5):
             keep = np.logical_or(cld_msk_ni == 2, cld_msk_ni == 3)
             data_ni = data_ni[keep]
         mean_ni = np.nanmean(data_ni)
-        data_ni -= mean_ni
-        std_ni = np.nanstd(data_ni)
         lo_ni = np.nanmin(data_ni)
         hi_ni = np.nanmax(data_ni)
+        data_ni -= mean_ni
+        std_ni = np.nanstd(data_ni)
         cnt_ni = np.sum(np.invert(np.isnan(data_ni)))
 
         no_icing_to_icing_ratio = cnt_ni/cnt_i