Skip to content
Snippets Groups Projects
Commit 8f455070 authored by rink's avatar rink
Browse files

add profile lon/lat to product output

parent d5bc7a35
No related branches found
No related tags found
No related merge requests found
......@@ -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
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment