fortran读文本文件—无法估计数据数目问题

        读文件的时候有时候会遇到无法事先估计数据数目的情况。这个问题最传统的解决方法是声明一个“超级巨大”的数组来保存数据，而这个数组在使用时大部分空间都是闲置的。虽然fortran可以支持可变大小的数组（allocatable），但是这种数组在使用之前也必须确定大小才能分配内存。

        C/C++里面可以通过指针来创建动态数组，现在fortran也支持指针，是不是也可以由类似的应用呢？答案是肯定的。fortran里面也可以通过指针来创建“串行结构”，可以配合需要向内存要求刚好的空间，只有一条数据时就要求一条数据的空间，有一千条数据时就要求一千条数据的空间，这样就可以有效地使用内存。

        下面通过一个例子来说明。要读取的数据放在一个叫data.txt的文本文件中，现在知道文件中有四列数据，但不知道文件中有几行数据。

        这是数据：



        下面是读取代码：

module linklist
implicit none
type::datalink
real*8::item(1:4) !声明一个1*4的数组用来存放一行数据
type(datalink),pointer::next !声明指针，用于指向下一条数据
end type
end module linklist

program readlinebyline
use linklist
implicit none
character*20 :: filename
integer::datasize,error,i

type(datalink),pointer::head,p

filename="data.txt"

open(10,file=filename,status="old",iostat=error)
if(error/=0) then
write(*,*) "Open file fail!"
stop
end if

allocate(head) !给head分配内存
nullify(head%next) !使head%next指向为空
p=>head

do while(.true.)
read(10,fmt=*,iostat=error) p%item  !读取一行数据赋给p%item
if(error/=0) exit
datasize=datasize+1 !datasize用于统计行数
allocate(p%next,stat=error) !添加下一个数据
if(error/=0) then
write(*,*) "Out of memory!"
stop
end if
p=>p%next !移动到下一行数据
nullify(p%next)

end do
close(10)

write(*,'("The total number of lines:",I3)') datasize

p=>head
do while(.true.)
write(*,*) p%item
if(.not. associated(p%next)) exit !如果p%next没有设置指向，就退出
p=>p%next
end do
end program readlinebyline

        程序中的注释应该写得很清楚了，按行读取数据，最后输出在屏幕上，结果如下：



        指针果然很好很强大啊，这只是一个很简单的小例子。但是串行方法也有一个缺点，就是只适合顺序读取数据，不适合作随机读取。比较好的编程方法是先用串行来读取文件，在获知文件大小之后将数据复制到一个可变大小的数组中，再把串行全部删除，这样就可以通过数组方便地进行查询了。