| 查看: 2293 | 回复: 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 ] |
回复此楼» 猜你喜欢
论文投稿，期刊推荐
已经有6人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有3人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
sun805869434
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 173.5
- 散金: 150
- 红花: 2
- 帖子: 57
- 在线: 26.9小时
- 虫号: 1666767
- 注册: 2012-03-05
- 性别: MM
- 专业: 传热传质学
12楼2013-06-14 09:04:45
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
2楼2013-06-06 13:38:25
学虫dj
金虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 892.4
- 散金: 63
- 帖子: 688
- 在线: 89.9小时
- 虫号: 1829966
- 注册: 2012-05-22
- 性别: GG
- 专业: 应用地球物理学
3楼2013-06-06 15:46:02
sun805869434
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 173.5
- 散金: 150
- 红花: 2
- 帖子: 57
- 在线: 26.9小时
- 虫号: 1666767
- 注册: 2012-03-05
- 性别: MM
- 专业: 传热传质学
4楼2013-06-06 18:37:32