diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py
index a2156bc009b70522df69fde6eca0d73af4334784..bd8edc3afb7322e8d8089be74f6c85729da0a835 100644
--- a/modules/amv/intercompare.py
+++ b/modules/amv/intercompare.py
@@ -418,7 +418,7 @@ def amv_raob_diff(amv_lat, amv_lon, amv_spd, amv_dir, amv_press, raob_dct):
     return spd_diffs, dir_diffs
 
 
-def run_best_fit(amvs, raob_dct, min_num_levs=50, test_a=True, test_b=True):
+def run_best_fit(amvs, raob_dct, min_num_levs=50):
     num_amvs = amvs.shape[0]
     bestfits = []
 
@@ -453,7 +453,7 @@ def run_best_fit(amvs, raob_dct, min_num_levs=50, test_a=True, test_b=True):
                 alat = amvs[j, amv_lat_idx]
                 alon = amvs[j, amv_lon_idx]
 
-                bf = best_fit(aspd, adir, aprs, alat, alon, rspd, rdir, rprs, test_a=test_a, test_b=test_b)
+                bf = best_fit(aspd, adir, aprs, alat, alon, rspd, rdir, rprs)
                 bf = (j, bf[0], bf[1], bf[2], bf[3])
                 bestfits.append(bf)
 
@@ -488,18 +488,18 @@ def run_amv_clavr_compare(amvs, clvr_filename):
     print('--------------')
 
 
-def run_best_fit_gfs_all(amvs, gfs_filename, test_a=True, test_b=True):
+def run_best_fit_gfs_all(amvs, gfs_filename):
     num_centers = len(amvs)
 
     bfs_s = []
     for k in range(num_centers):
-        bfs = run_best_fit_gfs(amvs[k], gfs_filename, test_a=test_a, test_b=test_b)
+        bfs = run_best_fit_gfs(amvs[k], gfs_filename)
         bfs_s.append(bfs)
 
     return bfs_s
 
 
-def run_best_fit_gfs(amvs, gfs_filename, test_a=True, test_b=True):
+def run_best_fit_gfs(amvs, gfs_filename):
     xr_dataset = xr.open_dataset(gfs_filename)
     num_amvs = amvs.shape[0]
     bestfits = []
@@ -523,7 +523,7 @@ def run_best_fit_gfs(amvs, gfs_filename, test_a=True, test_b=True):
         rspd = wspd[j]
         rdir = wdir[j]
 
-        bf = best_fit(aspd, adir, aprs, alat, alon, rspd, rdir, gfs_press, test_a=test_a, test_b=test_b)
+        bf = best_fit(aspd, adir, aprs, alat, alon, rspd, rdir, gfs_press)
         bf = (j, bf[0], bf[1], bf[2], bf[3])
         bestfits.append(bf)