| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1922 | 回復(fù): 10 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
xiaobear20金蟲(chóng) (小有名氣)
(⊙o⊙)…
|
[求助]
FORTRAN新手 求助主程序循環(huán)問(wè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 程序主要過(guò)程如下是,(計(jì)算的是一個(gè)二維網(wǎng)格點(diǎn)上的速度,即二維數(shù)組),先給迭代的矩陣賦初值,繼而進(jìn)行第一次迭代,每次迭代完后,把此次結(jié)果與上一次的結(jié)果進(jìn)行比較,逐點(diǎn)比較,選出最大誤差,再與收斂條件ESP比較,如滿足條件則輸出結(jié)果,如不滿足,再返回迭代,直到滿足收斂條件為止,(其中也對(duì)迭代的次數(shù)進(jìn)行了判斷) 我的問(wèn)題如下:不知道怎么回事,程序只進(jìn)行了一次計(jì)算就輸出了結(jié)果,沒(méi)有進(jìn)行第二次,還有誤差判斷CONVER(I,J)總是為0,我把子程序的結(jié)果也輸出看了一下,各個(gè)子程序還是運(yùn)行了一次。我就不知道問(wèn)題在哪了。是結(jié)構(gòu)問(wèn)題還是別的,求大神指道,分少別嫌棄啊,急等啊 子程序 CALL EDDY CALL CMOM CALL SOLV 在子程序EDDY中都用到了初始值,CMOM是矩陣系數(shù)的計(jì)算,SOLV是最后的結(jié)果計(jì)算,謝謝 |

至尊木蟲(chóng) (職業(yè)作家)
|
這個(gè)程序沒(méi)有辦法改,比如你的 W(I,J) 的值有更新嗎?在你的程序中,你的 WOLD 與 W 一直都是一樣的,那么 CONVER 當(dāng)然總是 0 了。因此,只迭代一次是很正常的。如果在 WOLD 之后,還有 W 更新,那么情況可能會(huì)不一樣。但這個(gè)沒(méi)有辦法看出來(lái),因?yàn)槟愕乃?subroutine 的返回值是啥都看不出來(lái)…… 另外,你似乎喜歡用全局變量作為 subroutine 的參數(shù)和返回值,很多時(shí)候,這個(gè)習(xí)慣并不好,呵呵…… |
金蟲(chóng) (小有名氣)
(⊙o⊙)…

至尊木蟲(chóng) (職業(yè)作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +6 | 小羊睡著了? 2026-03-23 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +8 | 碩星赴 2026-03-23 | 8/400 |
|
|
[考研] 收08調(diào)劑生 +6 | komorebi69 2026-03-18 | 6/300 |
|
|
[考研] 材料專(zhuān)碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
[考研] 一志愿陜師大生物學(xué)071000,298分,求調(diào)劑 +3 | SYA! 2026-03-23 | 3/150 |
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
[考研] 317求調(diào)劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 318求調(diào)劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
[考研] 一志愿深大,0703化學(xué),總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|