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

snapshot...

parent 429a7da5
Branches
No related tags found
No related merge requests found
......@@ -402,14 +402,14 @@ def compute_and_add_geo_hgt(path_to_amvs, path_to_gfs, amv_source, band='14', ou
amv_lons = amvs_nd[:, amv_lon_idx]
amv_lats = amvs_nd[:, amv_lat_idx]
# temp_prof_s, gfs_press = get_temp_prof_s(gfs_ds, nom_time, amv_lons, amv_lats)
temp_prof_s, gfs_press, sfc_press, sfc_temp, sfc_hgt, = get_temp_prof_s_intrp(nom_time, amv_lons, amv_lats)
temp_prof_s, gfs_press = get_temp_prof_s(gfs_ds, nom_time, amv_lons, amv_lats)
# temp_prof_s, gfs_press, sfc_press, sfc_temp, sfc_hgt, = get_temp_prof_s_intrp(nom_time, amv_lons, amv_lats)
alt_f = []
for k in range(num_amvs):
amv_press = amvs_nd[k, amv_pres_idx]
#alt = pressure_to_altitude(amv_press, None, gfs_press, temp_prof_s[k, :])
alt = pressure_to_altitude(amv_press, None, gfs_press, temp_prof_s[k, :], sfc_pres=sfc_press[k], sfc_temp=sfc_temp[k], sfc_elev=sfc_hgt[k])
alt = pressure_to_altitude(amv_press, None, gfs_press, temp_prof_s[k, :])
#alt = pressure_to_altitude(amv_press, None, gfs_press, temp_prof_s[k, :], sfc_pres=sfc_press[k], sfc_temp=sfc_temp[k], sfc_elev=sfc_hgt[k])
alt_s.append(alt)
alt_f.append(alt)
prs_s.append(amv_press)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment