diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 1c08d43e60881671ea868e2af634b9dfc4fde90e..df7832e6728274e31f00d50af319556f1f97b26b 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -473,15 +473,9 @@ def run_best_fit_all():
         out_list.append(dct)
         prd_dct = get_product_at_locs(m_d, ts, prd_files)
 
-        #analyze2(amvs_list, bf_list, raob_match_list, bf_gfs_list, amv_prod_list)
-
-        # amvs = get_amvs(amv_files, ts)
-        # amvs_list.append(amvs)
-        # bfs = run_best_fit_gfs(amvs, gfs_dir+gfs_files[k], amv_lat_idx=0, amv_lon_idx=1, amv_prs_idx=4, amv_spd_idx=5, amv_dir_idx=6)
-        # bf_list.append(bfs)
-
-    # bin_ranges, bin_pres, bin_spd, bin_dir = analyze2(amvs_list, bf_list, raob_match_list, bf_gfs_list, amv_prod_list)
-    # bin_ranges, bin_pres, bin_spd, bin_dir = compare_amvs_bestfit_all(amvs_list, bf_list, bin_size=100)
+        amvs = get_amvs(amv_files, ts)
+        bfs = run_best_fit_gfs(amvs, gfs_dir+gfs_files[k], amv_lat_idx=0, amv_lon_idx=1, amv_prs_idx=4, amv_spd_idx=5, amv_dir_idx=6)
+        out_list.append((amvs, np.array(bfs)))
 
     return out_list