| 查看: 1921 | 回復: 10 | ||
xiaobear20金蟲 (小有名氣)
(⊙o⊙)…
|
[求助]
FORTRAN新手 求助主程序循環(huán)問題
|
|
我現(xiàn)貼出主程序部分: 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 整數(shù) 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 程序主要過程如下是,(計算的是一個二維網(wǎng)格點上的速度,即二維數(shù)組),先給迭代的矩陣賦初值,繼而進行第一次迭代,每次迭代完后,把此次結果與上一次的結果進行比較,逐點比較,選出最大誤差,再與收斂條件ESP比較,如滿足條件則輸出結果,如不滿足,再返回迭代,直到滿足收斂條件為止,(其中也對迭代的次數(shù)進行了判斷) 我的問題如下:不知道怎么回事,程序只進行了一次計算就輸出了結果,沒有進行第二次,還有誤差判斷CONVER(I,J)總是為0,我把子程序的結果也輸出看了一下,各個子程序還是運行了一次。我就不知道問題在哪了。是結構問題還是別的,求大神指道,分少別嫌棄啊,急等啊 子程序 CALL EDDY CALL CMOM CALL SOLV 在子程序EDDY中都用到了初始值,CMOM是矩陣系數(shù)的計算,SOLV是最后的結果計算,謝謝 |

至尊木蟲 (職業(yè)作家)
金蟲 (小有名氣)
(⊙o⊙)…

至尊木蟲 (職業(yè)作家)
金蟲 (小有名氣)
(⊙o⊙)…

至尊木蟲 (職業(yè)作家)
至尊木蟲 (職業(yè)作家)
至尊木蟲 (職業(yè)作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料學碩333求調劑 +3 | 北道巷 2026-03-24 | 3/150 |
|
|---|---|---|---|---|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 5/250 |
|
|
[考研] 求調劑 +6 | 研研,接電話 2026-03-24 | 7/350 |
|
|
[考研] 材料與化工328分調劑 +4 | 。,。,。,。i 2026-03-23 | 4/200 |
|
|
[考研] 276求調劑。有半年電池和半年高分子實習經(jīng)歷 +9 | 材料學257求調劑 2026-03-23 | 10/500 |
|
|
[考研] 291求調劑 +8 | hhhhxn.. 2026-03-23 | 8/400 |
|
|
[考研] 303求調劑 +4 | 元夕元 2026-03-20 | 4/200 |
|
|
[考研] 一志愿西安交通大學材料工程專業(yè) 282分求調劑 +11 | 楓橋ZL 2026-03-18 | 13/650 |
|
|
[考研] 275求調劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 318求調劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
[考研] 求調劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
[考研] 303求調劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 材料求調劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
[考研] 材料 271求調劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 一志愿西南交大,求調劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
[考研] 308求調劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研] 330求調劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 304求調劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 一志愿武漢理工材料工程專碩調劑 +9 | Doleres 2026-03-19 | 9/450 |
|