读写错误，请帮忙。

VBS_CS · 发表于 2016-5-19 17:31:22

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

01  PROGRAM TEST_WRNF
02  IMPLICIT NONE
03  Real, Dimension(9,62)  ::EPa
04  CALL EllipseP(EPa)
05  END
06 
07       Subroutine EllipseP(EPa)
08       IMPLICIT NONE
09       Real,    Intent(OUT), Dimension(9,62) :: EPa
10       Integer :: I, J
11       EPa=0.
12       OPEN (103 ,FILE='CEP.txt', STATUS='UNKNOWN')
13       Call Notes103
14       Do I=1,62
15           Read(103,*) (EPa(J,I), J=1,9)
16       Enddo
17       CLOSE(103)
18       End  Subroutine  EllipseP
19 
20      Subroutine Notes103
21      CHARACTER*1 AX(80)
2211    Read(103,'(80A1)')(AX(I),I=1,80)
23      Do I=1,80
24        If(AX(I)==']') Return
25      Enddo
26      GOTO 11
27      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:48

你把很多问题想复杂了。

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

01Program test_wrnf
02  implicit none
03  real :: epa(9,5)
04  call EllipseP(epa)
05contains
06  subroutine EllipseP(epa)
07    implicit none
08    real,intent(out) :: epa(:,:)
09    integer :: i, j
10    epa=0.
11    open (103 ,file='cep.txt')
12    call notes103()
13    do i = 1 , size(epa,dim=2)
14      read(103,*) epa(:,i)
15    enddo
16    close(103)
17  end subroutine EllipseP
18  subroutine notes103
19    character(len=80) :: ax
20    Do
21      read(103,'(a80)') ax
22      if ( index( ax , ']' ) > 0 ) return
23    End Do
24  end subroutine notes103
25End Program test_wrnf

左耳北上南 · 发表于 2016-6-13 17:28:44

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

[文件读写] 读写错误，请帮忙。

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

规矩勋章

元老勋章

水王勋章