From 83f2629c6f268d04dddaf203e9d8e537c211aad6 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Tue, 13 Apr 2021 11:36:18 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 17 ++++++++++++++++-
1 file changed, 16 insertions(+), 1 deletion(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index f39a8686..e507ad1c 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
--
GitLab