diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 183c9333787c581c063e222e24e5c83adde5c9aa..c367d115b20053f1f38f3ea78773c52b7eda8239 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -488,6 +488,44 @@ def run_best_fit(raob_to_amv_dct, raob_dct, gfs_filename=None):
 
     return out_dct
 
+def run_compare():
+
+    raob_dir = '/home/rink/data/raob/'
+    raob_files = ['raob_soundings20191117_0000.cdf',
+                  'raob_soundings20191118_0000.cdf',
+                  'raob_soundings20191119_0000.cdf',
+                  'raob_soundings20191120_0000.cdf',
+                  'raob_soundings20191121_0000.cdf',
+                  'raob_soundings20191122_0000.cdf',
+                  'raob_soundings20191123_0000.cdf',
+                  'raob_soundings20191124_0000.cdf',
+                  'raob_soundings20191125_0000.cdf',
+                  'raob_soundings20191126_0000.cdf',
+                  'raob_soundings20191127_0000.cdf',
+                  'raob_soundings20191128_0000.cdf',
+                  'raob_soundings20191129_0000.cdf',
+                  'raob_soundings20191130_0000.cdf']
+
+    gfs_dir = '/home/rink/data/gfs/'
+    gfs_files = ['gfs.19111612_F012.h5',
+                 'gfs.19111712_F012.h5',
+                 'gfs.19111812_F012.h5',
+                 'gfs.19111912_F012.h5',
+                 'gfs.19112012_F012.h5',
+                 'gfs.19112112_F012.h5',
+                 'gfs.19112212_F012.h5',
+                 'gfs.19112312_F012.h5',
+                 'gfs.19112412_F012.h5',
+                 'gfs.19112512_F012.h5',
+                 'gfs.19112612_F012.h5',
+                 'gfs.19112712_F012.h5',
+                 'gfs.19112812_F012.h5',
+                 'gfs.19112912_F012.h5']
+
+    for k, file in enumerate(raob_files):
+        raob_dct, ts = get_raob_dict_cdf(raob_dir+file)
+        compare_gfs_raob(raob_dct, gfs_dir+gfs_files[k])
+
 
 def compare_gfs_raob(raob_dct, gfs_filename):
     keys = list(raob_dct.keys())
@@ -545,7 +583,6 @@ def compare_gfs_raob(raob_dct, gfs_filename):
         r_4 = np.where(raob_prs == 100)[0]
 
         spd_intrp = np.interp(gfs_press, raob_prs, raob_v)
-        #print(raob_v[r_0] - gfs_v[key_idx, g_0], raob_v[r_1] - gfs_v[key_idx, g_1], raob_v[r_3] - gfs_v[key_idx, g_3])
 
         if len(r_0) > 0:
             r_0 = r_0[0]
@@ -580,6 +617,7 @@ def compare_gfs_raob(raob_dct, gfs_filename):
     print(np.average(np.abs(lvl_2)))
     print(np.average(np.abs(lvl_3)))
     print(np.average(np.abs(lvl_4)))
+    print('**************************************')
 
     return None