From 8513343a6561fa31475634656a127841e4f31114 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 20 Jun 2022 12:18:26 -0500
Subject: [PATCH] enhance, enhance...
---
modules/aeolus/aeolus_amv.py | 43 +++++++++++++++++++++++++++---------
1 file changed, 33 insertions(+), 10 deletions(-)
diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 1b6a758b..07a0dcf4 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -2062,11 +2062,11 @@ def create_bestfit_file(filename, match_dct, raob_dct, gfs_at_raob_dct, bf_dct,
# the top level group for the output file
rootgrp = Dataset(filename, 'w', format='NETCDF4')
- dim_aparams = rootgrp.createDimension('num_raob_params', size=num_aparams)
- dim_alevs = rootgrp.createDimension('max_num_raob_levels', size=max_num_alevels)
- dim_amvs = rootgrp.createDimension('max_num_amvs', size=max_num_amvs)
- dim_num_aeolus_prof = rootgrp.createDimension('num_raob_profs', size=num_aprofs)
- dim_time = rootgrp.createDimension('time', size=1)
+ rootgrp.createDimension('num_raob_params', size=num_aparams)
+ rootgrp.createDimension('max_num_raob_levels', size=max_num_alevels)
+ rootgrp.createDimension('max_num_amvs', size=max_num_amvs)
+ rootgrp.createDimension('num_raob_profs', size=num_aprofs)
+ rootgrp.createDimension('time', size=1)
time = rootgrp.createVariable('time', 'f4', ['time'])
time.units = 'seconds since 1970-01-1 00:00:00'
@@ -2084,7 +2084,7 @@ def create_bestfit_file(filename, match_dct, raob_dct, gfs_at_raob_dct, bf_dct,
prf_spd = rootgrp.createVariable('raob_spd', 'f4', ['time', 'num_raob_profs', 'max_num_raob_levels'])
prf_spd.units = 'm s-1'
prf_lvl = rootgrp.createVariable('raob_levels', 'f4', ['time', 'num_raob_profs', 'max_num_raob_levels'])
- prf_lvl.units = 'millibar'
+ prf_lvl.units = 'hPa'
# ----- Product variables ----------------
nc4_vars = []
@@ -2096,8 +2096,13 @@ def create_bestfit_file(filename, match_dct, raob_dct, gfs_at_raob_dct, bf_dct,
amv_lat.units = 'degrees_north'
nc4_vars.append(amv_lat)
- nc4_vars.append(rootgrp.createVariable('amv_elem', 'i4', ['time', 'num_raob_profs', 'max_num_amvs']))
- nc4_vars.append(rootgrp.createVariable('amv_line', 'i4', ['time', 'num_raob_profs', 'max_num_amvs']))
+ amv_elem = rootgrp.createVariable('amv_elem', 'i4', ['time', 'num_raob_profs', 'max_num_amvs'])
+ amv_elem.long_name = 'FGF x coordinate Full Disk GOES16'
+ nc4_vars.append(amv_elem)
+
+ amv_line = rootgrp.createVariable('amv_line', 'i4', ['time', 'num_raob_profs', 'max_num_amvs'])
+ amv_line.long_name = 'FGF y coordinate Full Disk GOES16'
+ nc4_vars.append(amv_line)
var_s = rg_exmpl.variables
for pidx, param in enumerate(amv_params):
@@ -2119,19 +2124,37 @@ def create_bestfit_file(filename, match_dct, raob_dct, gfs_at_raob_dct, bf_dct,
# Bestfit variables ------------------
var_u = rootgrp.createVariable('bf_raob_u', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+ var_u.standard_name = 'eastward_wind'
+ var_u.long_name = 'wind at bestfit level with raob as background profile'
var_u.units = 'm s-1'
+
var_v = rootgrp.createVariable('bf_raob_v', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+ var_v.standard_name = 'northward_wind'
+ var_v.long_name = 'wind at bestfit level with raob as background profile'
var_v.units = 'm s-1'
+
var_p = rootgrp.createVariable('bf_raob_press', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
- var_p.units = 'millibar'
+ var_p.standard_name = 'air_pressure'
+ var_p.long_name = 'pressure level of bestfit with raob as background profile'
+ var_p.units = 'hPa'
+
var_flg = rootgrp.createVariable('bf_raob_flag', 'i4', ['time', 'num_raob_profs', 'max_num_amvs'])
var_u_g = rootgrp.createVariable('bf_gfs_u', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+ var_u_g.standard_name = 'eastward_wind'
+ var_u_g.long_name = 'wind at bestfit level with raob as background profile'
var_u_g.units = 'm s-1'
+
var_v_g = rootgrp.createVariable('bf_gfs_v', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+ var_v_g.standard_name = 'northward_wind'
+ var_v_g.long_name = 'wind at bestfit level with raob as background profile'
var_v_g.units = 'm s-1'
+
var_p_g = rootgrp.createVariable('bf_gfs_press', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
- var_p_g.units = 'millibar'
+ var_p_g.standard_name = 'air_pressure'
+ var_p_g.long_name = 'pressure level of bestfit with raob as background profile'
+ var_p_g.units = 'hPa'
+
var_flg_g = rootgrp.createVariable('bf_gfs_flag', 'i4', ['time', 'num_raob_profs', 'max_num_amvs'])
# Write data to file ---------------------
--
GitLab