Converting Fortran array to numpy array

Question

I am in the process of converting a Fortran 90 code to python. The Fortran code includes multidimensional array like:

integer fvpair(3,12,2)
integer jpair(3,2)
real*8 aread(3,12)

I wonder if the following numpy array are correct equivalent, assuming zero initialization:

fvpair = np.array([[np.zeros(3)],[np.zeros(12)],[np.zeros(2)]])
jpair = np.array([[np.zeros(3)],[np.zeros(2)]])
aread = np.array([[np.zeros(3)],[np.zeros(12)]])

You probably want fvpair = np.zeros((3,12,2), dtype=int).

tom10
– tom10

2020-11-28 03:50:54 +00:00
Commented Nov 28, 2020 at 3:50 — tom10
– tom10, Commented Nov 28, 2020 at 3:50

DYZ · Accepted Answer · 2020-11-28 05:11:29Z

3

If you want to preserve the original Fortran array storage order (column-major), do not forget to pass the order='F' flag!

fvpair = np.zeros((3,12,2), dtype=int, order='F')
jpair = np.zeros((3,2), dtype=int, order='F')
aread = np.zeros((3,2), dtype=float64, order='F')

answered Nov 28, 2020 at 5:11

DYZ

57.3k10 gold badges73 silver badges101 bronze badges

Sign up to request clarification or add additional context in comments.

Collectives™ on Stack Overflow

Converting Fortran array to numpy array

1 Answer 1

Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

Comments

Your Answer

Sign up or log in

Post as a guest

Related