24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 488  |  回复: 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的回帖

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的回帖
查看全部 3 个回答

gyf1017

金虫 (小有名气)

这个是做什么的程序?作为一个初学习fortran语言的入门者这个还是很难得
2楼2012-11-24 10:04:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿深大085601材料工程专业(专硕)300分可以调剂去哪 +7 10160315 2026-04-02 7/350 2026-04-02 22:48 by liveagain126
[考研] 085600调剂 +3 1amJJ 2026-04-02 3/150 2026-04-02 22:00 by 科研小专家
[考研] 化工调剂303分,过四级 +14 栖梧待风 2026-04-02 14/700 2026-04-02 21:58 by 科研小专家
[考研] 312求调剂 +5 小小墨123 2026-04-02 6/300 2026-04-02 20:12 by dongzh2009
[考研] 301求调剂 +14 骆驼男人 2026-04-02 14/700 2026-04-02 14:08 by baoball
[考研] 一志愿北京科技大学材料学硕328分求调剂 +6 1段时间 2026-03-31 7/350 2026-04-02 13:57 by 3041
[考研] 材料考研调剂 +9 Gs大王 2026-04-02 9/450 2026-04-02 10:24 by olim
[考研] 324分 085600材料与化工 +20 呆鹅oor 2026-03-27 20/1000 2026-04-02 10:13 by oooqiao
[考研] 274求调剂 +8 薛定谔的虎。 2026-04-01 8/400 2026-04-02 09:21 by 晴空210210
[考研] 367求调剂 +8 芋泥啵啵… 2026-03-28 8/400 2026-04-02 09:20 by olim
[考研] 301求调剂 +13 A_JiXing 2026-04-01 13/650 2026-04-02 09:01 by sanrepian
[考研] 295材料工程专硕求调剂 +19 1428151015 2026-03-27 19/950 2026-04-01 22:34 by peike
[考研] 0856,材料与化工321分求调剂 +13 大馋小子 2026-03-27 14/700 2026-04-01 15:58 by yanflower7133
[考研] 求调剂 +4 DADA怪 2026-03-31 4/200 2026-04-01 14:30 by ZXlzxl0425
[考研] 086000生物与医药298调剂求助 +4 元元青青 2026-03-31 6/300 2026-04-01 11:13 by syh9288
[考研] 环境工程 085701,267求调剂 +15 minht 2026-03-29 16/800 2026-04-01 10:13 by li_sujuan99
[考研] 286求调剂 +6 Faune 2026-03-30 6/300 2026-03-31 14:37 by jp9609
[考研] 一志愿浙江大学工科动力工程370,数一121,专业课135,现在能去哪里 +3 080700调剂 2026-03-30 4/200 2026-03-31 12:00 by KLMY666
[考研] 086000生物与医药调剂 +5 Feisty。 2026-03-28 9/450 2026-03-29 12:02 by longlotian
[考研] 调剂 +3 李嘉图·S·路 2026-03-27 3/150 2026-03-27 11:19 by wangjy2002
信息提示
请填处理意见