From 9d4f3c9eab7bc81047965cd4527ec3ce24b3b925 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 25 Apr 2022 14:31:33 -0500
Subject: [PATCH] snapshot...
---
modules/util/util.py | 12 +++++++++++-
1 file changed, 11 insertions(+), 1 deletion(-)
diff --git a/modules/util/util.py b/modules/util/util.py
index d14a0534..4384499e 100644
--- a/modules/util/util.py
+++ b/modules/util/util.py
@@ -274,6 +274,11 @@ def check_oblique(satzen, test_angle=70.0):
return True
+def get_average(tile_2d):
+ tile = tile_2d.flatten()
+ return np.nanmean(tile)
+
+
def get_grid_values_all(h5f, grid_name, scale_factor_name='scale_factor', add_offset_name='add_offset',
fill_value_name='_FillValue', range_name='actual_range', fill_value=None):
hfds = h5f[grid_name]
@@ -842,7 +847,7 @@ def write_icing_file(clvrx_str_time, output_dir, preds_dct, probs_dct, x, y, lon
h5f_out.close()
-def write_icing_file_nc4(clvrx_str_time, output_dir, preds_dct, probs_dct,
+def write_icing_file_nc4(clvrx_str_time, output_dir, preds_dct, probs_dct, bt_10_4,
x, y, lons, lats, elems, lines, satellite='GOES16', domain='CONUS',
has_time=False, use_nan=False, prob_thresh=0.5):
outfile_name = output_dir + 'icing_prediction_'+clvrx_str_time+'.nc'
@@ -928,6 +933,11 @@ def write_icing_file_nc4(clvrx_str_time, output_dir, preds_dct, probs_dct,
icing_pred_ds.setncattr('missing', -1)
icing_pred_ds[:,] = max_lvl
+ bt_ds = rootgrp.createVariable('bt_10_4', 'f4', [dim_1_name, dim_0_name])
+ bt_ds.setncattr('coordinates', geo_coords)
+ bt_ds.setncattr('grid_mapping', 'Projection')
+ bt_ds[:,] = bt_10_4
+
lon_ds = rootgrp.createVariable('longitude', 'f4', [dim_1_name, dim_0_name])
lon_ds.units = 'degrees_east'
lon_ds[:,] = lons
--
GitLab