diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 4c1b1fae26452f6b77c9104df6a0083a17d2bfe7..c24515bc0b3bf481b02ea8e06b5cf11361a8281b 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -469,13 +469,13 @@ def run_best_fit_all():
     for k, file in enumerate(raob_files):
         raob_dct, ts = get_raob_dict_cdf(raob_dir+file)
         m_d = match_amvs_to_raobs(raob_dct, ts, amv_files=amv_files)
-        bf_dct = run_best_fit(m_d, raob_dct, gfs_dir+gfs_files[k])
+        #bf_dct = run_best_fit(m_d, raob_dct, gfs_dir+gfs_files[k])
         prd_dct = get_product_at_locs(m_d, ts, prd_files)
-        out_list.append((bf_dct, prd_dct))
+        #out_list.append((bf_dct, prd_dct))
 
-        # 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)))
+        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), prd_dct))
 
     return out_list
 
@@ -548,7 +548,7 @@ def get_product_at_locs1x1(raob_to_amv_dct, ts, files, filepath=None):
             for k in range(num_amvs):
                 aaa[vidx, k] = var[ll[k], cc[k]].data
 
-        m_dct[key] = aaa
+        m_dct[key] = np.transpose(aaa, axes=[1,0])
 
     ds.close()