在Fortran中读取矩阵

Question

我尝试阅读并重写一个只有一系列特定列的矩阵。例如，我有一个11 x n矩阵，我想在2 x n矩阵中重写它，只有第2列和第3列。

program fluxXYlambda
implicit none

INTEGER*8::n,i
REAL*8,ALLOCATABLE::dati(:,:)

open(10,file='test.txt')
n=0

do
    read(10,*,end=100)
    n=n+1
end do
100 continue
rewind(10)

allocate(dati(11,n))

open(12,file='coordXY.txt')

do i=1,n
    read(10,*) dati(4,i),dati(5,i)
    write(12,*) dati(4,i),dati(5,i)
end do

end program flussiXYlambda

文件coordXY.txt已创建，但为空。

Answer 1

最后，由于@francescalus，我解决了这个问题。

命令read(10,*) dati(1:5,i)读取矩阵dati直到第5列，然后我可以保存它的特定列。非常感谢！