From 8f455070549c6d416111c90be7c7695abbe6e294 Mon Sep 17 00:00:00 2001
From: rink <rink@ssec.wisc.edu>
Date: Thu, 1 Oct 2020 13:33:04 -0500
Subject: [PATCH] add profile lon/lat to product output
---
modules/aeolus/aeolus_amv.py | 5 +++++
1 file changed, 5 insertions(+)
diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 97a4dd48..335c2efa 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -378,6 +378,8 @@ def create_file(filename, aeolus_to_amv_dct, aeolus_dct):
num_amvs_per_prof = rootgrp.createVariable('num_amvs_per_prof', 'i4', ['num_aeolus_profs'])
num_levs_per_prof = rootgrp.createVariable('num_levs_per_prof', 'i4', ['num_aeolus_profs'])
prof_time = rootgrp.createVariable('time', 'f4', ['num_aeolus_profs'])
+ prf_lon = rootgrp.createVariable('prof_longitude', 'f4', ['num_aeolus_profs'])
+ prf_lat = rootgrp.createVariable('prof_latitude', 'f4', ['num_aeolus_profs'])
prof_time.units = 'seconds since 1970-01-1 00:00:00'
prf_azm = rootgrp.createVariable('prof_azm', 'f4', ['profs'])
@@ -414,6 +416,9 @@ def create_file(filename, aeolus_to_amv_dct, aeolus_dct):
prf_spd[i_c:i_d] = prof[:, 6]
i_c += nlevs
+ prf_lat[idx] = prof[0, 0]
+ prf_lon[idx] = prof[0, 1]
+
num_amvs_per_prof[:] = num_amvs
num_levs_per_prof[:] = num_levs
prof_time[:] = times
--
GitLab