From e98c5a7ed35905a2cec70b29a36577c680cb28a4 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Wed, 21 Sep 2022 12:14:10 -0500
Subject: [PATCH] snapshot...
---
modules/amv/caliop_clavrx_amv.py | 9 +++++++++
1 file changed, 9 insertions(+)
diff --git a/modules/amv/caliop_clavrx_amv.py b/modules/amv/caliop_clavrx_amv.py
index 57c92573..657ceddd 100644
--- a/modules/amv/caliop_clavrx_amv.py
+++ b/modules/amv/caliop_clavrx_amv.py
@@ -417,6 +417,15 @@ def compute_and_add_geo_hgt(path_to_amvs, path_to_gfs, amv_source, band='14', ou
alt_f = np.reshape(alt_f, [num_amvs, 1])
amvs_nd = np.hstack([amvs_nd, alt_f])
+ amv_params = amv_ds.get_parameters()
+ all_params = [amv_ds.lon_name, amv_ds.lat_name, amv_ds.elem_name, amv_ds.line_name] + amv_params
+ coords = {'num_params': all_params}
+ dims = ['num_amvs', 'num_params']
+ amvs_da = xr.DataArray(amvs_nd, coords=coords, dims=dims)
+ new_amv_fname = os.path.split(fname)[1]
+ new_amv_fname = os.path.splitext(new_amv_fname)[0] + 'nc4'
+ amvs_da.to_netcdf(new_amv_fname, format='NETCDF4')
+
new_amv_fname = os.path.split(fname)[1]
np.savetxt(new_amv_fname, amvs_nd, fmt=num_fmts, header=new_header, delimiter=',')
--
GitLab