snapshot...

modules/amv/caliop_clavrx_amv.py

@@ -404,13 +404,13 @@ def compute_and_add_geo_hgt(path_to_amvs, path_to_gfs, amv_source, band='14', ou
         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 = get_temp_prof_s_intrp(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])
+            alt = pressure_to_altitude(amv_press, None, gfs_press, temp_prof_s[k, :], sfc_pres=sfc_press[k], sfc_temp=sfc_temp[k], sfc_hgt=sfc_hgt[k])
             alt_s.append(alt)
             alt_f.append(alt)
             prs_s.append(amv_press)