作者py94g (...)
看板Fortran
標題[請益] FORTRAN讀檔案的問題
時間Mon May 7 15:19:04 2012
如果把檔案給分開,應該就不會有上一個迴圈如何寫的問題了。
然而分成多個檔案,每一個檔案是3行10000列。
檔名是test_000001 到 test_000010 有10個檔,只有數值沒Header沒有附檔名。
我寫的程式可以Compiler但是就是有錯...不知錯在哪裡,
想和迴圈問題一起請教大家,
以下是我寫的幾行:
Program read
Dimension POINT(3, 10000) ! 設定矩陣維度,3行N列。
! 讀test_000001 ...到 test_000010 有10個檔
character(len=11) inputname ! 檔案有11個字元
do i = 1,10 ! 我要讀入10個檔案
write(inputname,"(a5,i6)")'test_',000000+i
C 請問大家,我這一行有寫錯嗎?
open(1,file=inputname)
read(1,*) POINT ! 讀進
write(*,10) POINT ! 秀在螢幕
enddo
10 FORMAT(1X, 3(3F10.5)) ! 輸出的檔案格式 Raw data。
Stop
END
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.20.237
1F:→ callmei:我還沒看內文的CODE 只是想先跟你說 05/07 15:40
2F:→ callmei:大多時候用DEBUG模式編譯的程式下去跑 可以知道問題在哪裡 05/07 15:41
3F:→ py94g:我覺得問題在檔名的000001~000010,我再試試看改檔名 05/07 16:43
4F:→ py94g:test_000001改成test_1...test_000010改成test_10 05/07 16:49
5F:→ py94g:原 write(inputname,"(a5,i6)")'test_',000000+i 05/07 16:50
6F:→ py94g:改成write(inputname,"(a5,i1)")'test_',+i 05/07 16:50
7F:→ py94g:這樣可以跑,但是只能跑到test_9,而test_10出錯。 05/07 16:52
8F:→ py94g:原因是write(inputname,"(a5,i1)")'test_',+i中的i1 05/07 16:53
9F:→ py94g:檔案超過9就不行了,我想這是沒人會同意這種寫法的。 05/07 16:54