diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index f94713e0a5e9ed8c7a2bd164f76754eb30e51e07..95e4377e8be5160e1165021179798529f9c5ff5e 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -511,7 +511,7 @@ def get_product_at_locs(raob_to_amv_dct, ts, files, filepath=None):
         for k in range(num_amvs):
             aaa[:, k] = param_nd[:, ll[k], cc[k]]
 
-        m_dct[key] = aaa
+        m_dct[key] = aaa.transpose(aaa, axes=(1, 0))
 
     ds.close()
 
@@ -1064,6 +1064,7 @@ def compare_amvs_bestfit_all(all_list, bin_size=200):
 
     amvs = np.concatenate(amvs_list)
     bfs = np.concatenate(bfs_list)
+    print(amvs.shape, bfs.shape)
     prd = np.concatenate(prd_list)
 
     print(amvs.shape, bfs.shape, prd.shape)