diff --git a/modules/amv/caliop_clavrx_amv.py b/modules/amv/caliop_clavrx_amv.py index 57c925731648ae2588dd6072be4f0da665b026b1..657cedddb2a0b6525cff9cd1cba5dd2f8d528b94 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=',')