From ebdac2565ea451ba60eb6b52959c82cb8968c274 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sat, 20 Mar 2021 10:40:45 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 9 ++++++---
1 file changed, 6 insertions(+), 3 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index a48f6ca9..19baeaf6 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -101,7 +101,7 @@ def get_grid_values(h5f, grid_name, j_c, i_c, half_width, scale_factor_name='sca
return grd_vals
-def create_file(filename, data_dct, ds_list, lon_c, lat_c, time_s, fl_alt_s, icing_intensity=None):
+def create_file(filename, data_dct, ds_list, lon_c, lat_c, time_s, fl_alt_s, icing_intensity):
h5f_expl = h5py.File(a_clvr_file, 'r')
h5f = h5py.File(filename, 'w')
@@ -168,6 +168,7 @@ def run(pirep_dct, outfile=None, outfile_l1b=None, dt_str_start=None, dt_str_end
lat_c = []
time_s = []
fl_alt_s = []
+ ice_int_s = []
for idx, time in enumerate(time_keys):
if t_start is not None:
if time < t_start:
@@ -233,6 +234,7 @@ def run(pirep_dct, outfile=None, outfile_l1b=None, dt_str_start=None, dt_str_end
lat_c.append(lat_s[0])
time_s.append(time)
fl_alt_s.append(fl)
+ ice_int_s.append(I)
if len(time_s) == 0:
return
@@ -247,10 +249,11 @@ def run(pirep_dct, outfile=None, outfile_l1b=None, dt_str_start=None, dt_str_end
lat_c = np.array(lat_c)
time_s = np.array(time_s)
fl_alt_s = np.array(fl_alt_s)
+ ice_int_s = np.array(ice_int_s)
if outfile is not None:
outfile = add_time_range_to_filename(outfile, t_start, t_end)
- create_file(outfile, data_dct, ds_list, lon_c, lat_c, time_s, fl_alt_s)
+ create_file(outfile, data_dct, ds_list, lon_c, lat_c, time_s, fl_alt_s, ice_int_s)
data_dct = {}
for ds_name in l1b_ds_list:
@@ -258,7 +261,7 @@ def run(pirep_dct, outfile=None, outfile_l1b=None, dt_str_start=None, dt_str_end
if outfile_l1b is not None:
outfile_l1b = add_time_range_to_filename(outfile_l1b, t_start, t_end)
- create_file(outfile_l1b, data_dct, l1b_ds_list, lon_c, lat_c, time_s, fl_alt_s)
+ create_file(outfile_l1b, data_dct, l1b_ds_list, lon_c, lat_c, time_s, fl_alt_s, ice_int_s)
def analyze2(ice_dct):
--
GitLab