| 查看: 447 | 回复: 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 |
回复此楼» 猜你喜欢
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
孩子确诊有中度注意力缺陷
已经有12人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
请教 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