我正在将Fortran 90代码转换为python。Fortran代码包含多维数组,例如:
integer fvpair(3,12,2)
integer jpair(3,2)
real*8 aread(3,12)
我不知道下面的numpy数组是否正确,假设初始化为零:
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)]])
如果要保留原始的Fortran阵列存储顺序(主要列),请不要忘记传递该order='F'
标志!
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')