[请益] FORTRAN读档案的问题

楼主: py94g (...)   2012-05-07 15:19:04
如果把档案给分开,应该就不会有上一个循环如何写的问题了。
然而分成多个档案,每一个档案是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
作者: callmei (挨)   2012-05-07 15:40:00
我还没看内文的CODE 只是想先跟你说大多时候用DEBUG模式编译的程式下去跑 可以知道问题在哪里
楼主: py94g (...)   2012-05-07 16:43:00
我觉得问题在档名的000001~000010,我再试试看改档名test_000001改成test_1...test_000010改成test_10原 write(inputname,"(a5,i6)")'test_',000000+i改成write(inputname,"(a5,i1)")'test_',+i这样可以跑,但是只能跑到test_9,而test_10出错。原因是write(inputname,"(a5,i1)")'test_',+i中的i1档案超过9就不行了,我想这是没人会同意这种写法的。
作者: gilocustom   2012-05-07 17:19:00
#18pYwxH_

Links booklink

Contact Us: admin [ a t ] ucptt.com