diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index ffe4fa0d0c308ebcae6684d8176c06833c4d9613..05e2ff61143f1b8c638b07e74f71f7e85dfd64f9 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -2117,13 +2117,19 @@ def create_bestfit_file(filename, match_dct, raob_dct, gfs_at_raob_dct, bf_dct,
 
     # Bestfit variables ------------------
     var_u = rootgrp.createVariable('bf_raob_u', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    var_u.units = 'm s-1'
     var_v = rootgrp.createVariable('bf_raob_v', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    var_v.units = 'm s-1'
     var_p = rootgrp.createVariable('bf_raob_press', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    var_p.units = 'millibar'
     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_u_g.units = 'm s-1'
     var_v_g = rootgrp.createVariable('bf_gfs_v', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    var_v_g.units = 'm s-1'
     var_p_g = rootgrp.createVariable('bf_gfs_press', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    var_p_g.units = 'millibar'
     var_flg_g = rootgrp.createVariable('bf_gfs_flag', 'i4', ['time', 'num_raob_profs', 'max_num_amvs'])
 
     # Write data to file  ---------------------