diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index edd3174fcc44e7aaf69218525f221fc3aff8d811..1019f6870a88d6ddcf4ba856ac4a596280f4506b 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -123,7 +123,10 @@ def match_amvs_to_raobs(raob_dict, raob_time, amv_ds, time_window=10, filepath=N
         vld = good
 
     param_nd = np.vstack(param_s)
+    print(param_nd.shape)
     param_nd = param_nd[:, vld]
+    print(param_nd.shape)
+    print(filepath, '----------------')
     cc = param_nd[2, :]
     ll = param_nd[3, :]
 
@@ -555,6 +558,10 @@ def run_best_fit(raob_to_amv_dct, raob_dct, gfs_at_raob_dct):
     return out_dct
 
 
+# def analyze_bestfit_product_file(file):
+#     xr_ds = xr.load_dataset(file)
+
+
 def analyze2(amvs, bfs, raob_match, bfs_gfs, amv_prod):
 
     good_amvs = amvs