diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 6af7a69441fa9024fb6f700a499fbb3f17464354..1b136861d5030addcb75dbaef3547c8d345dbfa0 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -474,9 +474,10 @@ def run_best_fit_all():
         amvs_list, bf_list, raob_match_list, bf_gfs_list = run_best_fit(m_d, raob_dct, gfs_dir+gfs_files[k],
                                     amvs_list=amvs_list, bf_list=bf_list, raob_match_list=raob_match_list, bf_gfs_list=bf_gfs_list)
         # analyze2(amvs_list, bf_list, raob_match_list, bf_gfs_list)
-        prd_dct = get_product_at_locs(m_d, ts, prd_files, amv_prod_list=amv_prod_list)
-        # bfs_dct = run_best_fit2(m_d, raob_dct, gfs_dir+gfs_files[k])
-    analyze2(amvs_list, bf_list, raob_match_list, bf_gfs_list, amv_prod_list)
+        # prd_dct = get_product_at_locs(m_d, ts, prd_files, amv_prod_list=amv_prod_list)
+    bin_ranges, bin_pres, bin_spd, bin_dir = analyze2(amvs_list, bf_list, raob_match_list, bf_gfs_list, amv_prod_list)
+
+    return bin_ranges, bin_pres, bin_spd, bin_dir
 
 
 def get_product_at_locs(raob_to_amv_dct, ts, files, amv_prod_list=[], filepath=None):