[Numpy-discussion] example reading binary Fortran file
Sturla Molden
sturla@molden...
Fri Jan 30 10:23:45 CST 2009
On 1/30/2009 5:03 PM, David Froger wrote:
> I think it will be a good idea to put the Fortran writting-arrays code
> and the Python reading-array script in the cookbook and maybe a page to
> help people comming from Fortran to start with Python ?
If you want to be completely safe, read the file in Fortran, then send
it as an array to Python (use f2py). Aside from that, assuming your
compiler only writes the raw bytes in Fortran order to the file:
If you have
real(4),dimension(nx,ny):: ux
in Fortran, and writes ux to disk. It could be retrieved like this in NumPy:
ux = np.fromfile(nx*ny, dtype=np.float32).view((nx,ny), order='F')
assuming real(kind=4) is equivalent to np.float32.
Sturla Molden
More information about the Numpy-discussion
mailing list