Прочитанный двоичный файл на Фортране на Python
Я пытаюсь прочитать вывод двоичного файла из кода Fortran ниже, но результаты не совпадают с выходным файлом.
Код Фортрана 77:
program test
implicit none
integer i,j,k,l
real*4 pcp(2,3,4)
open(10, file='pcp.bin', form='unformatted')
l = 0
do i=1,2
do j=1,2
do k=1,2
print*,k+l*2
pcp(i,j,k)=k+l*2
l = l + 1
enddo
enddo
enddo
do k=1,4
write(10)pcp(:,:,k)
enddo
close(10)
stop
end
Я пытаюсь использовать код Python ниже:
from scipy.io import FortranFile
f = FortranFile('pcp.bin', 'r')
a = f.read_reals(dtype=float)
print(a)