Fortran 代码中与 EOF 命令相关的错误

Error related to EOF command in fortran code

我是 Fortran 的初学者,正在尝试使用 gfortran 编译 Fixed-Term Fortran 代码。我有一堆错误,我可以修复它们。但是,有一个与 "EOF" 相关的错误,我无法解决。有什么办法可以解决这个问题吗? ("EOF" 两行分别是第 40 行和第 121 行。)

 37         OPEN(4,FILE="ABCE.Pn")
 38 
 39         OPEN(5,FILE="../sta.txt")
 40         DO WHILE (.not.EOF(5))
 41                 N=N+1
 42                 READ(5,*)STA(N)%COD,STA(N)%NAME,STA(N)%LAT,
 43      $          STA(N)%LON,STA(N)%H
 44         ENDDO
 45         NSTA=N
 46         CLOSE(5)`

 ......

121         DO WHILE (.not.EOF(1))
122                 READ(1,'(A60)',ERR=999) TIT
123 C               IF(IYEAR.GE.2008.OR. 
      (IYEAR.EQ.2007.AND.MONTH.GE.11)) 
124 C     $         TIT=TIT(2:60)
125                 IF(TIT(1:60).EQ.'')THEN         ! NEW EARTHQUAKE`

错误:

      DO WHILE (.not.EOF(5))                                            
                     1
Error: Operand of .not. operator at (1) is REAL(4)
ReadP2Pn.for:121.21:

      DO WHILE (.not.EOF(1))                                            
                     1
Error: Operand of .not. operator at (1) is REAL(4)

EOF(5) 是 non-standard。您应该在读取语句中检查 EOF(遗憾的是它看起来像一个 goto):

 40         DO WHILE (.true.)
 41                 N=N+1
 42                 READ(5,*,end=990)STA(N)%COD,STA(N)%NAME,STA(N)%LAT,
 43      $          STA(N)%LON,STA(N)%H
 44         ENDDO
 45    990  NSTA=N