24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 486  |  回复: 2

zouji1983

金虫 (正式写手)

[求助] 请帮忙看一个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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gyf1017

金虫 (小有名气)

这个是做什么的程序?作为一个初学习fortran语言的入门者这个还是很难得
2楼2012-11-24 10:04:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zouji1983

金虫 (正式写手)

引用回帖:
2楼: Originally posted by gyf1017 at 2012-11-24 10:04:21
这个是做什么的程序?作为一个初学习fortran语言的入门者这个还是很难得

就是用Newton-Raphon解非线性方程的,我照着网上的例子编的,但是这个错误我不晓得怎么回事啊
3楼2012-11-26 18:06:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zouji1983 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 285求调剂 +3 AZMK 2026-04-02 4/200 2026-04-02 23:58 by wwytracy
[考研] 085602 找调剂 +3 逆时针快乐 2026-04-02 3/150 2026-04-02 21:23 by dongzh2009
[考研] 土木304求调剂 +4 兔突突突, 2026-04-02 5/250 2026-04-02 21:16 by 兔突突突,
[基金申请] esi高被引论文是不是能对中标有所加分和帮助呢 +4 redcom 2026-04-01 5/250 2026-04-02 20:48 by redcom
[考研] 319求调剂 +16 太容易1018 2026-04-01 16/800 2026-04-02 20:12 by seattle40
[考研] 一志愿华东理工大学,080500学硕,317分,求调剂 +13 s1145 2026-03-31 14/700 2026-04-02 18:55 by 哦哦嗯哈
[考研] 一志愿南昌大学324求调剂 +12 hanamiko 2026-04-01 12/600 2026-04-02 14:51 by 5896
[考研] 275求调剂 +13 jjjjjjjjjjl 2026-03-27 13/650 2026-04-02 13:07 by yulian1987
[考研] 求调剂,一志愿 南京航空航天大学 ,080500材料科学与工程学硕,总分289分 +11 @taotao 2026-03-29 11/550 2026-04-02 10:04 by realme321
[考研] 一志愿北交大材料工程,总分358 +4 cs0106 2026-04-01 4/200 2026-04-02 07:42 by 尚水阁主
[考研] 材料求调剂 +8 呢呢妮妮 2026-04-01 8/400 2026-04-02 07:13 by yjolah
[考研] 求调剂0703 +5 周嘉尧 2026-03-31 8/400 2026-04-01 20:32 by ltltkkk
[考研] 求调剂 +4 DADA怪 2026-03-31 4/200 2026-04-01 14:30 by ZXlzxl0425
[考研] 267求调剂 +13 uiybh 2026-03-31 13/650 2026-04-01 10:25 by 探123
[考研] 301求调剂 +8 axibli 2026-04-01 8/400 2026-04-01 09:51 by 我的船我的海
[考研] 335求调剂 +3 321* 2026-03-31 4/200 2026-04-01 00:00 by 321*
[考研] 抱歉 +3 田洪有 2026-03-30 3/150 2026-03-30 19:11 by 迷糊CCPs
[考研] 药学105500求调剂 +3 Ssun。。 2026-03-28 3/150 2026-03-28 11:24 by lxf170613
[考研] 085405 考的11408求各位老师带走 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 09:19 by 乐呵呵的追梦人
[考研] 352分 化工与材料 +5 海纳百川Ly 2026-03-27 5/250 2026-03-28 03:39 by fmesaito
信息提示
请填处理意见