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

snaphot...

parent 79628868
Branches
No related tags found
No related merge requests found
......@@ -52,6 +52,8 @@ def match_calipso_clavrx_to_amvs(calipso_clavrx_ds, calipso_clavrx_file, amv_fil
nom_time = calipso_clavrx_ds.get_datetime(calipso_clavrx_file).timestamp()
calipso_clavrx_nc4 = Dataset(calipso_clavrx_file)
print('start processing: ', calipso_clavrx_file)
amv_fname, ftime, f_idx = amv_files.get_file(nom_time)
if f_idx is None:
return None
......@@ -64,7 +66,6 @@ def match_calipso_clavrx_to_amvs(calipso_clavrx_ds, calipso_clavrx_file, amv_fil
gfs_xr = xr.open_dataset(gfs_fname)
gfs_press = gfs_xr['pressure levels']
gfs_press = gfs_press.values
#gfs_press = gfs_press[::-1]
match_dict[nom_time] = []
......@@ -132,11 +133,10 @@ def match_calipso_clavrx_to_amvs(calipso_clavrx_ds, calipso_clavrx_file, amv_fil
temp_prof = get_vert_profile_s(gfs_xr, ['temperature'], lons, lats, method='nearest')
temp_prof = temp_prof.values
temp_prof_s = temp_prof[0, :, :]
#temp_prof_s = temp_prof_s[:, ::-1]
for k in range(num_amvs):
alt = pressure_to_altitude(param_nd[4, k], None, gfs_press, temp_prof_s[k, :])
print(param_nd[4, k], alt)
# print(param_nd[4, k], alt)
return match_dict
......
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment