作者moneyyoyo (钱阿宥)
看板Fortran
标题[问题] 想请问读大量档案的问题
时间Sun Aug 16 01:42:06 2009
因为模拟的关系，有很多笔资料要整理
其实也只是整理出资料的平均跟均方根
之前因为模拟的时间并不长（产生１０００笔资料）
所以没有发现问题
可是现在加长模拟时间
所以资料数也变多了（约８０００笔）
但是却在读取第７０００笔资料完後出现错误
forrtl:sever (24): end-of-file during read, unit 7000,
file E:\48\data\windprofile6-8_20b(uctr)-7100
Image PC Routine Line Source
ctr_profile.exe 00412669 Unknown Unknown Unknown
ctr_profile.exe 004124C7 Unknown Unknown Unknown
ctr_profile.exe 004116A4 Unknown Unknown Unknown
ctr_profile.exe 00411AD9 Unknown Unknown Unknown
ctr_profile.exe 0040CA15 Unknown Unknown Unknown
ctr_profile.exe 00401273 Unknown Unknown Unknown
ctr_profile.exe 0044D3A9 Unknown Unknown Unknown
ctr_profile.exe 00432839 Unknown Unknown Unknown
kernel132.dll 7C816FE7 Unknown Unknown Unknown
Incrementally linked image--PC correlation disabled.
其实我心里是觉得不是程式问题而是矩阵的问题
爬文之後又觉得难道是我矩阵太大无法负担吗？
所以想请问板上高手是否能指点迷津
感谢~
－－以下为程式－－－－－
program aa
implicit none
character:: onput_file*10 ,input_file*10
CHARACTER FNM*40,FNM1*40,FNM2*40,word(20)*40,wordaa,temp1(20)*40
integer i,ii,jj,a,j,a1,a2,kk,jjj,n,data_n
real point,x(50),y(50),z(50)
real filedata(19,10,50),temp(19,10,50)
real velocity(19,8000,50),mean_V(19,50),rms_V(19,50)
real velocityx(19,8000,50),mean_Vx(19,50),rms_Vx(19,50)
real velocityy(19,8000,50),mean_Vy(19,50)
real velocityz(19,8000,50),mean_Vz(19,50)
write(*,*)"资料点数目:"
data_n=19 !有19条线要读取
OPEN(7000,FILE="FILENAME(plctr).txt")
OPEN(8000,FILE="mean_V(plctr).dat")
OPEN(8500,FILE="V(plctr).xls")
!OPEN(9000,FILE="mean_dV(plctr).dat")
!----------------------------------------
!***************** DATA *****************
!----------------------------------------
READ(7000,*) KK
!--------每笔资料读取--------------------
DO J=1,KK
READ(7000,'(a)') FNM
fnm1=fnm(1:len_trim(fnm))
OPEN(J,FILE=FNM1)
!plane-ctr point 9520
DO n=1,data_n
read(J,*)word(n),wordaa,point
read(J,*)
do i=1,point
read(J,*)filedata(n,1,i)
end do
read(J,*)
read(J,*)
do i=1,point
read(J,*)filedata(n,2,i)
end do
read(J,*)
read(J,*)
do i=1,point
read(J,*)filedata(n,3,i)
end do
read(J,*)
read(J,*)
!风速
do i=1,point
read(J,*)filedata(n,4,i)
end do
read(J,*)
read(J,*)
do i=1,point
read(J,*)filedata(n,5,i)
end do
read(J,*)
read(J,*)
do i=1,point
read(J,*)filedata(n,6,i)
end do
read(J,*)
read(J,*)
do i=1,point
read(J,*)filedata(n,7,i)
end do
read(J,*)
read(J,*)
.........(略)
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 61.229.169.186
1F：推 yangzi:出现什麽错误讯息 08/16 11:06
※ 编辑: moneyyoyo 来自: 61.229.183.215 (08/16 12:12)
2F：→ UDP:语法没错，执行错误，是酱子吗? 08/16 19:40
3F：推 yangzi:检查一下第7000笔资料有无异常或错误 08/16 22:39
4F：→ yangzi:找不出来也可以尝试改顺序 把第7000笔先拿走 看其他能否ok 08/16 22:40
5F：→ yangzi:如果不是资料异常的话 就要检查其他问题了XD 08/16 22:40
哈哈~不好意思我找到错误了....
我错在很笨的地方
OPEN(7000,FILE="FILENAME(plctr).txt")
改成
OPEN(70000,FILE="FILENAME(plctr).txt")
至於为什麽...我不会解释>///////////<
谢谢大家罗
※ 编辑: moneyyoyo 来自: 163.13.138.187 (08/17 13:08)