From 74b7acdf47340f9da91bbf193989889f7fea7026 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Wed, 10 Mar 2021 21:21:55 -0600
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 29 +++++++++++++++++++++++++++--
1 file changed, 27 insertions(+), 2 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index d5faeec0..38164b7e 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -28,6 +28,8 @@ ds_list = ['temp_8_5um_nom', 'temp_10_4um_nom', 'temp_11_0um_nom', 'temp_13_3um_
'cloud_mask']
ds_grd_dct = {name: [] for name in ds_list}
+a_clvr_file = '/Users/tomrink/data/clavrx/clavrx_OR_ABI-L1b-RadC-M3C01_G16_s20190020002186.level2.nc'
+
def setup():
ice_dict, no_ice_dict = pirep_icing(pirep_file)
@@ -89,6 +91,17 @@ 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):
+ #h5f_expl = h5py.File(a_clvr_file, 'r')
+ hf5 = h5py.File(filename, 'w')
+
+ for ds_name in ds_list:
+ data = data_dct[ds_name]
+ hf5.create_dataset(ds_name, data=data)
+
+ hf5.close()
+
+
def run(pirep_dct, outfile=None):
time_keys = list(pirep_dct.keys())
@@ -144,8 +157,13 @@ def run(pirep_dct, outfile=None):
cnt += 1
+ data_dct = {}
+ for ds_name in ds_list:
+ data_dct[ds_name] = np.array(ds_grd_dct[ds_name])
+
print('num images: ', cnt, len(time_keys), miss_a, miss_b)
+ create_file(outfile, data_dct, ds_list)
def analyze(ice_dct, no_ice_dct):
@@ -210,5 +228,12 @@ def analyze(ice_dct, no_ice_dct):
print(dt_str[2:])
-def create_file(filename, ds_list, ds_types):
- pass
\ No newline at end of file
+def create_file(filename, data_dct, ds_list):
+ #h5f_expl = h5py.File(a_clvr_file, 'r')
+ hf5 = h5py.File(filename, 'w')
+
+ for ds_name in ds_list:
+ data = data_dct[ds_name]
+ hf5.create_dataset(ds_name, data=data)
+
+ hf5.close()
\ No newline at end of file
--
GitLab