diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index d739d785dd6e2e7653fbb6096a917d9e09f2674f..f090c3d7ea1a0ca42d3c2837e142fb0d5e6e6265 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -556,7 +556,7 @@ def subset_by_time(match_file, dt_str_0, dt_str_1, amv_var_names):
     mx_nlevs = np.max(n_levs_per_prof[time_idxs[0]:time_idxs[0]+n_times])
 
     amvs = np.zeros((n_times, mx_namvs, nvars))
-    profs = np.zeros((n_times, mx_nlevs, nvars))
+    profs = np.zeros((n_times, mx_nlevs, 4))
     amvs.fill(np.nan)
     profs.fill(np.nan)