唉，回家慢慢想吧。

一行一行地读就可以了，不过要加上iostat=nstat子句，否则最后一行会跳出的：）

既然知道数据总数N，每行的数据数M个，就知道完整N/M行，最后一行剩余R=mod(N,M)个。然后不就可以控制读出么？

昨晚回去想了一下，今天又弄了一个早晨，后处理程序差不多完成了。
昨天造成数据一直没有读成功的原因是程序的语法错误，我好多时候都犯了语法错误，郁闷！
我定义了一个数组读数据，将数据保存在数组中。read(1,*) (a(i),i=1,np),等到需要的时候，写完数据后，再释放内存。这样做的缺点是，如果总节点数np比较大的话，占用内存就大了。我一时也想不到什么好办法。
我原本也想采用mod(N,M)函数控制的，但是又觉得代码多了一些。
还有一个问题，请指导一下。我定义了一个字符串数组，读入每行的字符串，但是如果程序这样写的话，无法读入数据，一直报错：
    read(1,'(a)',iostat=error)   strN(i)
这样写的话，就可以读入数据。
    read(1,'(a)',iostat=error)   astr
        strN(i)=astr
为什么产生这样的问题？[em10][em10][em18]

read(1,'(a)',iostat=error)   strN(i:i)

是这样读啊,晕死了[em10]

呵呵，是的，不过建议还是一贴一个问题比较好些：）
且解决后尽快结贴免得给广告机会：）