| 查看: 2320 | 回复: 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 ] |
回复此楼» 猜你喜欢
最近几年招的学生写论文不引自己组发的文章
已经有4人回复
-
职称评审没过，求安慰
已经有54人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- fortran 程序有write和没有write两种情况计算的结果为什么会不一样?
已经有14人回复
- 一百个FORTRAN任务怎么能一次提交SSH计算
已经有3人回复
- 帮忙把公式转化为FORTRAN子程序
已经有6人回复
- 一个主程序，一个subroutine程序，要如何将两个程序组合在fortran里运行
已经有4人回复
- 又来求助啦，fortran里使用多文件，其中的子程序文件怎么观察参数值啊？
已经有3人回复
- 请教 fortran 运行错误的原因
已经有13人回复
- FORTRAN新手 求助主程序循环问题
已经有10人回复
- 大家帮我看一个fortran的程序，我总是计算不出正确的结果
已经有20人回复
- FORTRAN 基础知识讲解
已经有21人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- 有关fortran的一次而问题，希望大家能帮帮忙，谢谢
已经有4人回复
- 【求助】f90子程序怎么编译错误？请大家帮忙，谢谢
已经有3人回复
- 【求助】Intel Visual Fortran与Compaq Visual Fortran的区别？
已经有18人回复
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
8楼2013-06-06 21:27:39
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
10楼2013-06-07 09:28:31
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
sun805869434
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 173.5
- 散金: 150
- 红花: 2
- 帖子: 57
- 在线: 26.9小时
- 虫号: 1666767
- 注册: 2012-03-05
- 性别: MM
- 专业: 传热传质学
5楼2013-06-06 18:50:26
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
6楼2013-06-06 20:32:17
sun805869434
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 173.5
- 散金: 150
- 红花: 2
- 帖子: 57
- 在线: 26.9小时
- 虫号: 1666767
- 注册: 2012-03-05
- 性别: MM
- 专业: 传热传质学
7楼2013-06-06 21:01:21
学虫dj
金虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 892.4
- 散金: 63
- 帖子: 688
- 在线: 89.9小时
- 虫号: 1829966
- 注册: 2012-05-22
- 性别: GG
- 专业: 应用地球物理学
9楼2013-06-07 08:45:29