scripts/gray: fix read_eqdsk bug
The poloidal flux matrix is written in fortran order!
This commit is contained in:
parent
eacfa2efe2
commit
9442356797
@ -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
|
||||
|
Loading…
Reference in New Issue
Block a user