24小时热门版块排行榜    

查看: 2867  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

sun805869434

金虫 (小有名气)

[求助] fortran里一个子程序不运行是怎么回事

整个fortran程序编译完,运行后就显示内存不能读取,切换到Debug,按F5后,可以发现浓度值C(J,K)计算结果不对,但是无论更改程序中该子程序的任何计算式,得到的结果都是一样的,好像跟我设置的计算式没有关系,这是没有读取相应的子程序吗?还是读取其他部分了?


SUBROUTINE COM
!--------------------------------------------------------
PARAMETER(NJ=120,NT=120)
COMMON U(NJ,NT),C(NJ,NT)
COMMON DER(NT),Y(NJ),DYL(NJ)
COMMON K,LITIM,DT
COMMON DLN,DS,RHOF,RHOS,DD,GR,M1,YDN,DXL,YL,DY,C0       
!--------------------------------------------------------
ENTRY VECTORY

U(J,0)=0
DO K=1,LITIM
DO J=1,M1
IF(J.NE.M1)THEN
U(J,K)=(GR*DT+U(J,K-1))/(DD*DT+1)
ELSE
U(J,K)=0
ENDIF
ENDDO
ENDDO

RETURN

!------------------------------------------------------------
ENTRY CONCENTRATE       

C(J,0)=C0
DO K=1,LITIM
DO J=1,M1
IF(J.EQ.1)THEN
C(J,K)=(U(J,K)*C(J,K-1))/U(J,K+1)
ELSEIF (J.GT.1.AND.J.LT.M1)THEN
C(J,K)=((U(J,K)*C(J,K-1))/U(J,K+1))+C(J-1,K-1)-C(J-1,K)
ELSE
C(J,K)=0
ENDIF
ENDDO
ENDDO

RETURN
这是子程序中的两个部分,第一个子程序运行后U(J,K)的结果是正确的,但是进行第二个子程序时,要用到上一个的结果U(J,K),运行结果C(J,K)是错误的,并且不管我怎么改C(J,K)的计算式,得出的结果不变,还是原来错误的结果,不知道原因出在哪。

[ Last edited by sun805869434 on 2013-6-6 at 10:14 ]
回复此楼

» 猜你喜欢

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

优秀是一种习惯
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

引用回帖:
7楼: Originally posted by sun805869434 at 2013-06-06 21:01:21
我的程序比较多,复制比较不方面,我把附件发上来,您能帮我看看吗?这是一个主程序,和一个子程序。谢谢。...

太多不行,你需要自己先简化。不能指望别人读1000条程序之类。我可以回答你的问题,帮助你提高技巧,但不是帮你干活
8楼2013-06-06 21:27:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

感谢参与,应助指数 +1
你的问题有关运行结果,所以需要你调用部分的程序。这样说说是无法判断的。
2楼2013-06-06 13:38:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

学虫dj

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
我感觉是你定义的全局变量的原因,为啥要定义那么多全局变量嫩
穷则独善其事,达则兼济天下
3楼2013-06-06 15:46:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sun805869434

金虫 (小有名气)

引用回帖:
3楼: Originally posted by 学虫dj at 2013-06-06 15:46:02
我感觉是你定义的全局变量的原因,为啥要定义那么多全局变量嫩

我有很多个子程序,这只是其中的两个,其他的子程序中能用到的我都定义了。要是全局变量的原因的话,为什么第一个能运行并得到正确的结果,而第二个就不能了呢?跟第二个程序中需要调用第一个运行结果有关吗?
优秀是一种习惯
4楼2013-06-06 18:37:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见