我对Python还不太熟悉,我从头开始用Python编写可视化代码(以避免使用昂贵的专有程序,如IDL)。到目前为止,我一直使用IDL和gnuplot。我想做的是:
我使用fortran将二维数组写入未格式化的直接访问文件,我希望能够用python读取这些文件。下面给出了准确的测试代码。实际的代码是一个巨大的并行代码,但数据输出几乎是完全相同的格式。program binary_out
implicit none
integer :: i,j,t,rec_array
double precision, dimension(100,100) :: fn
double precision, parameter :: p=2*3.1415929
INQUIRE(IOLENGTH=rec_array) fn
open(unit=10,file='test',status='new',form='unformatted',access='direct',recl=rec_array)
fn=0
write(10,rec=1) fn
do t=1,3
do i=1,100
do j=1,100
fn(i,j)=sin(i*p*t/100)*cos(j*p*t/100)
enddo
enddo
write(10,rec=t+1) fn
enddo
close(10)
end program binary_out
对于t=1,程序应该给我零,对于t值的增