 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1942 | 回復(fù): 10 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
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ù)組），先給迭代的矩陣賦初值，繼而進(jìn)行第一次迭代，每次迭代完后，把此次結(jié)果與上一次的結(jié)果進(jìn)行比較，逐點比較，選出最大誤差，再與收斂條件ESP比較，如滿足條件則輸出結(jié)果，如不滿足，再返回迭代，直到滿足收斂條件為止，（其中也對迭代的次數(shù)進(jìn)行了判斷） 我的問題如下：不知道怎么回事，程序只進(jìn)行了一次計算就輸出了結(jié)果，沒有進(jìn)行第二次，還有誤差判斷CONVER(I,J)總是為0，我把子程序的結(jié)果也輸出看了一下，各個子程序還是運(yùn)行了一次。我就不知道問題在哪了。是結(jié)構(gòu)問題還是別的，求大神指道，分少別嫌棄啊，急等啊 子程序 CALL EDDY CALL CMOM CALL SOLV 在子程序EDDY中都用到了初始值，CMOM是矩陣系數(shù)的計算，SOLV是最后的結(jié)果計算，謝謝 |
回復(fù)此樓» 猜你喜歡
347求調(diào)劑
已經(jīng)有9人回復(fù)
-
070300化學(xué)求調(diào)劑
已經(jīng)有8人回復(fù)
-
070300一志愿211，312分求調(diào)劑院校
已經(jīng)有8人回復(fù)
-
考研調(diào)劑求助
已經(jīng)有3人回復(fù)
-
362求調(diào)劑
已經(jīng)有5人回復(fù)
-
085600 295分求調(diào)劑
已經(jīng)有11人回復(fù)
-
290分調(diào)劑求助
已經(jīng)有5人回復(fù)
-
311求調(diào)劑一志愿合肥工業(yè)大學(xué)
已經(jīng)有8人回復(fù)
-
085602 307分 求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南昌大學(xué)324求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 【求助】Fortran語言賦值問題？
已經(jīng)有3人回復(fù)
- 【求助】fortran編程問題
已經(jīng)有3人回復(fù)
- 【求助】調(diào)整系統(tǒng)時間后fortran報錯問題
已經(jīng)有7人回復(fù)
- 【求助】如何寫FORTRAN程序?qū)崿F(xiàn)求平均最近鄰距離
已經(jīng)有34人回復(fù)
- 【求助】用fortran怎么描述下面這樣的情況，我是新手，請大家?guī)兔?/a>
已經(jīng)有8人回復(fù)
- 【求助】Fortran編譯目標(biāo)文件后的鏈接問題
已經(jīng)有4人回復(fù)
- 【求助】求指點，在Fortran里面，怎樣實現(xiàn)“數(shù)組維度可調(diào)”的數(shù)組?
已經(jīng)有5人回復(fù)
snoopyzhao
至尊木蟲 (職業(yè)作家)
- 程序強(qiáng)帖: 16
- 應(yīng)助: 157 (高中生)
- 貴賓: 0.02
- 金幣: 18844.7
- 紅花: 29
- 帖子: 3803
- 在線: 1422.4小時
- 蟲號: 183750
- 注冊: 2006-02-13
- 專業(yè): 污染生態(tài)化學(xué)
7樓2012-04-07 16:30:15
snoopyzhao
至尊木蟲 (職業(yè)作家)
- 程序強(qiáng)帖: 16
- 應(yīng)助: 157 (高中生)
- 貴賓: 0.02
- 金幣: 18844.7
- 紅花: 29
- 帖子: 3803
- 在線: 1422.4小時
- 蟲號: 183750
- 注冊: 2006-02-13
- 專業(yè): 污染生態(tài)化學(xué)
【答案】應(yīng)助回帖
★
感謝參與，應(yīng)助指數(shù) +1
xzhdty: 金幣+1, 歡迎常來 2012-04-07 14:48:19
感謝參與，應(yīng)助指數(shù) +1
xzhdty: 金幣+1, 歡迎常來 2012-04-07 14:48:19
|
這個程序沒有辦法改，比如你的 W(I,J) 的值有更新嗎？在你的程序中，你的 WOLD 與 W 一直都是一樣的，那么 CONVER 當(dāng)然總是 0 了。因此，只迭代一次是很正常的。如果在 WOLD 之后，還有 W 更新，那么情況可能會不一樣。但這個沒有辦法看出來，因為你的所有 subroutine 的返回值是啥都看不出來…… 另外，你似乎喜歡用全局變量作為 subroutine 的參數(shù)和返回值，很多時候，這個習(xí)慣并不好，呵呵…… |
2樓2012-04-07 12:55:00
xiaobear20
金蟲 (小有名氣)
(⊙o⊙)…
- 應(yīng)助: 0 (幼兒園)
- 金幣: 652.8
- 帖子: 106
- 在線: 25.8小時
- 蟲號: 1252476
- 注冊: 2011-04-01
- 性別: GG
- 專業(yè): 傳熱傳質(zhì)學(xué)
3樓2012-04-07 13:35:37
snoopyzhao
至尊木蟲 (職業(yè)作家)
- 程序強(qiáng)帖: 16
- 應(yīng)助: 157 (高中生)
- 貴賓: 0.02
- 金幣: 18844.7
- 紅花: 29
- 帖子: 3803
- 在線: 1422.4小時
- 蟲號: 183750
- 注冊: 2006-02-13
- 專業(yè): 污染生態(tài)化學(xué)
4樓2012-04-07 14:37:59
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 調(diào)劑求院校招收 +7 | 鶴鯨鴿 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 070300化學(xué)279求調(diào)劑 +7 | 哈哈哈^_^ 2026-03-31 | 9/450 |
|
|
|
[考研] 0817化工學(xué)碩調(diào)劑 +7 | 努力上岸中！ 2026-03-31 | 7/350 |
|
|
|
[考研] 0856調(diào)劑 +6 | 曲聽筠 2026-03-30 | 6/300 |
|
|
|
[考研] 吉大生物學(xué)326分求調(diào)劑 +3 | sunnyupup 2026-03-31 | 3/150 |
|
|
|
[考研]
 一志愿大連理工大學(xué)材料求調(diào)劑
+6
|
Gymno 2026-03-30 | 6/300 |
|
|
|
[考研] 本科211總分289，08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
[考研] 297求調(diào)劑 +17 | 田洪有 2026-03-26 | 18/900 |
|
|
|
[考研] 286求調(diào)劑 +5 | Faune 2026-03-30 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
|
[考研] 生物技術(shù)與工程 +7 | 1294608413 2026-03-25 | 8/400 |
|
|
|
[考研] 349求調(diào)劑 +6 | 李木子啊哈哈 2026-03-25 | 6/300 |
|
|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿吉大071010，316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
|
[考研]
材料調(diào)劑
 5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|