diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 69d1a92d3f3d8e8fe8bc25d451ff41db91b0b987..e3867acad58ae1676e2e7b92e72ace3c86042b29 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -2069,7 +2069,7 @@ def create_bestfit_file(filename, match_dct, raob_dct, gfs_at_raob_dct, bf_dct,
     rootgrp.createDimension('time', size=1)
 
     time = rootgrp.createVariable('time', 'f4', ['time'])
-    time.units = 'seconds since 1970-01-1 00:00:00'
+    time.units = 'seconds since 1970-01-01 00:00:00'
     time.calendar = 'standard'
     #amv_file_names = rootgrp.createVariable('amv_file_names', str, ['num_amv_files'])
 
@@ -2104,6 +2104,9 @@ def create_bestfit_file(filename, match_dct, raob_dct, gfs_at_raob_dct, bf_dct,
     amv_line.long_name = 'FGF y coordinate Full Disk GOES16'
     nc4_vars.append(amv_line)
 
+    dist_to_raob = rootgrp.createVariable('dist_to_raob', 'f4', ['time', 'num_raob_profs', 'max_num_amvs'])
+    dist_to_raob.units = 'km'
+
     var_s = rg_exmpl.variables
     for pidx, param in enumerate(amv_params):
         v = var_s[param]
@@ -2187,6 +2190,8 @@ def create_bestfit_file(filename, match_dct, raob_dct, gfs_at_raob_dct, bf_dct,
         var_v_g[0, idx, 0:num_amvs] = bf_gfs[:, 1]
         var_p_g[0, idx, 0:num_amvs] = bf_gfs[:, 2]
 
+        dist_to_raob[0, idx, 0:num_amvs] = haversine_np(key.lon, key.lat, amvs[:, 0], amvs[:, 1])
+
     rg_exmpl.close()
     rootgrp.close()