From 877999e61c1999260d8a12b874169234966f36a0 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Tue, 9 Mar 2021 21:30:30 -0600
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 22 ++--------------------
1 file changed, 2 insertions(+), 20 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 5bf904ce..d8bf21b4 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -19,7 +19,8 @@ goes_ds_dct = {}
ice_dict = None
no_ice_dict = None
time_keys = None
-pirep_file = '/home/rink/data/pireps/pireps_2019010000_2019063023.csv'
+#pirep_file = '/home/rink/data/pireps/pireps_2019010000_2019063023.csv'
+pirep_file = '/home/rink/data/pireps/pireps_20160101_20190630.csv'
ds_list = ['temp_8_5um_nom', 'temp_10_4um_nom', 'temp_11_0um_nom', 'temp_13_3um_nom', 'cld_height_acha',
'supercooled_cloud_fraction', 'cld_temp_acha', 'cld_opd_acha', 'cloud_phase', 'solar_zenith_angle',
@@ -147,25 +148,6 @@ def run(pirep_dct, outfile=None):
def analyze(ice_dct, no_ice_dct):
- ice_times = np.array(list(ice_dct.keys()))
- no_ice_times = np.array(list(no_ice_dct.keys()))
-
- dt_str = '201601010000'
- dto = datetime.datetime.strptime(dt_str, '%Y%m%d%H%M').replace(tzinfo=timezone.utc)
- t0 = dto.timestamp()
-
- dt_str = '201906300000'
- dto = datetime.datetime.strptime(dt_str, '%Y%m%d%H%M').replace(tzinfo=timezone.utc)
- t1 = dto.timestamp()
-
- tline = np.arange(t0, t1, 600)
- hist_a, edges = np.histogram(no_ice_times, bins=tline)
- print(len(np.nonzero(hist_a)[0]))
-
- hist_b, edges = np.histogram(ice_times, bins=tline)
- print(len(np.nonzero(hist_b)[0]))
-
- print(np.sum(np.logical_and(hist_a > 0, hist_b > 0)))
last_file = None
ice_files = []
--
GitLab