diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 1b6a758bea73e75bed65e37051a13db93917aab9..07a0dcf4aa061bdbbedb45fff47572e8ffd1f465 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -2062,11 +2062,11 @@ def create_bestfit_file(filename, match_dct, raob_dct, gfs_at_raob_dct, bf_dct,
     # the top level group for the output file
     rootgrp = Dataset(filename, 'w', format='NETCDF4')
 
-    dim_aparams = rootgrp.createDimension('num_raob_params', size=num_aparams)
-    dim_alevs = rootgrp.createDimension('max_num_raob_levels', size=max_num_alevels)
-    dim_amvs = rootgrp.createDimension('max_num_amvs', size=max_num_amvs)
-    dim_num_aeolus_prof = rootgrp.createDimension('num_raob_profs', size=num_aprofs)
-    dim_time = rootgrp.createDimension('time', size=1)
+    rootgrp.createDimension('num_raob_params', size=num_aparams)
+    rootgrp.createDimension('max_num_raob_levels', size=max_num_alevels)
+    rootgrp.createDimension('max_num_amvs', size=max_num_amvs)
+    rootgrp.createDimension('num_raob_profs', size=num_aprofs)
+    rootgrp.createDimension('time', size=1)
 
     time = rootgrp.createVariable('time', 'f4', ['time'])
     time.units = 'seconds since 1970-01-1 00:00:00'
@@ -2084,7 +2084,7 @@ def create_bestfit_file(filename, match_dct, raob_dct, gfs_at_raob_dct, bf_dct,
     prf_spd = rootgrp.createVariable('raob_spd', 'f4', ['time', 'num_raob_profs', 'max_num_raob_levels'])
     prf_spd.units = 'm s-1'
     prf_lvl = rootgrp.createVariable('raob_levels', 'f4', ['time', 'num_raob_profs', 'max_num_raob_levels'])
-    prf_lvl.units = 'millibar'
+    prf_lvl.units = 'hPa'
 
     # ----- Product variables ----------------
     nc4_vars = []
@@ -2096,8 +2096,13 @@ def create_bestfit_file(filename, match_dct, raob_dct, gfs_at_raob_dct, bf_dct,
     amv_lat.units = 'degrees_north'
     nc4_vars.append(amv_lat)
 
-    nc4_vars.append(rootgrp.createVariable('amv_elem', 'i4', ['time', 'num_raob_profs', 'max_num_amvs']))
-    nc4_vars.append(rootgrp.createVariable('amv_line', 'i4', ['time', 'num_raob_profs', 'max_num_amvs']))
+    amv_elem = rootgrp.createVariable('amv_elem', 'i4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    amv_elem.long_name = 'FGF x coordinate Full Disk GOES16'
+    nc4_vars.append(amv_elem)
+
+    amv_line = rootgrp.createVariable('amv_line', 'i4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    amv_line.long_name = 'FGF y coordinate Full Disk GOES16'
+    nc4_vars.append(amv_line)
 
     var_s = rg_exmpl.variables
     for pidx, param in enumerate(amv_params):
@@ -2119,19 +2124,37 @@ 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.standard_name = 'eastward_wind'
+    var_u.long_name = 'wind at bestfit level with raob as background profile'
     var_u.units = 'm s-1'
+
     var_v = rootgrp.createVariable('bf_raob_v', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    var_v.standard_name = 'northward_wind'
+    var_v.long_name = 'wind at bestfit level with raob as background profile'
     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_p.standard_name = 'air_pressure'
+    var_p.long_name = 'pressure level of bestfit with raob as background profile'
+    var_p.units = 'hPa'
+
     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.standard_name = 'eastward_wind'
+    var_u_g.long_name = 'wind at bestfit level with raob as background profile'
     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.standard_name = 'northward_wind'
+    var_v_g.long_name = 'wind at bestfit level with raob as background profile'
     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_p_g.standard_name = 'air_pressure'
+    var_p_g.long_name = 'pressure level of bestfit with raob as background profile'
+    var_p_g.units = 'hPa'
+
     var_flg_g = rootgrp.createVariable('bf_gfs_flag', 'i4', ['time', 'num_raob_profs', 'max_num_amvs'])
 
     # Write data to file  ---------------------