利用FORTRAN 读取文件行数的小代码

topmad

每次读文件，都要先看看有多少行，烦得很啊，有木有啊，要是直接打开多好啊，可是自己太笨，想不出好办法

直到有一天，在自己面前堆了几千个文件，每个文件的行数都不同，瞬间石化了{:soso_e136:}{:soso_e136:}

必须找一个好的方法了，感谢baidu，google，还有万能的清风{:soso_e102:}，热心的深深版主{:soso_e151:}.....

终于，找到了{:soso__1886989460682607670_4:}

不看不知道  看了就悔恨了，原来这么短的小代码，就能实现

还在用老方法看行数？ 你out了。。。。{:soso_e122:}{:soso_e122:}{:soso_e122:}{:soso_e122:}{:soso_e122:}  （虽然3分钟之前，我还在out行列{:soso_e133:}{:soso_e133:}{:soso_e133:}{:soso_e133:}

不多说啦，直接上

implicit none
real line
integer i
i=0  !清风提醒，特加此行
open(10,file='topmad.txt',status='old')
do while (.true.)
read(10,*,end=100) line
i=i+1
enddo
100 continue
print*,i

end

然后 i 就是行数啦  
不信你试试

何永利

idl中直接有现成的函数，file_lines，不过估计算法是一样的

mofangbao

我还真希望我是万能的，哈哈
额。。。还是发现了一个小漏洞，I的初始值应该为
integer i=0

topmad

何永利 发表于 2011-7-14 23:45

                               
登录/注册后可看大图

idl中直接有现成的函数，file_lines，不过估计算法是一样的

何版主也来捧场啦，感谢感谢  没怎么接触过idl  可能后面的工作用得到 还望多多指教

何永利

topmad 发表于 2011-7-15 09:10

                               
登录/注册后可看大图

何版主也来捧场啦，感谢感谢  没怎么接触过idl  可能后面的工作用得到 还望多多指教

:lol，我也要用fortran的，当然要向你们学习啊！

新小子

这个蛮好，有没有处理列数不一样的哦，呵呵……

topmad

新小子 发表于 2011-7-15 10:21

                               
登录/注册后可看大图

这个蛮好，有没有处理列数不一样的哦，呵呵……

估计我也要碰到同样的问题了  如果研究出来了  就发上来

言深深

mofangbao 发表于 2011-7-15 07:52

                               
登录/注册后可看大图

我还真希望我是万能的，哈哈
额。。。还是发现了一个小漏洞，I的初始值应该为
integer i=0

哈哈，其实fortran里面变量不定义的话，默认的全是0。然后我比较喜欢用I-N规则，连implicit none都可以省略啦！
不知，TOP的这个程序挺不错，很赞啊

artmunich

我以前用这样计算的
read(10,"(A160)")tempstr
j=0
statnum=0
do while(.not.eof(10))
read(10,"(I5)")i
if(i/=j)then
j=i
statnum=statnum+1
end if
end do

topmad

artmunich 发表于 2011-7-15 22:20

                               
登录/注册后可看大图

我以前用这样计算的
read(10,"(A160)")tempstr
j=0

这个方法我也看到了，我感觉那个更加直观，就放上去了，呵呵，a版好勤奋啊