文件读取问题

于鑫常青藤 · 发表于 2016-4-12 14:25:11

尝试先把文件内容读入字符串，然后再从字符串里读取信息，但是没有成功。请高手指点！
1  0.1000E+01
2  0.1414E+01
3  0.1732E+01
4  0.2000E+01
5  0.2236E+01
6  0.2449E+01
7  0.2646E+01
8  0.2828E+01
9  0.3000E+01
10  0.3162E+01
这是待读取的文件数据

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

01program ex33
02    implicit none
03    integer :: a, b, fileid, k, error
04    character (len=512):: str
05    character (len=512) :: filename="ex32.txt"
06    open (newunit=fileid, file=filename, form="formatted", access="direct", recl=512, status="old")
07    do while (.true.)
08        read (*,*) k
09        read (fileid, "(a512)", rec=k, iostat=error) str
10        if (error/=0) exit
11        read (str, *) a, b
12        write (*,*) a, b
13    end do
14    close (fileid)
15end

这是读取的代码（尝试先读入字符串，再读取）但是没有成功，不知道错在哪里，请高手指点！

vvt · 发表于 2016-4-12 15:13:51

你的文件每行的宽度不一样，无法使用直接方式读取。请用顺序方式读取

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

01program ex33
02  implicit none
03  integer :: i , a, b, fileid, k, error
04  character (len=512):: str
05  character (len=512) :: filename="ex32.txt"
06  open (newunit=fileid, file=filename , status="old") !//采用有格式顺序读写
07  do
08    read (*,*) k
09    rewind(fileid)
10    do i = 1 , k - 1
11      read( fileid , * )
12    end do
13    read ( fileid, "(a512)", iostat=error ) str
14    if (error/=0) exit
15    read (str, *) a, b
16    write (*,*) a, b
17  end do
18  close (fileid)
19end

于鑫常青藤 · 发表于 2016-4-12 15:31:33

open (newunit=fileid, file=filename, status="new", access="sequential", form="formatted", recl=16)
do i=1,10,1
!write (fileid,"(i2, 2xE10.4)") i, sqrt(real(i))
end do
close(fileid)
带读取的文件是按照统一的格式输出到“ex32.txt”中的，宽度应该一样吧？

于鑫常青藤 · 发表于 2016-4-12 15:33:36

vvt 发表于 2016-4-12 15:13
你的文件每行的宽度不一样，无法使用直接方式读取。请用顺序方式读取

[mw_shl_code=fortran,true]program ...

open (newunit=fileid, file=filename, status="new", access="sequential", form="formatted", recl=16)
do i=1,10,1
!write (fileid,"(i2, 2xE10.4)") i, sqrt(real(i))
end do
close(fileid)
带读取的文件是按照统一的格式输出到“ex32.txt”中的，宽度应该一样吧？