 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1937 | 回復(fù): 10 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
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 程序主要過程如下是，（計(jì)算的是一個(gè)二維網(wǎng)格點(diǎn)上的速度，即二維數(shù)組），先給迭代的矩陣賦初值，繼而進(jìn)行第一次迭代，每次迭代完后，把此次結(jié)果與上一次的結(jié)果進(jìn)行比較，逐點(diǎn)比較，選出最大誤差，再與收斂條件ESP比較，如滿足條件則輸出結(jié)果，如不滿足，再返回迭代，直到滿足收斂條件為止，（其中也對(duì)迭代的次數(shù)進(jìn)行了判斷） 我的問題如下：不知道怎么回事，程序只進(jìn)行了一次計(jì)算就輸出了結(jié)果，沒有進(jìn)行第二次，還有誤差判斷CONVER(I,J)總是為0，我把子程序的結(jié)果也輸出看了一下，各個(gè)子程序還是運(yùn)行了一次。我就不知道問題在哪了。是結(jié)構(gòu)問題還是別的，求大神指道，分少別嫌棄啊，急等啊 子程序 CALL EDDY CALL CMOM CALL SOLV 在子程序EDDY中都用到了初始值，CMOM是矩陣系數(shù)的計(jì)算，SOLV是最后的結(jié)果計(jì)算，謝謝 |
回復(fù)此樓» 猜你喜歡
一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑!
已經(jīng)有4人回復(fù)
-
298求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料學(xué)碩333求調(diào)劑
已經(jīng)有11人回復(fù)
-
0856求調(diào)劑
已經(jīng)有7人回復(fù)
-
356求調(diào)劑
已經(jīng)有3人回復(fù)
-
332求調(diào)劑
已經(jīng)有7人回復(fù)
-
348求調(diào)劑
已經(jīng)有4人回復(fù)
-
305求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑
已經(jīng)有13人回復(fù)
-
調(diào)劑310
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【求助】Fortran語言賦值問題？
已經(jīng)有3人回復(fù)
- 【求助】fortran編程問題
已經(jīng)有3人回復(fù)
- 【求助】調(diào)整系統(tǒng)時(shí)間后fortran報(bào)錯(cuò)問題
已經(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ù)
- 【求助】求指點(diǎn)，在Fortran里面，怎樣實(shí)現(xiàn)“數(shù)組維度可調(diào)”的數(shù)組?
已經(jīng)有5人回復(fù)
xiaobear20
金蟲 (小有名氣)
(⊙o⊙)…
- 應(yīng)助: 0 (幼兒園)
- 金幣: 652.8
- 帖子: 106
- 在線: 25.8小時(shí)
- 蟲號(hào): 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小時(shí)
- 蟲號(hào): 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
|
這個(gè)程序沒有辦法改，比如你的 W(I,J) 的值有更新嗎？在你的程序中，你的 WOLD 與 W 一直都是一樣的，那么 CONVER 當(dāng)然總是 0 了。因此，只迭代一次是很正常的。如果在 WOLD 之后，還有 W 更新，那么情況可能會(huì)不一樣。但這個(gè)沒有辦法看出來，因?yàn)槟愕乃?subroutine 的返回值是啥都看不出來…… 另外，你似乎喜歡用全局變量作為 subroutine 的參數(shù)和返回值，很多時(shí)候，這個(gè)習(xí)慣并不好，呵呵…… |
2樓2012-04-07 12:55:00
snoopyzhao
至尊木蟲 (職業(yè)作家)
- 程序強(qiáng)帖: 16
- 應(yīng)助: 157 (高中生)
- 貴賓: 0.02
- 金幣: 18844.7
- 紅花: 29
- 帖子: 3803
- 在線: 1422.4小時(shí)
- 蟲號(hào): 183750
- 注冊: 2006-02-13
- 專業(yè): 污染生態(tài)化學(xué)
4樓2012-04-07 14:37:59
xiaobear20
金蟲 (小有名氣)
(⊙o⊙)…
- 應(yīng)助: 0 (幼兒園)
- 金幣: 652.8
- 帖子: 106
- 在線: 25.8小時(shí)
- 蟲號(hào): 1252476
- 注冊: 2011-04-01
- 性別: GG
- 專業(yè): 傳熱傳質(zhì)學(xué)
5樓2012-04-07 15:34:19
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856求調(diào)劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研] 309求調(diào)劑 +6 | 誰不是少年 2026-03-29 | 6/300 |
|
|
[考研] 289求調(diào)劑 +13 | 新時(shí)代材料 2026-03-27 | 13/650 |
|
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+7
|
nnnnnnn5 2026-03-25 | 7/350 |
|
|
|
[考研] 295求調(diào)劑 +4 | wei-5 2026-03-26 | 4/200 |
|
|
|
[考研] 322求調(diào)劑 +7 | 宋明欣 2026-03-27 | 7/350 |
|
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4 | 小張zxy張 2026-03-26 | 8/400 |
|
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9 | 丹青奶蓋 2026-03-26 | 10/500 |
|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
|
[考研] 329求調(diào)劑 +5 | 1() 2026-03-22 | 5/250 |
|
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。�！ +4 | Acaciad 2026-03-25 | 4/200 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
|
[考研] 各位老師您好：本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|