Fortran：读取行数未知的文件

Question

我正在尝试打开行数未知的 ASCII 文件（但每行中的条目数固定 - 请参阅下面的示例文件）。它转到文件末尾，但之后失败。

选项 1：使用 iostat 语句：使用“Use, intrinsic :: iso_fortran_env, Only : iostat_end”或简单地将 io 视为整数时，我得到相同的错误。

选项 2：使用“end= ..”选项。这是我之前使用的（F77）。

这两种方法原则上都应该有效，但两者都没有。任何解决此问题的帮助将不胜感激。

谢谢， 皮纳基。

计划：

!=======================================!
program read
!=======================================!
!  Use, intrinsic :: iso_fortran_env, Only : iostat_end
  implicit none

  integer :: i,io,n     

  n=0
  open(10,file='a.dat',status='old',action='read')

!=================================!
! option 1
!=================================!  
  do 
     read(10,*,iostat=io)
     write(*,*)io
!     if (io.eq.iostat_end) exit
     if (io.ne.0) exit
     n=n+1
     write(*,*)'n=',n
  enddo
!=================================!

!=================================!
! option 2
!=================================!
  do 10 i=1,1000000
     read(10,*,end=10)
     n=n+1
     write(*,*)n
10   continue
!================================!

  close(10)
  write(*,*)'n=',n
  
end program read
!======================!

使用“gfortran --std=f2003 -o read.out read.f90”编译

错误信息：

===============================================

dyld: lazy symbol binding failed: Symbol not found: ___emutls_get_address
  Referenced from: /usr/local/opt/gcc/lib/gcc/11/libgfortran.5.dylib
  Expected in: /usr/lib/libSystem.B.dylib

dyld: Symbol not found: ___emutls_get_address
  Referenced from: /usr/local/opt/gcc/lib/gcc/11/libgfortran.5.dylib
  Expected in: /usr/lib/libSystem.B.dylib

程序收到信号SIGABRT：进程中止信号。

此错误的回溯： #0 0x105364f8e #1 0x10536419d #2 0x7fff6eed05fc zsh: 中止 ./read.out

================================================ ======

我正在阅读的文件：

================================================ ======

 1    -1.1559859375    2.0399371094    0.1686166667    0.8242152778
  2    -1.1618015625    2.1375250000    0.1765231481    0.8105046296
  3    -1.1696710937    2.2325417969    0.1860513889    0.7936782407
  1    -1.1730312500    2.3271382813    0.1975254630    0.7718773148
  2    -1.1767945313    2.3942726563    0.2113162500    0.7446933333
  3    -1.1694437500    2.4738000000    0.2281966667    0.7099266667
  1    -1.1566140625    2.5312164063    0.2494466667    0.6636841667
  2    -1.1293765625    2.5746707031    0.2747766667    0.6066154167
  3    -1.0836390625    2.5938144531    0.3026616667    0.5403733333
  1    -1.0380632812    2.5721433594    0.3302462121    0.4727159091

================================================ =========================

Answer 1

我测试了您的代码，第一个选项有效。如果您从do 10 i=1,1000000 中删除标签10 并在10 continue 之前添加end do，您的第二个选项也将起作用。您还必须在第二个选项的 do 循环之前添加 rewind(10) 或 close(10); open(10,file='a.dat',status='old',action='read')。但更有趣的是，在 Fortran 2008 中使用检查文件结尾错误代码出现的 is_iostat_end() 内部函数很容易解决这个问题。这是一个现代的实现，

program read
    implicit none
    integer :: i, n, iostat, fileUnit
    n = 0
    open(newunit = fileUnit, file = 'a.dat', status = 'old', action = 'read')
    do
         read(fileUnit,*,iostat = iostat)
         if (is_iostat_end(iostat)) exit
         write(*,*) iostat
         n = n + 1
         write(*,*)'n = ', n
    end do
    close(fileUnit)
end program read

使用以下命令编译并运行它，

gfortran --std=f2008 -o read.out read.f90
./read.out

尽管问题很小，但您通过提供完整代码、数据文件、编译代码的方法以及遇到的错误提出了一个非常好的问题。它值得一票。