Read(u=UNIT,*) 在 fortran 中只读取一半文件

Question

所以我试图创建一个 Fortran 子例程，从我在另一个子例程中创建的文件中读取大小为 125x125x125 的矩阵（或我猜的张量），但由于某种原因它不起作用。我已经使用几乎相同的程序成功完成了它，但对于不同大小的矩阵（70x200x70），并且出于某种原因，当我更改数组的大小和其他一些东西以适应新数据时，程序只读取一半数据文件，并抛出一个错误，显示 At line 133 of file Busqueda.f90 (unit = 2, file = 'test2.txt') Fortran Runtime error: End of file。所以我的工作代码是：

Subroutine DataSorting(datasort,indexsort,datosden,datostemp)
use globals
implicit None
real*8, dimension(70,200,70)              :: datosdentemp
real*8, dimension(70,200,70), intent(out) :: datostemp, datosden
real*8, dimension(980000), intent(out)    :: datasort
integer,dimension(3,980000),intent(out) :: indexsort
integer, dimension(3) :: Result
integer               :: i
rewind(2)
read(2,*) datosden
datosdentemp = datosden
rewind(3)
read(3,*) datostemp
do i = 1, 10
    Result = MAXLOC(datosdentemp)
    datasort(i) = datosden(Result(1),Result(2),Result(3))
    indexsort(1,i) = Result(1)
    indexsort(2,i) = Result(2)
    indexsort(3,i) = Result(3)
    datosdentemp(Result(1),Result(2),Result(3)) = 0
end do
End Subroutine DataSorting

不起作用的代码是：

Subroutine DataSorting(datasort,indexsort,datosden,datostemp)
use globals
implicit None
real*8, dimension(125,125,125)              :: datosdentemp
real*8, dimension(125,125,125), intent(out) :: datostemp, datosden
real*8, dimension(1953125), intent(out)    :: datasort
integer,dimension(3,1953125),intent(out) :: indexsort
integer, dimension(3) :: Result
integer               :: i
rewind(2)
read(2,*) datosden
datosdentemp = datosden
rewind(3)
read(3,*) datostemp
do i = 1, 10
    Result = MAXLOC(datosdentemp)
    datasort(i) = datosden(Result(1),Result(2),Result(3))
    indexsort(1,i) = Result(1)
    indexsort(2,i) = Result(2)
    indexsort(3,i) = Result(3)
    datosdentemp(Result(1),Result(2),Result(3)) = 0
end do
End Subroutine DataSorting

在非工作程序的情况下，我读取数据的文件中的数据仅排列为 125 列和 15625 行，而在非工作程序的情况下，数据排列为 70 列和工作行有 14000 行，文件就是这样，我的意思是在这两种情况下它们实际上都是这样构造的。我对文件 2 和 3 使用了临时文件，并将它们更改为“.txt”文件，以查看输入文件是否有问题，但没有。然后我构建一个新文件并使用 do-loop 进行读取，看看这是否是问题所在，但不是，但至少这帮助我意识到代码正在读取 15625 的一半，但我不知道知道为什么会发生这种情况，我是 Fortran 新手，我的论文导师告诉我用 Fortran 编写这段代码，因为必须处理的数据量会让我的 Julia 代码变慢，所以说实话，任何帮助都会很大赞赏。