Fortran Coder

标题: 读写错误，请帮忙。 [打印本页]

作者: VBS_CS 时间: 2016-5-19 17:31
标题: 读写错误，请帮忙。

[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

作者: fcode 时间: 2016-5-19 21:12
你把很多问题想复杂了。

[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

作者: 左耳北上南 时间: 2016-6-13 17:28
第14行，你这里是有5行数据可以读入，不能写成62行。这里报错。
还有二楼的程序比你的简洁而且适用性广。

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)