| 查看: 1832 | 回复: 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是最后的结果计算，谢谢 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有10人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】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
- 专业: 污染生态化学
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
xiaobear20
金虫 (小有名气)
(⊙o⊙)…
- 应助: 0 (幼儿园)
- 金币: 652.8
- 帖子: 106
- 在线: 25.8小时
- 虫号: 1252476
- 注册: 2011-04-01
- 性别: GG
- 专业: 传热传质学
5楼2012-04-07 15:34:19
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
- 专业: 污染生态化学
7楼2012-04-07 16:30:15
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
8楼2012-04-07 16:32:15
请问你解决了吗，我也是在调用子程序的时候，虽然子程序有做运算，但感觉最后的值没返回去。我也是每次循环的结果都一样，不知问题在哪里。 9楼2012-06-28 15:51:57
10楼2012-06-28 15:52:30
