diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 19eabafae0cb5fefb27e8255444503e559a201a0..fe4ef424cfddbb416bab0755256aef77c7da3e1c 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -524,7 +524,7 @@ def subset_by_time(match_file, dt_str_0, dt_str_1):
         accum_amv += n_amvs
         accum_prf += n_levs
 
-    prof_times_locs = np.hstack((times[time_idxs], lons[time_idxs], lats[time_idxs]))
+    prof_times_locs = np.column_stack([times[time_idxs], lons[time_idxs], lats[time_idxs]])
 
     return profs, amvs, prof_times_locs