diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 1ae8c11b53aa766a888da808d19533f83d304d1c..6046651db5aa655620f5eeed0f0719441dafff66 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -2105,6 +2105,11 @@ def create_bestfit_file(filename, match_dct, raob_dct, gfs_at_raob_dct, bf_dct,
     var_p = rootgrp.createVariable('bf_raob_press', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
     var_flg = rootgrp.createVariable('bf_raob_flag', 'i4', ['time', 'num_raob_profs', 'max_num_amvs'])
 
+    var_u_g = rootgrp.createVariable('bf_gfs_u', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    var_v_g = rootgrp.createVariable('bf_gfs_v', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    var_p_g = rootgrp.createVariable('bf_gfs_press', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    var_flg_g = rootgrp.createVariable('bf_gfs_flag', 'i4', ['time', 'num_raob_profs', 'max_num_amvs'])
+
     # Write data to file  ---------------------
     prf_lon[0, :] = alons
     prf_lat[0, :] = alats
@@ -2131,6 +2136,9 @@ def create_bestfit_file(filename, match_dct, raob_dct, gfs_at_raob_dct, bf_dct,
         var_p[0, idx, 0:num_amvs] = bf_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]
+        var_p_g[0, idx, 0:num_amvs] = bf_gfs[:, 2]
 
     rg_exmpl.close()
     rootgrp.close()