[问题] Fortran如何抓取指定间隔列的资料

楼主: fermmy (羃幂)   2017-05-18 20:12:14
大家好
想问大家一个很基本的可是已经困扰我很久的问题
想要让Fortran读一个记事本,然后几行输出结果,类似excel里面的offset指令的功能,
其实也是可以用excel做可是因为记事本很大太费时,想请问应该如何用Fortran写。
范例档案.txt
15 76
85 69
52 80
25 31
65 41
我写的
program text
REAL a1, a2,l
parameter (M=6)
DIMENSION a1(M), a2(M)
open(11,file='file.txt')
k=6
do i=1,k
read(11,*) a1(i), a2(i)
l=i*2 !两行想取一次
if(l<k) then
write(*,*) l !这里有跑出来l的值可是下面write部分就全部都是0了
write(*,*) a1(l), a2(l)
endif
enddo
end
想要输出档案为
15 76
52 80
65 41
因为好几年没有碰Fortran了...知道问题可能有点笨...不好意思,麻烦大家了,谢谢大
家m(_ _)m
作者: rex0707 (我就是我)   2017-05-19 01:29:00
a1向量大小为2 可是i=1~6 那a1(i)不就会超出你设的大小而且看不太懂你的问题 是想要每两行输出一次资料看有没有错吗?
作者: jubilee2 (3321)   2017-05-19 09:34:00
2解 1:if(mod(i,2)==1) 2:拆成两个do循环(一个读一个写)
作者: Cypresslin (啊哈哈~~)   2017-05-19 11:55:00
如果你资料总列数为偶数x,那就是用一个i=1,x/2的循环里面放两个read,一个读a1, a2一个不放变量题外话,可以用shell的话 sed -n '1~2p' 指令一行搞定
作者: rex0707 (我就是我)   2017-05-19 12:20:00
应该是每行资料都要读 只是每两行输出一次检查do i = 1 , M / read a1 a2 / if (mod==1) write / enddo改一下 才能输出偶数行资料 (mod(i,2)==0)
作者: blc (Anemos)   2017-05-19 17:52:00
do; read; read; write; end do
作者: janhon (海牛)   2017-05-21 04:37:00
real ? integer 吧?!

Links booklink

Contact Us: admin [ a t ] ucptt.com