文件读取问题

于鑫常青藤 · 发表于 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] 纯文本查看 复制代码

program ex33
    implicit none
    integer :: a, b, fileid, k, error
    character (len=512):: str
    character (len=512) :: filename="ex32.txt"
    open (newunit=fileid, file=filename, form="formatted", access="direct", recl=512, status="old")
    do while (.true.)
        read (*,*) k
        read (fileid, "(a512)", rec=k, iostat=error) str
        if (error/=0) exit
        read (str, *) a, b
        write (*,*) a, b
    end do
    close (fileid)
end

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

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

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

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

program ex33
  implicit none
  integer :: i , a, b, fileid, k, error
  character (len=512):: str
  character (len=512) :: filename="ex32.txt"
  open (newunit=fileid, file=filename , status="old") !//采用有格式顺序读写
  do
    read (*,*) k
    rewind(fileid)
    do i = 1 , k - 1
      read( fileid , * )
    end do
    read ( fileid, "(a512)", iostat=error ) str
    if (error/=0) exit
    read (str, *) a, b
    write (*,*) a, b
  end do
  close (fileid)
end

于鑫常青藤 · 发表于 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”中的，宽度应该一样吧？