diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py
index a99d92d3d14697d011c41f3d4865cd36bb9db7ba..c4674f93360dd9ddb8b50dd909bb1df8d9afb982 100644
--- a/modules/amv/intercompare.py
+++ b/modules/amv/intercompare.py
@@ -579,7 +579,7 @@ def run_best_fit_gfs_all(amvs, gfs_filename):
 
 
 def run_best_fit_gfs(amvs, gfs_filename, amv_lat_idx=amv_lat_idx, amv_lon_idx=amv_lon_idx, amv_prs_idx=amv_prs_idx,
-                     amv_spd_idx=amv_spd_idx, amv_dir_idx=amv_dir_idx):
+                     amv_spd_idx=amv_spd_idx, amv_dir_idx=amv_dir_idx, method='nearest'):
     xr_dataset = xr.open_dataset(gfs_filename)
     num_amvs = amvs.shape[0]
     bestfits = []
@@ -588,7 +588,7 @@ def run_best_fit_gfs(amvs, gfs_filename, amv_lat_idx=amv_lat_idx, amv_lon_idx=am
     gfs_press = gfs_press.data
     gfs_press = gfs_press[::-1]
 
-    uv_wind = get_vert_profile_s(xr_dataset, ['u-wind', 'v-wind'], amvs[:, amv_lon_idx], amvs[:, amv_lat_idx], method='nearest')
+    uv_wind = get_vert_profile_s(xr_dataset, ['u-wind', 'v-wind'], amvs[:, amv_lon_idx], amvs[:, amv_lat_idx], method=method)
     uv_wind = uv_wind.data
     wspd, wdir = spd_dir_from_uv(uv_wind[0,:,:], uv_wind[1,:,:])
     wspd = wspd.magnitude