| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 2473 | 回復(fù): 13 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
sun805869434金蟲(chóng) (小有名氣)
|
[求助]
fortran里一個(gè)子程序不運(yùn)行是怎么回事
|
|
|
整個(gè)fortran程序編譯完,運(yùn)行后就顯示內(nèi)存不能讀取,切換到Debug,按F5后,可以發(fā)現(xiàn)濃度值C(J,K)計(jì)算結(jié)果不對(duì),但是無(wú)論更改程序中該子程序的任何計(jì)算式,得到的結(jié)果都是一樣的,好像跟我設(shè)置的計(jì)算式?jīng)]有關(guān)系,這是沒(méi)有讀取相應(yīng)的子程序嗎?還是讀取其他部分了? SUBROUTINE COM !-------------------------------------------------------- PARAMETER(NJ=120,NT=120) COMMON U(NJ,NT),C(NJ,NT) COMMON DER(NT),Y(NJ),DYL(NJ) COMMON K,LITIM,DT COMMON DLN,DS,RHOF,RHOS,DD,GR,M1,YDN,DXL,YL,DY,C0 !-------------------------------------------------------- ENTRY VECTORY U(J,0)=0 DO K=1,LITIM DO J=1,M1 IF(J.NE.M1)THEN U(J,K)=(GR*DT+U(J,K-1))/(DD*DT+1) ELSE U(J,K)=0 ENDIF ENDDO ENDDO RETURN !------------------------------------------------------------ ENTRY CONCENTRATE C(J,0)=C0 DO K=1,LITIM DO J=1,M1 IF(J.EQ.1)THEN C(J,K)=(U(J,K)*C(J,K-1))/U(J,K+1) ELSEIF (J.GT.1.AND.J.LT.M1)THEN C(J,K)=((U(J,K)*C(J,K-1))/U(J,K+1))+C(J-1,K-1)-C(J-1,K) ELSE C(J,K)=0 ENDIF ENDDO ENDDO RETURN 這是子程序中的兩個(gè)部分,第一個(gè)子程序運(yùn)行后U(J,K)的結(jié)果是正確的,但是進(jìn)行第二個(gè)子程序時(shí),要用到上一個(gè)的結(jié)果U(J,K),運(yùn)行結(jié)果C(J,K)是錯(cuò)誤的,并且不管我怎么改C(J,K)的計(jì)算式,得出的結(jié)果不變,還是原來(lái)錯(cuò)誤的結(jié)果,不知道原因出在哪。 [ Last edited by sun805869434 on 2013-6-6 at 10:14 ] |

金蟲(chóng) (正式寫(xiě)手)

鐵桿木蟲(chóng) (著名寫(xiě)手)
工程和科學(xué)數(shù)值計(jì)算咨詢
金蟲(chóng) (正式寫(xiě)手)

金蟲(chóng) (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿上海海洋大學(xué)083200食品學(xué)碩,求調(diào)劑,接受其他專業(yè) +6 | what張 2026-04-01 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 材料專碩322分 +9 | 哈哈哈吼吼吼哈 2026-04-02 | 9/450 |
|
|
[考研] 08工科求調(diào)劑290分 +5 | 1314捧花 2026-04-02 | 8/400 |
|
|
[考研] 262求調(diào)劑 +5 | 勵(lì)志一定發(fā)文章 2026-04-02 | 6/300 |
|
|
[考研] 一志愿廈門(mén)大學(xué)材料工程專碩354找調(diào)劑。! +7 | 貝唄鋇鋇 2026-03-30 | 7/350 |
|
|
[考研] 材料工程322分 +8 | 哈哈哈吼吼吼哈 2026-04-01 | 8/400 |
|
|
[考研] 311求調(diào)劑 +16 | zchqwer 2026-04-01 | 18/900 |
|
|
[考研] 354求調(diào)劑 +4 | lxb598 2026-03-31 | 5/250 |
|
|
[考研] 805600專碩材料與化工348分求調(diào)劑 +6 | 上學(xué)啦! 2026-04-01 | 6/300 |
|
|
[考研] 材料求調(diào)劑 +8 | 呢呢妮妮 2026-04-01 | 8/400 |
|
|
[考研] 建環(huán),能源,土木老師路過(guò)看一看!! +4 | 嘿嘿uu 2026-04-01 | 4/200 |
|
|
[考研] 一志愿西安交大材料學(xué)碩(英一數(shù)二)347,求調(diào)劑到高分子/材料相關(guān)專業(yè) +7 | zju51 2026-03-31 | 9/450 |
|
|
[考研] 英一數(shù)二生物信息學(xué)287分,本科生物科學(xué),求調(diào)劑 +3 | 碧水xyz 2026-03-29 | 4/200 |
|
|
[考研] 一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +7 | @taotao 2026-03-30 | 7/350 |
|
|
[考研] 0856調(diào)劑 +7 | 曲聽(tīng)筠 2026-03-30 | 7/350 |
|
|
[考研] 299求調(diào)劑 +8 | 嗯嗯嗯嗯2 2026-03-27 | 8/400 |
|
|
[考研] 274求調(diào)劑 +6 | xiao愛(ài)同學(xué) 2026-03-30 | 6/300 |
|
|
[考研]
|
Gymno 2026-03-30 | 6/300 |
|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|
[考研] 266求調(diào)劑 +11 | 陽(yáng)陽(yáng)哇塞 2026-03-27 | 12/600 |
|