snapshot...

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, ')')