From a2f160ceda37b2ebe11cc998ae83ce4e4ca948d5 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Wed, 14 Apr 2021 16:35:08 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index bc94aea8..1239d075 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, ')')
--
GitLab