 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1920 | 回復: 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ù)組），先給迭代的矩陣賦初值，繼而進行第一次迭代，每次迭代完后，把此次結(jié)果與上一次的結(jié)果進行比較，逐點比較，選出最大誤差，再與收斂條件ESP比較，如滿足條件則輸出結(jié)果，如不滿足，再返回迭代，直到滿足收斂條件為止，（其中也對迭代的次數(shù)進行了判斷） 我的問題如下：不知道怎么回事，程序只進行了一次計算就輸出了結(jié)果，沒有進行第二次，還有誤差判斷CONVER(I,J)總是為0，我把子程序的結(jié)果也輸出看了一下，各個子程序還是運行了一次。我就不知道問題在哪了。是結(jié)構(gòu)問題還是別的，求大神指道，分少別嫌棄啊，急等啊 子程序 CALL EDDY CALL CMOM CALL SOLV 在子程序EDDY中都用到了初始值，CMOM是矩陣系數(shù)的計算，SOLV是最后的結(jié)果計算，謝謝 |
回復此樓» 猜你喜歡
284求調(diào)劑
已經(jīng)有10人回復
-
一志愿山東大學藥學學碩求調(diào)劑
已經(jīng)有4人回復
-
07化學280分求調(diào)劑
已經(jīng)有4人回復
-
298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑
已經(jīng)有12人回復
-
求材料，環(huán)境專業(yè)調(diào)劑
已經(jīng)有3人回復
-
335求調(diào)劑
已經(jīng)有5人回復
-
求調(diào)劑
已經(jīng)有7人回復
-
一志愿吉大化學322求調(diào)劑
已經(jīng)有4人回復
-
環(huán)境學碩288求調(diào)劑
已經(jīng)有8人回復
-
341求調(diào)劑(一志愿湖南大學070300)
已經(jīng)有6人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 【求助】Fortran語言賦值問題？
已經(jīng)有3人回復
- 【求助】fortran編程問題
已經(jīng)有3人回復
- 【求助】調(diào)整系統(tǒng)時間后fortran報錯問題
已經(jīng)有7人回復
- 【求助】如何寫FORTRAN程序?qū)崿F(xiàn)求平均最近鄰距離
已經(jīng)有34人回復
- 【求助】用fortran怎么描述下面這樣的情況，我是新手，請大家?guī)兔?/a>
已經(jīng)有8人回復
- 【求助】Fortran編譯目標文件后的鏈接問題
已經(jīng)有4人回復
- 【求助】求指點，在Fortran里面，怎樣實現(xiàn)“數(shù)組維度可調(diào)”的數(shù)組?
已經(jīng)有5人回復
snoopyzhao
至尊木蟲 (職業(yè)作家)
- 程序強帖: 16
- 應(yīng)助: 157 (高中生)
- 貴賓: 0.02
- 金幣: 18844.7
- 紅花: 29
- 帖子: 3803
- 在線: 1422.4小時
- 蟲號: 183750
- 注冊: 2006-02-13
- 專業(yè): 污染生態(tài)化學
【答案】應(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 當然總是 0 了。因此，只迭代一次是很正常的。如果在 WOLD 之后，還有 W 更新，那么情況可能會不一樣。但這個沒有辦法看出來，因為你的所有 subroutine 的返回值是啥都看不出來…… 另外，你似乎喜歡用全局變量作為 subroutine 的參數(shù)和返回值，很多時候，這個習慣并不好，呵呵…… |
2樓2012-04-07 12:55:00
xiaobear20
金蟲 (小有名氣)
(⊙o⊙)…
- 應(yīng)助: 0 (幼兒園)
- 金幣: 652.8
- 帖子: 106
- 在線: 25.8小時
- 蟲號: 1252476
- 注冊: 2011-04-01
- 性別: GG
- 專業(yè): 傳熱傳質(zhì)學
3樓2012-04-07 13:35:37
snoopyzhao
至尊木蟲 (職業(yè)作家)
- 程序強帖: 16
- 應(yīng)助: 157 (高中生)
- 貴賓: 0.02
- 金幣: 18844.7
- 紅花: 29
- 帖子: 3803
- 在線: 1422.4小時
- 蟲號: 183750
- 注冊: 2006-02-13
- 專業(yè): 污染生態(tài)化學
4樓2012-04-07 14:37:59
xiaobear20
金蟲 (小有名氣)
(⊙o⊙)…
- 應(yīng)助: 0 (幼兒園)
- 金幣: 652.8
- 帖子: 106
- 在線: 25.8小時
- 蟲號: 1252476
- 注冊: 2011-04-01
- 性別: GG
- 專業(yè): 傳熱傳質(zhì)學
5樓2012-04-07 15:34:19
snoopyzhao
至尊木蟲 (職業(yè)作家)
- 程序強帖: 16
- 應(yīng)助: 157 (高中生)
- 貴賓: 0.02
- 金幣: 18844.7
- 紅花: 29
- 帖子: 3803
- 在線: 1422.4小時
- 蟲號: 183750
- 注冊: 2006-02-13
- 專業(yè): 污染生態(tài)化學
6樓2012-04-07 16:28:19
snoopyzhao
至尊木蟲 (職業(yè)作家)
- 程序強帖: 16
- 應(yīng)助: 157 (高中生)
- 貴賓: 0.02
- 金幣: 18844.7
- 紅花: 29
- 帖子: 3803
- 在線: 1422.4小時
- 蟲號: 183750
- 注冊: 2006-02-13
- 專業(yè): 污染生態(tài)化學
7樓2012-04-07 16:30:15
snoopyzhao
至尊木蟲 (職業(yè)作家)
- 程序強帖: 16
- 應(yīng)助: 157 (高中生)
- 貴賓: 0.02
- 金幣: 18844.7
- 紅花: 29
- 帖子: 3803
- 在線: 1422.4小時
- 蟲號: 183750
- 注冊: 2006-02-13
- 專業(yè): 污染生態(tài)化學
8樓2012-04-07 16:32:15
請問你解決了嗎，我也是在調(diào)用子程序的時候，雖然子程序有做運算，但感覺最后的值沒返回去。我也是每次循環(huán)的結(jié)果都一樣，不知問題在哪里。 9樓2012-06-28 15:51:57
10樓2012-06-28 15:52:30
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xiaobear20 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +10 | Zhao anqi 2026-03-22 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 求材料，環(huán)境專業(yè)調(diào)劑 +3 | 18567500178 2026-03-18 | 3/150 |
|
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +7 | 紙魚ly 2026-03-21 | 8/400 |
|
|
[考研] 考研化學308分求調(diào)劑 +7 | 你好明天你好 2026-03-23 | 8/400 |
|
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達問津 2026-03-23 | 5/250 |
|
|
|
[考研] 291 求調(diào)劑 +4 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 5/250 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-17 | 10/500 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 9/450 |
|
|
|
[考研] 北科281學碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
|
[考研] 材料學碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
|
[考研] 材料學碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
|
[考研] 083200學碩321分一志愿暨南大學求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿南昌大學，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研]
 295復試調(diào)劑
+8
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
|
[考研] 0703化學調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|