diff --git a/bfps/DNS.py b/bfps/DNS.py
index 6a81508021135470a516dbd1f0f73916d0e4a607..2e18b51c794639ffa0cc070a7a3716c5496be2d6 100644
--- a/bfps/DNS.py
+++ b/bfps/DNS.py
@@ -879,8 +879,8 @@ class DNS(_code):
         # the case of increasing/decreasing by the same factor in all directions.
         # in principle we could write something more generic, but i'm not sure
         # how complicated that would be
-        dst_shape = (self.parameters['nz'],
-                     self.parameters['ny'],
+        dst_shape = (self.parameters['ny'],
+                     self.parameters['nz'],
                      (self.parameters['nx']+2) // 2,
                      3)
         src_file = h5py.File(src_file_name, 'r')