| 查看: 522 | 回复: 2 | ||
[求助]
请帮忙看一个fortran小程序
|
|
请高手帮忙看看这个程序有什么问题，运行的时候显示：array bounds exceeds. 非常感谢！ subroutine convert_composition(ix,iy) integer :: ix, iy, k real(8) :: tol,f(xp(1,1,ix,iy)),df(xp(1,1,ix,iy)),xp0(1,1,ix,iy) f(xp(1,1,ix,iy))=1.06366481d10+7.94800172d9+1.72992961d9*log(xp(1,1,ix,iy))-1.72992961d9*log(1-xp(1,1,ix,iy))-1.182988286d10*(1-2*xp(1,1,ix,iy))+2.75714286d9*(6*xp(1,1,ix,iy)-6*xp(1,1,ix,iy)*xp(1,1,ix,iy)-1)+3.083731429d9*(24*xp(1,1,ix,iy)*xp(1,1,ix,iy)-10*xp(1,1,ix,iy)-16*xp(1,1,ix,iy)*xp(1,1,ix,iy)*xp(1,1,ix,iy)+1))-dmudx(1,1,2,ix,iy)*(x(1,ix,iy,old)-phi(1,ix,iy)*xp(1,1,ix,iy))/phi(2,ix,iy) + 9.229804981d10 df(xp(1,1,ix,iy))=2*1.182988286d10+1.72992961d9/xp(1,1,ix,iy)/(1-xp(1,1,ix,iy))+2.75714286d9*(6-12*xp(1,1,ix,iy))+3.083731429d9*(48*xp(1,1,ix,iy)-10-48*xp(1,1,ix,iy)*xp(1,1,ix,iy))+dmudx(1,1,2,ix,iy)*phi(1,ix,iy)/phi(2,ix,iy) tol=1e-6 xp0(1,1,ix,iy)=0.2 k=0 10 k=k+1 xp(1,1,ix,iy)=xp0(1,1,ix,iy)-f(xp0(1,1,ix,iy))/df(xp0(1,1,ix,iy)) if(abs(xp(1,1,ix,iy)-xp0(1,1,ix,iy)).lt.tol)goto 20 xp0(1,1,ix,iy)=xp(1,1,ix,iy) goto 10 20 write(*,*) xp(1,1,ix,iy) xp(1,2,ix,iy)=(x(1,ix,iy,old)-phi(1,ix,iy)*xp(1,1,ix,iy))/phi(2,ix,iy) write(*,*) xp(1,2,ix,iy) end subroutine convert_composition |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
请教 fortran 运行错误的原因
已经有13人回复
- 请热心人能提供一个带imsl库的intel visual fortran编译器下载链接
已经有6人回复
- 大家帮我看一个fortran的程序，我总是计算不出正确的结果
已经有20人回复
- 我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？
已经有4人回复
- 【已完结】求助：哪位帮忙写个简单的fortran小程序？谢谢。
已经有6人回复
- 求助高手帮忙编一个Fortran的小程序，重金悬赏......
已经有15人回复
- 求助计算机高人，编写个Fortran程序
已经有13人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- 请教关于in cell形式向original形式坐标转化的fortran算法
已经有3人回复
- 有关fortran的一次而问题，希望大家能帮帮忙，谢谢
已经有4人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】FORTRAN写文件时每行开头都多了一个空格？
已经有6人回复
- 【求助】用fortran怎么描述下面这样的情况，我是新手，请大家帮忙
已经有8人回复
- 【求助】请问一下fortran的可视化编程
已经有7人回复
- 【求助】这样的公式，在fortran里该怎么写了，谢谢【求助成功】
已经有4人回复
- 【求助】请问有拟合分子力场参数的FORTRAN程序吗？
已经有13人回复
- 【交流】Fortran语言答疑专帖
已经有157人回复
2楼2012-11-24 10:04:21
3楼2012-11-26 18:06:45
