diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index e96e7d4a07094499a873028fa89b48ddc2e1d025..929388214a55a7d3860b5d654f4a135981c97966 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -419,6 +419,7 @@ def run_best_fit_all():
         raob_dct, ts = get_raob_dict_cdf(raob_dir+file, raob_dtsts[k])
         m_d = match_amvs_to_raobs(raob_dct, ts, amv_files=amv_files)
         amvs_list, bf_list, raob_match_list, bf_gfs_list = run_best_fit(m_d, raob_dct, gfs_dir+gfs_files[k])
+        analyze2(None, None, amvs_list, bf_list, raob_match_list, bf_gfs_list)
 
 
 def run_best_fit(raob_to_amv_dct, raob_dct, gfs_filename=None):