diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index bc94aea806698e2fc918afb93b6d477f62b6cbdd..1239d07510b6135cb3af0e96914d10b5222ca18f 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1131,7 +1131,7 @@ def tile_extract(outfile='/home/rink/tiles_out.h5', train_params=train_params_da
     h5f_expl.close()
 
 
-def run_mean_std(check_cloudy=False):
+def run_mean_std(check_cloudy=False, no_icing_to_icing_ratio=5):
     ds_list = ['cld_height_acha', 'cld_geo_thick', 'cld_press_acha',
                'supercooled_cloud_fraction', 'cld_temp_acha', 'cld_opd_acha',
                'cld_reff_acha', 'cld_reff_dcomp', 'cld_reff_dcomp_1', 'cld_reff_dcomp_2', 'cld_reff_dcomp_3',
@@ -1190,8 +1190,8 @@ def run_mean_std(check_cloudy=False):
         data_ni -= mean_ni
         std_ni = np.nanstd(data_ni)
 
-        mean = (mean_i + 5*mean_ni)/6
-        std = (std_i + 5*std_ni)/6
+        mean = (mean_i + no_icing_to_icing_ratio*mean_ni)/(no_icing_to_icing_ratio + 1)
+        std = (std_i + no_icing_to_icing_ratio*std_ni)/(no_icing_to_icing_ratio + 1)
 
         print(dname,': (', mean, mean_i, mean_ni, ') (', std, std_i, std_ni, ')')