作者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