diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index df838a49cbe6bdafaef7be31ccc00fa66d5bec88..6a28b4dc69c47a1ffa516f873529d4332719415d 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -2169,6 +2169,21 @@ def create_bestfit_file(filename, match_dct, raob_dct, raob_ts, gfs_at_raob_dct,
 
     var_flg_g = rootgrp.createVariable('bf_gfs_flag', 'i4', ['time', 'num_raob_profs', 'max_num_amvs'])
 
+    var_u_cp = rootgrp.createVariable('closest_press_u', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    var_u_cp.standard_name = 'eastward_wind'
+    var_u_cp.long_name = 'wind at closest raob pressure level to amv within +/- 60mb'
+    var_u_cp.units = 'm s-1'
+
+    var_v_cp = rootgrp.createVariable('closest_press_v', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    var_v_cp.standard_name = 'northward_wind'
+    var_v_cp.long_name = 'wind at closest raob pressure level to amv within +/- 60mb'
+    var_v_cp.units = 'm s-1'
+
+    var_p_cp = rootgrp.createVariable('closest_press', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    var_p_cp.standard_name = 'air_pressure'
+    var_p_cp.long_name = 'closest raob pressure level to amv within +/- 60mb'
+    var_p_cp.units = 'hPa'
+
     # Write data to file  ---------------------
     time[0] = raob_ts
     prf_lon[0, :] = alons
@@ -2196,6 +2211,11 @@ def create_bestfit_file(filename, match_dct, raob_dct, raob_ts, gfs_at_raob_dct,
         var_v[0, idx, 0:num_amvs] = bf_raob[:, 1]
         var_p[0, idx, 0:num_amvs] = bf_raob[:, 2]
 
+        cp_raob = bf_tup[2]
+        var_u_cp[0, idx, 0:num_amvs] = cp_raob[:, 0]
+        var_v_cp[0, idx, 0:num_amvs] = cp_raob[:, 1]
+        var_p_cp[0, idx, 0:num_amvs] = cp_raob[:, 2]
+
         bf_gfs = bf_tup[3]
         var_u_g[0, idx, 0:num_amvs] = bf_gfs[:, 0]
         var_v_g[0, idx, 0:num_amvs] = bf_gfs[:, 1]