Skip to content
Snippets Groups Projects
Commit a5470e32 authored by tomrink's avatar tomrink
Browse files

add new match file creator

parent 9891d473
No related branches found
No related tags found
No related merge requests found
...@@ -1455,7 +1455,8 @@ def match_aeolus_to_clavrx(aeolus_dict, clvrx_files): ...@@ -1455,7 +1455,8 @@ def match_aeolus_to_clavrx(aeolus_dict, clvrx_files):
cc, ll = nav.earth_to_lc(lon, lat) cc, ll = nav.earth_to_lc(lon, lat)
if cc is None or ll is None: if cc is None or ll is None:
continue continue
if cc - 4 < 0 or ll - 4 < 0 or cc + 5 >= num_elems or ll + 5 >= num_lines:
continue
# c_rng, l_rng = get_search_box(nav, lon, lat) # c_rng, l_rng = get_search_box(nav, lon, lat)
# if c_rng is None: # if c_rng is None:
# continue # continue
...@@ -1545,15 +1546,18 @@ def create_file_2(match_dct, filename, clvrx_params=['cld_height_acha', 'cld_pre ...@@ -1545,15 +1546,18 @@ def create_file_2(match_dct, filename, clvrx_params=['cld_height_acha', 'cld_pre
var.units = u var.units = u
nc4_vars.append(var) nc4_vars.append(var)
# Write data to file... # Write data to file ---------------------
prf_lon[:] = alons
prf_lat[:] = alats
prf_elem[:] = elems
prf_line[:] = lines
prf_time[:] = atimes
for idx, key in enumerate(keys): for idx, key in enumerate(keys):
tup = match_dct.get(key) tup = match_dct.get(key)
prof = tup[4] prof = tup[4]
param_nd = tup[5] param_nd = tup[5]
for pidx, param in enumerate(clvrx_params):
nc4_vars[pidx][idx, :, :] = param_nd[pidx, :, :]
nlevs = prof.shape[0] nlevs = prof.shape[0]
for k in range(nlevs): for k in range(nlevs):
prf_spd[idx,k] = prof[k,6] prf_spd[idx,k] = prof[k,6]
...@@ -1561,11 +1565,9 @@ def create_file_2(match_dct, filename, clvrx_params=['cld_height_acha', 'cld_pre ...@@ -1561,11 +1565,9 @@ def create_file_2(match_dct, filename, clvrx_params=['cld_height_acha', 'cld_pre
prf_hht[idx,k] = prof[k,3] prf_hht[idx,k] = prof[k,3]
prf_hhb[idx,k] = prof[k,4] prf_hhb[idx,k] = prof[k,4]
prf_lon[:] = alons for pidx, param in enumerate(clvrx_params):
prf_lat[:] = alats nc4_vars[pidx][idx, :, :] = param_nd[pidx, :, :]
prf_elem[:] = elems
prf_line[:] = lines
prf_time[:] = atimes
rootgrp.close() rootgrp.close()
... ...
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment