diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 3689822c0ea7c92fb19dea73bd8ab2aa7ef94c7c..dcb9dbdf8a095cb88c1e2a6b957e797badfc7e7b 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -112,7 +112,6 @@ def match_amvs_to_raobs(raob_dict, raob_time, amv_ds, time_window=10, filepath=N
             param_s.append(data.data[1])
         else:
             param_s.append(data.data)
-    param_s.append(np.arange(num_amvs))  # The absolute dataset index of amvs
 
     # filter
     qc_name = amv_ds.get_qc_params()
@@ -2138,9 +2137,6 @@ def create_bestfit_file(filename, match_dct, raob_dct, raob_ts, gfs_at_raob_dct,
     dist_to_raob.units = 'km'
     dist_to_raob[:,] = np.nan
 
-    amv_source_index = rootgrp.createVariable('amv_source_index', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
-    amv_source_index[:,] = -9999
-
     var_s = rg_exmpl.variables
     for pidx, param in enumerate(amv_params):
         v = var_s[param]
@@ -2258,12 +2254,6 @@ def create_bestfit_file(filename, match_dct, raob_dct, raob_ts, gfs_at_raob_dct,
 
         dist_to_raob[0, idx, 0:num_amvs] = haversine_np(key.lon, key.lat, amvs[:, 0], amvs[:, 1])
 
-        param_s = match_dct.get(key)
-        idx_idx = param_s.shape[0] - 1
-        amv_idxs = param_s[idx_idx, ]
-
-        amv_source_index[0, idx, 0:num_amvs] = amv_idxs[:]
-
     rg_exmpl.close()
     rootgrp.close()