diff --git a/scripts/gray/__init__.py b/scripts/gray/__init__.py
index 91ddaa5..8922480 100644
--- a/scripts/gray/__init__.py
+++ b/scripts/gray/__init__.py
@@ -58,7 +58,7 @@ def read_eqdsk(filepath: Path) -> Eqdsk:
     skip(file, nr * 4)
 
     # read flux and safety factor
-    psi = read(file, nr * nz).reshape(nr, nz)
+    psi = read(file, nr * nz).reshape((nr, nz), order='F').T
     q = read(file, nr)
 
     # normalise flux