| 查看: 1756 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
xiaobear20金虫 (小有名气)
(⊙o⊙)…
|
[求助]
FORTRAN新手 求助主程序循环问题
|
|
|
我现贴出主程序部分： PROGRAM MAIN USE SAME IMPLICIT NONE N=120 M=241 ITMAX=60 ITER=0 CALL INPUT CALL GRID DO I=1,M !!!!!!初始条件!!!!!!!!!! DO J=2,N-1 W(I,J)=1-EXP(-ETA(J)) ENDDO ENDDO 60 ITER=ITER+1 SELECT CASE(ITER-ITMAX) CASE(:-1) ! <= -1 整数 GO TO 70 CASE(0) ! =0 WRITE(6,2500) GO TO 70 CASE(1  ! >=1WRITE(6,2600) GO TO 90 END SELECT 70 DO I=1,M DO J=2,N-1 WOLD(I,J)=W(I,J) ENDDO ENDDO CALL EDDY CALL CMOM CALL SOLV DO I=2,M DO J=2,N-1 CONVER(I,J)=ABS(1-(WOLD(I,J)/W(I,J))) ENDDO ENDDO MX=MAXVAL(CONVER) ！选取误差中的最大值 IF(MX.LE.EPS)THEN ! <=0.0001 !收敛判断 EPS=0.0001 误差 CONVER GO TO 90 ELSE GOTO 60 ENDIF 90 CALL OUTPUT END PROGRAM 程序主要过程如下是，（计算的是一个二维网格点上的速度，即二维数组），先给迭代的矩阵赋初值，继而进行第一次迭代，每次迭代完后，把此次结果与上一次的结果进行比较，逐点比较，选出最大误差，再与收敛条件ESP比较，如满足条件则输出结果，如不满足，再返回迭代，直到满足收敛条件为止，（其中也对迭代的次数进行了判断） 我的问题如下：不知道怎么回事，程序只进行了一次计算就输出了结果，没有进行第二次，还有误差判断CONVER(I,J)总是为0，我把子程序的结果也输出看了一下，各个子程序还是运行了一次。我就不知道问题在哪了。是结构问题还是别的，求大神指道，分少别嫌弃啊，急等啊 子程序 CALL EDDY CALL CMOM CALL SOLV 在子程序EDDY中都用到了初始值，CMOM是矩阵系数的计算，SOLV是最后的结果计算，谢谢 |
回复此楼» 猜你喜欢
求个博导看看
已经有18人回复
-
自荐读博
已经有6人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有5人回复
-
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】Fortran语言赋值问题？
已经有3人回复
- 【求助】fortran编程问题
已经有3人回复
- 【求助】调整系统时间后fortran报错问题
已经有7人回复
- 【求助】如何写FORTRAN程序实现求平均最近邻距离
已经有34人回复
- 【求助】用fortran怎么描述下面这样的情况，我是新手，请大家帮忙
已经有8人回复
- 【求助】Fortran编译目标文件后的链接问题
已经有4人回复
- 【求助】求指点，在Fortran里面，怎样实现“数组维度可调”的数组?
已经有5人回复
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
6楼2012-04-07 16:28:19
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
2楼2012-04-07 12:55:00
xiaobear20
金虫 (小有名气)
(⊙o⊙)…
- 应助: 0 (幼儿园)
- 金币: 652.8
- 帖子: 106
- 在线: 25.8小时
- 虫号: 1252476
- 注册: 2011-04-01
- 性别: GG
- 专业: 传热传质学
3楼2012-04-07 13:35:37
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
4楼2012-04-07 14:37:59