作者fermmy (羃冪)
看板Fortran
標題[問題] Fortran如何抓取指定間隔列的資料
時間Thu May 18 20:12:14 2017
大家好
想問大家一個很基本的可是已經困擾我很久的問題
想要讓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
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 42.73.157.183
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Fortran/M.1495109537.A.E86.html
1F：→ rex0707: a1向量大小為2 可是i=1~6 那a1(i)不就會超出你設的大小05/19 01:29
2F：→ rex0707: 而且看不太懂你的問題 是想要每兩行輸出一次資料看有沒有05/19 01:30
3F：→ rex0707: 錯嗎?05/19 01:30
謝謝r大的回復~我有再修改一下文章，我文中的M寫錯了應該是6。對，是希望透過這個簡
單的範例每2行輸出資料，看程式有沒有寫錯。
※ 編輯: fermmy (223.137.251.225), 05/19/2017 08:57:50
4F：→ jubilee2: 2解 1:if(mod(i,2)==1) 2:拆成兩個do迴圈(一個讀一個寫)05/19 09:34
5F：→ Cypresslin: 如果你資料總列數為偶數x，那就是用一個i=1,x/2的迴圈05/19 11:55
6F：→ Cypresslin: 裡面放兩個read，一個讀a1, a2一個不放變數05/19 11:56
7F：→ Cypresslin: 題外話，可以用shell的話 sed -n '1~2p' 指令一行搞定05/19 11:59
8F：→ rex0707: 應該是每行資料都要讀 只是每兩行輸出一次檢查05/19 12:20
9F：→ rex0707: do i = 1 , M / read a1 a2 / if (mod==1) write / enddo05/19 12:21
10F：→ rex0707: 改一下 才能輸出偶數行資料 (mod(i,2)==0)05/19 13:27
11F：→ blc: do; read; read; write; end do05/19 17:52
謝謝大家的解惑～!!!已經用mod指令寫出來了，再次謝謝各位m(_ _)m
※ 編輯: fermmy (114.36.71.33), 05/19/2017 23:36:43
12F：→ janhon: real ? integer 吧?! 05/21 04:37