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

snapshot...

parent f65306f8
No related branches found
No related tags found
No related merge requests found
...@@ -1557,6 +1557,7 @@ def create_aeolus_clavrx_match_file(match_dct, filename, clvrx_params): ...@@ -1557,6 +1557,7 @@ def create_aeolus_clavrx_match_file(match_dct, filename, clvrx_params):
atimes = [] atimes = []
elems = [] elems = []
lines = [] lines = []
fidxs = []
# scan to get max num levels # scan to get max num levels
keys = list(match_dct.keys()) keys = list(match_dct.keys())
...@@ -1572,6 +1573,7 @@ def create_aeolus_clavrx_match_file(match_dct, filename, clvrx_params): ...@@ -1572,6 +1573,7 @@ def create_aeolus_clavrx_match_file(match_dct, filename, clvrx_params):
atimes.append(key) atimes.append(key)
elems.append(tup[0]) elems.append(tup[0])
lines.append(tup[1]) lines.append(tup[1])
fidxs.append(tup[2])
nlevs = prof.shape[0] nlevs = prof.shape[0]
if nlevs > max_num_alevels: if nlevs > max_num_alevels:
...@@ -1582,6 +1584,7 @@ def create_aeolus_clavrx_match_file(match_dct, filename, clvrx_params): ...@@ -1582,6 +1584,7 @@ def create_aeolus_clavrx_match_file(match_dct, filename, clvrx_params):
atimes = np.array(atimes) atimes = np.array(atimes)
elems = np.array(elems) elems = np.array(elems)
lines = np.array(lines) lines = np.array(lines)
fidxs = np.array(fidxs)
# Sample file to retrieve and copy variable attributes # Sample file to retrieve and copy variable attributes
rg_exmpl = Dataset('/home/rink/data/clavrx/clavrx_OR_ABI-L1b-RadF-M6C01_G16_s20192930000343.level2.nc', 'r') rg_exmpl = Dataset('/home/rink/data/clavrx/clavrx_OR_ABI-L1b-RadF-M6C01_G16_s20192930000343.level2.nc', 'r')
...@@ -1733,8 +1736,9 @@ def match_amvs_to_aeolus_fast(aeolus_dict, amv_files_path, amv_source='OPS', ban ...@@ -1733,8 +1736,9 @@ def match_amvs_to_aeolus_fast(aeolus_dict, amv_files_path, amv_source='OPS', ban
ds.close() ds.close()
for prof in profs: for prof in profs:
lat = prof[0, 0] profv = prof.values
lon = prof[0, 1] lat = profv[0, 0]
lon = profv[0, 1]
cc_prf, ll_prf = nav.earth_to_lc(lon, lat) cc_prf, ll_prf = nav.earth_to_lc(lon, lat)
if cc_prf is None or ll_prf is None: if cc_prf is None or ll_prf is None:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment