 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1258 | 回復(fù): 3 | ||
okagh鐵蟲 (初入文壇)
|
[求助]
菜鳥求助關(guān)于fortran if循環(huán)的問題，經(jīng)驗(yàn)豐富的前輩肯定一眼能看出問題出在哪
|
|
！ 我把源程序出問題的那部分發(fā)了上來，這個(gè)主要是一個(gè)if循環(huán)，clock在0~100和100~200區(qū)間程序一切正常debug也沒問題，但是在200~300區(qū)間Ac的輸出就不正常了，后來我debug模式調(diào)試，在200步之后發(fā)現(xiàn)循環(huán)是正常，但是在第十一行單步調(diào)試提示“無可用源”，程序運(yùn)行到第十二行時(shí)查看局部變量中的blf ，看到blf的值經(jīng)過十一行的運(yùn)算沒有改變，值還是上一個(gè)循環(huán)結(jié)尾的blf的值，因此在200~300區(qū)間程序把第十一行跳了過去，運(yùn)行時(shí)又不顯示出錯(cuò)，但結(jié)果是錯(cuò)的，請(qǐng)教大神能不能指出我的錯(cuò)誤在哪？ 一 DO 1000 CLOCK=1,10000 二 DO 100 i=((CLOCK-1)*10+1),CLOCK*10 三 if ((CLOCK.LT.100 .AND. CLOCK.GT.0) .OR. (CLOCK.LT.300 .AND. CLOCK.GT.200)) then！ 四 DO 110 j=1,mm 五 Aa(j)=Aa0*Ac(j) 六 ok(j)=Ke0*Ac(j)/Ad(j+50) 七 f(j)=ok(j)*dt1/(Aa(j)*dx*dx) 八 g(j)=g0*Ad(j+50)/Ac(j) 九 alf(j)=g(j)**2*dt1*dt1/(Ab*Aa(j)*(1+g(j)*dt1/Ab))-1-2*f(j)-g(j)*dt1/Aa(j) 十 alf(j)=alf(j)/f(j) 十一 blf(j)=-Ac(j)-(W/Aa(j))*exp(-alpha*j*dx)*exp(-4*0.69315*((i-11000)/1000.0)**2)-g(j)*dt1*Ad(j+50)/Aa(j) ！出問題的地方 出問題的地方 出問題的地方 十二 blf(j)=blf(j)/f(j) 110 CONTINUE blf(1)=blf(1)-Ac(1) blf(mm)=blf(mm)-Ac(mm) DO 120 k=2,mm alf(k)=alf(k)-1/alf(k-1) blf(k)=blf(k)-blf(k-1)/alf(k-1) 120 CONTINUE Ac(mm)=blf(mm)/alf(mm) DO 130 k=mm-1,1,-1 Ac(k)=(blf(k)-Ac(k+1))/alf(k) 130 CONTINUE else DO 140 j=1,mm Aa(j)=Aa0*Ac(j) ok(j)=Ke0*Ac(j)/Ad(j+50) f(j)=ok(j)*dt1/(Aa(j)*dx*dx) g(j)=g0*Ad(j+50)/Ac(j) alf(j)=g(j)**2*dt1*dt1/(Ab*Aa(j)*(1+g(j)*dt1/Ab))-1-2*f(j)-g(j)*dt1/Aa(j) alf(j)=alf(j)/f(j) blf(j)=-Ac(j)-g(j)*dt1*Ad(j+50)/Aa(j) blf(j)=blf(j)/f(j) 140 CONTINUE blf(1)=blf(1)-Ac(1) blf(mm)=blf(mm)-Ac(mm) DO 150 k=2,mm alf(k)=alf(k)-1/alf(k-1) blf(k)=blf(k)-blf(k-1)/alf(k-1) 150 CONTINUE Ac(mm)=blf(mm)/alf(mm) DO 160 k=mm-1,1,-1 Ac(k)=(blf(k)-Ac(k+1))/alf(k) 160 CONTINUE end if 100 CONTINUE write(*,*)Ac(5) 1000 CONTINUE |
回復(fù)此樓» 猜你喜歡
0856材料專碩353求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有6人回復(fù)
-
293求調(diào)劑
已經(jīng)有4人回復(fù)
-
306求0703調(diào)劑一志愿華中師范
已經(jīng)有10人回復(fù)
-
【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】
已經(jīng)有7人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有12人回復(fù)
-
材料與化工328分調(diào)劑
已經(jīng)有4人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有20人回復(fù)
-
一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271
已經(jīng)有6人回復(fù)
-
085404電子信息284分求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 菜鳥求助fortran數(shù)組聲明及子程序的一些問題啊
已經(jīng)有4人回復(fù)
- Fortran輸出二進(jìn)制數(shù)據(jù)，用matlab讀出，結(jié)果出錯(cuò)，求助問題何在？
已經(jīng)有4人回復(fù)
- 求助：運(yùn)行Fortran程序 提示出現(xiàn)的問題?
已經(jīng)有16人回復(fù)
- 菜鳥求助fortran編寫二階微分方程
已經(jīng)有5人回復(fù)
- fortran循環(huán)求助
已經(jīng)有3人回復(fù)
- 求助：intel visual fortran的格式問題
已經(jīng)有14人回復(fù)
- 求助：win7裝fortran軟件需要注意哪些問題，有沒有針對(duì)性的軟件？
已經(jīng)有6人回復(fù)
- FORTRAN新手 求助主程序循環(huán)問題
已經(jīng)有10人回復(fù)
- 用fortran遇到了一個(gè)問題，來求助了，謝謝先
已經(jīng)有3人回復(fù)
- 【求助】同樣一個(gè)問題用matlab和fortran編程計(jì)算，結(jié)果相差比較大，不知道為什么
已經(jīng)有4人回復(fù)
- 【求助】求高手關(guān)于Fortran數(shù)據(jù)讀寫的問題�！疽淹杲Y(jié)】
已經(jīng)有12人回復(fù)
- 【求助】一點(diǎn)FORTRAN的問題【已完結(jié)】
已經(jīng)有3人回復(fù)
- 【求助】Fortran編譯目標(biāo)文件后的鏈接問題
已經(jīng)有4人回復(fù)
qchem
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 280 (大學(xué)生)
- 金幣: 6871.6
- 散金: 1411
- 紅花: 43
- 帖子: 2586
- 在線: 669.2小時(shí)
- 蟲號(hào): 183011
- 注冊: 2006-02-12
- 專業(yè): 理論和計(jì)算化學(xué)
2樓2014-10-16 13:45:31
okagh
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 31
- 帖子: 18
- 在線: 24.9小時(shí)
- 蟲號(hào): 2819104
- 注冊: 2013-11-22
- 性別: GG
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
3樓2014-10-16 18:04:40
okagh
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 31
- 帖子: 18
- 在線: 24.9小時(shí)
- 蟲號(hào): 2819104
- 注冊: 2013-11-22
- 性別: GG
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
4樓2014-10-16 21:25:16
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 okagh 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
|
[考研] 327求調(diào)劑 +5 | prayer13 2026-03-23 | 5/250 |
|
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +3 | 你好明天你好 2026-03-23 | 3/150 |
|
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +5 | jiajunX 2026-03-22 | 5/250 |
|
|
|
[考研] 招08考數(shù)學(xué) +6 | laoshidan 2026-03-20 | 14/700 |
|
|
|
[考研]
求調(diào)劑材料學(xué)碩080500，總分289分
 5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +11 | 楓橋ZL 2026-03-18 | 13/650 |
|
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 268求調(diào)劑 +9 | 簡單點(diǎn)0 2026-03-17 | 9/450 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|