fortran程序数组越界?急求助……
老师给的一段fortran……本人根本没学过,硬着头皮改,由于程序较长,我就截一段出来,跪求大神帮忙!1--H2,2--CO,3--CO2,4--CH3OH,5--H2O...
老师给的一段fortran……本人根本没学过,硬着头皮改,由于程序较长,我就截一段出来,跪求大神帮忙
! 1--H2,2--CO,3--CO2,4--CH3OH,5--H2O,6--DME
real ccat(100),aa(8),bb(8),af1(8),af2(8),y0(8),cg(100,8)&
,cl(100,8),s(8),ps(100,8),par(100,8),ys(100,8),rbz(100,8),&
fm(100,8),akla(100),zgfa(100,8),zgfb(100,8),zgfc(100,8),&
zgfd(100,8),zgfs(8),gg(100,8),q(8),zgfc1(100,8),zgfd1(100,8)&
,yout(8),cph(100,8),fz(100,8),YI(8)
real k1,vin,rm1,rd1
real k2,k3,kco,kco2
real kp1,kp2,kp3
real kwh,km
character abc*8,def*8,xyz*8
省略好多……
! 计算催化剂床层分布
i=0
do 200 k=1,N
i=i+1
z=i*al/N
ccat(0)=(aMcat/vj)*a/(exp(a)-1.0)
ccat(i)=ccat(0)*exp(a*z/al)
200 continue
do 202 i=0,N,12
write(2,*)'ccat(',i,')',ccat(i)
202 continue
出错的就是计算催化剂床层的那一块,后面也有类似的错误,错误提示是
错误 19 error #5561: Subscript #1 of the array CCAT has value 0 which is less than the lower bound of 1
百度后有人说是数组越界,但是完全不知道如何改……跪求,万分感谢 展开
! 1--H2,2--CO,3--CO2,4--CH3OH,5--H2O,6--DME
real ccat(100),aa(8),bb(8),af1(8),af2(8),y0(8),cg(100,8)&
,cl(100,8),s(8),ps(100,8),par(100,8),ys(100,8),rbz(100,8),&
fm(100,8),akla(100),zgfa(100,8),zgfb(100,8),zgfc(100,8),&
zgfd(100,8),zgfs(8),gg(100,8),q(8),zgfc1(100,8),zgfd1(100,8)&
,yout(8),cph(100,8),fz(100,8),YI(8)
real k1,vin,rm1,rd1
real k2,k3,kco,kco2
real kp1,kp2,kp3
real kwh,km
character abc*8,def*8,xyz*8
省略好多……
! 计算催化剂床层分布
i=0
do 200 k=1,N
i=i+1
z=i*al/N
ccat(0)=(aMcat/vj)*a/(exp(a)-1.0)
ccat(i)=ccat(0)*exp(a*z/al)
200 continue
do 202 i=0,N,12
write(2,*)'ccat(',i,')',ccat(i)
202 continue
出错的就是计算催化剂床层的那一块,后面也有类似的错误,错误提示是
错误 19 error #5561: Subscript #1 of the array CCAT has value 0 which is less than the lower bound of 1
百度后有人说是数组越界,但是完全不知道如何改……跪求,万分感谢 展开
1个回答
展开全部
! 1--H2,2--CO,3--CO2,4--CH3OH,5--H2O,6--DME
real ccat(100)
改为
! 1--H2,2--CO,3--CO2,4--CH3OH,5--H2O,6--DME
real ccat(0:100)
real ccat(100)
改为
! 1--H2,2--CO,3--CO2,4--CH3OH,5--H2O,6--DME
real ccat(0:100)
更多追问追答
追问
call YII(tt,P,Y0,YI)同错,求改
cg(0,1)=YI(2)*y0(1)*pp/(r*tt)
cg(0,2)=YI(1)*y0(2)*pp/(r*tt)
cg(0,3)=YI(6)*y0(3)*pp/(r*tt)
Cg(0,4)=YI(3)*y0(4)*pp/(r*tt)
cg(0,5)=YI(5)*y0(5)*pp/(r*tt)
cg(0,6)=YI(4)*y0(6)*pp/(r*tt)
追答
同改
本回答被提问者采纳
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询