From 1b6cd85d531633301e7042e52be2fc86ed4675b4 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Thu, 14 Jul 2022 12:51:22 -0500
Subject: [PATCH] fix some cf issues
---
modules/aeolus/aeolus_amv.py | 20 ++++++++++++++++++++
1 file changed, 20 insertions(+)
diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index df838a49..6a28b4dc 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -2169,6 +2169,21 @@ def create_bestfit_file(filename, match_dct, raob_dct, raob_ts, gfs_at_raob_dct,
var_flg_g = rootgrp.createVariable('bf_gfs_flag', 'i4', ['time', 'num_raob_profs', 'max_num_amvs'])
+ var_u_cp = rootgrp.createVariable('closest_press_u', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+ var_u_cp.standard_name = 'eastward_wind'
+ var_u_cp.long_name = 'wind at closest raob pressure level to amv within +/- 60mb'
+ var_u_cp.units = 'm s-1'
+
+ var_v_cp = rootgrp.createVariable('closest_press_v', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+ var_v_cp.standard_name = 'northward_wind'
+ var_v_cp.long_name = 'wind at closest raob pressure level to amv within +/- 60mb'
+ var_v_cp.units = 'm s-1'
+
+ var_p_cp = rootgrp.createVariable('closest_press', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+ var_p_cp.standard_name = 'air_pressure'
+ var_p_cp.long_name = 'closest raob pressure level to amv within +/- 60mb'
+ var_p_cp.units = 'hPa'
+
# Write data to file ---------------------
time[0] = raob_ts
prf_lon[0, :] = alons
@@ -2196,6 +2211,11 @@ def create_bestfit_file(filename, match_dct, raob_dct, raob_ts, gfs_at_raob_dct,
var_v[0, idx, 0:num_amvs] = bf_raob[:, 1]
var_p[0, idx, 0:num_amvs] = bf_raob[:, 2]
+ cp_raob = bf_tup[2]
+ var_u_cp[0, idx, 0:num_amvs] = cp_raob[:, 0]
+ var_v_cp[0, idx, 0:num_amvs] = cp_raob[:, 1]
+ var_p_cp[0, idx, 0:num_amvs] = cp_raob[:, 2]
+
bf_gfs = bf_tup[3]
var_u_g[0, idx, 0:num_amvs] = bf_gfs[:, 0]
var_v_g[0, idx, 0:num_amvs] = bf_gfs[:, 1]
--
GitLab