作者HAKUKE (亚流师)
看板Fortran
标题[问题] Error
时间Tue Jan 26 23:29:10 2010
forrtl severe 24 end-of-file during read ,unit 11
forrtl severe 161 Program Exception - array bounds exceeded
请问有人知道这两个错误是什麽意思呢
网路上的我还是看不太懂该怎麽改正
code是先从许多数据中读取
然後算内插
第二个error发生时是有15038笔数据
他的意思是超出电脑记忆体的容量吗 还是code本身的问题勒
谢谢!
---------------
implicit none
integer*4::np,ng,b(2)
integer*4, allocatable::cvg_neib(:,:,:),ng_neib(:,:)
integer*4::i,j,nn
real*8,
allocatable::xp(:),yp(:),up(:),vp(:),pp(:),xg(:),yg(:),ug(:),vg(:),pg(:)
real*8::dx_ini,dy_ini,sh,mass,dt
real*8::x_max,x_min,y_max,y_min
real*8::w,wx,wy,wsum
real*8::a(10)
character*100::filename1,filename2
common/LL/x_max,x_min,y_max,y_min
common/all2/dx_ini,dy_ini,sh,mass,dt
! np = no. of particles
np=1600
! ng = no. of grids
ng=40
! entire domain including ghost paricles
x_max=2.000000E-03
x_min=-1.000000E-03
y_max=2.000000E-03
y_min=-1.000000E-03
! initial particle spacing
dx_ini=0.000025
dy_ini=0.000025
mass=0.62500006E-00
! smoothing length
sh=0.25000001E-04
dt=5.e-5
allocate(cvg_neib(2,ng,200),ng_neib(2,ng))
allocate(xp(np),yp(np),up(np),vp(np),pp(np),xg(ng),yg(ng),ug(ng),vg(ng)
,pg(ng))
write(*,*)'Keyin the filename of the input file'
read(*,*)filename1
write(*,*)
write(*,*)'Keyin the filename of the output file'
read(*,*)filename2
write(*,*)
open(11,file=filename1)
open(21,file=filename2)
! for data of particles
do i=1,11
read(11,*)
end do
do i=1,np
read(11,*)(a(j),j=1,10),b(1),b(2)
xp(i)=a(2)
yp(i)=a(3)
up(i)=a(4)
vp(i)=a(5)
pp(i)=a(8)
end do
! for data of grids
do i=1,ng-1
xg(i)=0.5
yg(i)=dy_ini*i
end do
call grid_interpolation(xp,yp,up,vp,pp,np,xg,yg,ug,vg,pg,ng)
write(21,*)'variables = "y"'
write(21,*)' "u" '
do i=1,ng
write(21,'(3(es12.5,1x))')yg(i),ug(i),vg(i)
end do
pause 'ok'
end
--
我要读的档案 有1600笔 每笔里面有x y u v等共11个物理量
标颜色那 是我看不懂的地方
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.250.123
※ 编辑: HAKUKE 来自: 140.112.250.123 (01/26 23:29)
1F：推 laa7352:方便的话程式码贴上来看看 01/27 00:07
2F：推 egggg:第一个，你读的档案内容比你程式想读的要少 01/27 00:13
3F：→ egggg:第二个 你叫用的Array位置，比你宣告的要大 01/27 00:16
4F：→ egggg:例如 integer:: A(100) 你却在程式中叫了A(101)这种情况 01/27 00:16
※ 编辑: HAKUKE 来自: 140.112.250.123 (01/27 01:01)
5F：→ HAKUKE:先谢谢楼上 01/27 01:06
※ 编辑: HAKUKE 来自: 140.112.250.123 (01/27 01:07)