diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 68099c36527eaa7e83bd048cc793d36bc533e749..4e5bb2d41c17b7f44aaaa865b7b4a5a0e8e3dff3 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -431,6 +431,7 @@ def analyze2(raob_to_amv_dct, raob_dct, gfs_filename=None):
     amvs = np.transpose(amvs, axes=[1, 0])
     bfs = np.stack(bf_list, axis=0)
     raob_match = np.stack(raob_match_list, axis=0)
+    bfs_gfs = np.stack(bf_gfs_list, axis=0)
 
     good_amvs = amvs
     num_good = good_amvs.shape[0]