读写错误，请帮忙。

VBS_CS

[Fortran] 纯文本查看 复制代码

    PROGRAM TEST_WRNF
    IMPLICIT NONE
    Real, Dimension(9,62)  ::EPa
    CALL EllipseP(EPa)
    END

         Subroutine EllipseP(EPa)
         IMPLICIT NONE
         Real,    Intent(OUT), Dimension(9,62) :: EPa
         Integer :: I, J
         EPa=0.
         OPEN (103 ,FILE='CEP.txt', STATUS='UNKNOWN') 
         Call Notes103
         Do I=1,62 
             Read(103,*) (EPa(J,I), J=1,9)
         Enddo
         CLOSE(103)
         End  Subroutine  EllipseP

        Subroutine Notes103
        CHARACTER*1 AX(80)
  11    Read(103,'(80A1)')(AX(I),I=1,80)
        Do I=1,80
          If(AX(I)==']') Return
        Enddo
        GOTO 11
        END

以下为文本信息：
  [ρ/Rw  β       m     n       mn     a/b    Ae      Be       AeBe]
   0.05   12.84   5.22   0.352   1.84   14.83  0.290   0.0195   5.66
   0.06   14.07   5.00   0.361   1.81   13.85  0.295   0.0213   6.27
   0.07   15.20   4.79   0.369   1.77   12.98  0.297   0.0229   6.81
   0.08   16.26   4.57   0.378   1.73   12.09  0.297   0.0245   7.27
   0.09   17.25   4.36   0.386   1.68   11.30  0.294   0.0261   7.67

左耳北上南

第14行，你这里是有5行数据可以读入，不能写成62行。这里报错。
还有二楼的程序比你的简洁而且适用性广。

fcode

你把很多问题想复杂了。

[Fortran] 纯文本查看 复制代码

Program test_wrnf
  implicit none
  real :: epa(9,5)
  call EllipseP(epa)
contains
  subroutine EllipseP(epa)
    implicit none
    real,intent(out) :: epa(:,:)
    integer :: i, j
    epa=0.
    open (103 ,file='cep.txt') 
    call notes103()
    do i = 1 , size(epa,dim=2)
      read(103,*) epa(:,i)
    enddo
    close(103)
  end subroutine EllipseP
  subroutine notes103
    character(len=80) :: ax
    Do
      read(103,'(a80)') ax
      if ( index( ax , ']' ) > 0 ) return
    End Do
  end subroutine notes103
End Program test_wrnf