diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index d4232f343408c98dfcf3ed3bf4abe30c49168f79..1a378b6eb926c2451859d138d8af8383592edf8b 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -1686,6 +1686,7 @@ def match_amvs_to_aeolus_fast(aeolus_dict, amv_files_path, amv_source='OPS', ban
         amv_files = get_datasource(amv_files_path, amv_source, band=band)
     nav = amv_files.get_navigation()
     amv_params = amv_files.get_parameters()
+    all_params = [amv_files.lon_name, amv_files.lat_name, amv_files.elem_name, amv_files.line_name] + amv_params
     match_dict = {}
 
     keys = list(aeolus_dict.keys())
@@ -1748,7 +1749,9 @@ def match_amvs_to_aeolus_fast(aeolus_dict, amv_files_path, amv_source='OPS', ban
             # dist = haversine_np(lon, lat, amv_lons[in_box], amv_lats[in_box])
             param_nd = np.vstack(param_s)
             param_nd = param_nd[:, in_box]
-            match_dict[key].append((cc_prf, ll_prf, f_idx, prof, param_nd))
+            coords = {'num_params': amv_var_names}
+            amvs_da = xr.DataArray(param_nd, coords=coords, dims=['num_params', 'num_amvs'])
+            match_dict[key].append((cc_prf, ll_prf, f_idx, prof, amvs_da))
 
     return match_dict