文件数据读写

伊芸贾

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

program shujv1
  implicit none
  real,allocatable::poro(:)
  integer::k,stat
  character(len=160) :: line
  open(unit=4,file='.\poro.dat',iostat=stat)
  if(stat/=0) stop
  k=0
  Do
    read(4,'(a)',iostat=stat) line
    if (stat/=0)exit
    if ((line(1:2)/='--').and.(line/=' ').and.(line(1:6)/='poro')) k=k+1
  End Do
  allocate(poro(k))
  rewind(4)
  k = 0
  Do
    read(4,'(a)',iostat=stat) line
    if (stat/=0)exit
    if ((line(1:2)/='--').and.(line/=' ').and.(line(1:6)/='poro')) then
      k = k + 1
      read(line,*,iostat=stat) poro(k)
    end if
  End Do
  write(*,'(f10.4)') poro
end program shujv1