版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

bohdshen

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 275.7
帖子: 14
在線: 4.7小時
蟲號: 918750
注冊: 2009-12-02
專業(yè): 應(yīng)用地球物理學(xué)

[交流] 【求助】求助程序我問題解決送80金幣已有3人參與

這個非常簡單是我太笨請幫忙看下如我寫有問題請大俠幫我改寫
輸入DEN AC已知。計算POR，初值PORI=0.029
公式如下：

問題：只能計算一小段的POR 再往下就不能算函數(shù)關(guān)系表達式無誤
提示錯誤：
run-time error M6201: MATH
- **: DOMAIN error
未加粗體是該程序的頭尾固定部分，不用改寫。粗體部分為我加的程序。

   BLOCK DATA
   IMPLICIT INTEGER*4(I-N)
   INTEGER*2 NUMLOG,IZ, IDC, IZ2
   Character*4 NameI,NameO,NameC
   COMMON /INPC/NOI,NAMEI(2)
   COMMON /OUTC/NOO,NAMEO(1)
   COMMON /CONC/NOC,NAMEC(1)
   COMMON /INP /AC  ,DEN
   COMMON /OUTP/POR
   COMMON /CON /ZZ(1)
   COMMON /UNIT/IU,JU,KU,LU,MU
   COMMON /HD  /DEP,DEPO,SDEP,EDEP,ECC,ICO(20),IWELL(20),
   & NUMLOG,IZ,LOG1(20),STDEP,ENDEP,RLEV,IDC,IZ2,FEET,
   & LOG2(20),IDUM1(118)
   DATA IU,JU,KU,LU,MU/1,2,4,3,0/
   DATA NoI,NoO,NoC/2,1,1/
   DATA NAMEI/'AC  ','DEN '/
   DATA NAMEO/'POR '/
   DATA NAMEC/'PORI'/
   DATA ZZ /0.029/
   END

   Program MyProgram

   IMPLICIT INTEGER*4(I-N)
   INTEGER*2 NUMLOG,IZ, IDC, IZ2
   Character*4 NameI,NameO,NameC
   COMMON /INPC/NOI,NAMEI(2)
   COMMON /OUTC/NOO,NAMEO(1)
   COMMON /CONC/NOC,NAMEC(1)
   COMMON /INP /AC  ,DEN
   COMMON /OUTP/POR
   COMMON /CON /PORI
   COMMON /INPB/BUFI(25600)
   COMMON /OUTB/BUFO(25600)
   COMMON /UNIT/IU,JU,KU,LU,MU
   COMMON /HD  /DEP,DEPO,SDEP,EDEP,ECC,ICO(20),IWELL(20),
   & NUMLOG,IZ,LOG1(20),STDEP,ENDEP,RLEV,IDC,IZ2,FEET,
   & LOG2(20),IDUM1(118)
   CALL PART
   CALL RDFLNM

   CALL CONST
   CALL IN
   CALL OUT
1    CALL IN
   AS=92414*DEN/(AC**2)
   CALL WASAT(AS,PORI,POR)


   CALL OUT
   GOTO 1
   END
   SUBROUTINE WASAT(AS,PORI,POR)
9    AA=(1-PORI)**(3/(1-PORI))
   AM=10.17*((1-AA)**2)
   AN=AS-47.317*AA
   AT=(AM/AN)-PORI
   IF(ABS(AT) .LE. 0.001) THEN
   POR=PORI*100
   ELSE
   PORI=PORI+0.002
   GOTO 9
   ENDIF
   RETURN
   END

回復(fù)此樓

» 猜你喜歡

299求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑已經(jīng)有6人回復(fù)
300求調(diào)劑，材料科學(xué)英一數(shù)二已經(jīng)有8人回復(fù)
招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑已經(jīng)有5人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有7人回復(fù)
289求調(diào)劑已經(jīng)有13人回復(fù)
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分已經(jīng)有9人回復(fù)
081200-11408-276學(xué)碩求調(diào)劑已經(jīng)有6人回復(fù)
調(diào)劑求院校招收已經(jīng)有5人回復(fù)
調(diào)劑310 已經(jīng)有8人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【求助】畢業(yè)論文的問題（有關(guān)酰胺合成）-能解決愿傾家蕩產(chǎn)發(fā)金幣已經(jīng)有31人回復(fù)
【求助】求代碼對應(yīng)的物質(zhì)結(jié)構(gòu)，解決問題金幣全送已經(jīng)有6人回復(fù)
【求助】關(guān)于win7睡眠設(shè)置的問題，謝謝。有金幣送哦~~ 已經(jīng)有10人回復(fù)

1樓 2010-05-09 20:51:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

★
resonant(金幣+1):感謝回復(fù)：-） 2010-05-10 13:54:40

你的公式與你的子程序表達的似乎不是一樣的東西呀……

你是求那個公式的值還是別的？

贊一下(1人)

回復(fù)此樓

2樓2010-05-10 12:46:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

★
resonant(金幣+1):感謝建設(shè)性指導(dǎo)：-） 2010-05-10 14:11:18

我猜問題可能出在 WASAT 這個子程中，也就是 PORI = PORI + 0.002 這一步，如果 PORI 等于或大于 1 時， AA 這個表達式還有意義嗎？

在這一句后面加個判斷好了，比如:

IF (PORI .GE. 1.0) STOP

另外，你的 WASAT 這個子程中，修改一下下面兩句話：

   AA=(1.0-PORI)**(3.0/(1.0-PORI))
   AM=10.17*((1.0-AA)**2.0)

贊一下(1人)

回復(fù)此樓

3樓2010-05-10 14:00:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bohdshen

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 275.7
帖子: 14
在線: 4.7小時
蟲號: 918750
注冊: 2009-12-02
專業(yè): 應(yīng)用地球物理學(xué)

nono2009:建議通過“引用回復(fù)該帖”，以便別人收到你的message. 2010-05-10 15:12:05

您好
這個PORI是初始孔隙度
最后POR應(yīng)該在0~1之間
我做的這段深度內(nèi)孔隙度POR實際值已經(jīng)知道在0.09到0.25之間
我開是沒用wasat直接寫問題也一樣是做了1米就不做了
無加深處是FoRWARD測井軟件的固定部分
它是順序結(jié)構(gòu) 逐點深度計算后來別人說我直接寫怕有沖突就讓我用wasat
我現(xiàn)在因為下午要去開會
所以晚上來探討您的回復(fù)
謝謝

贊一下(1人)

回復(fù)此樓

4樓2010-05-10 14:53:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bohdshen

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 275.7
帖子: 14
在線: 4.7小時
蟲號: 918750
注冊: 2009-12-02
專業(yè): 應(yīng)用地球物理學(xué)

resonant:nano區(qū)長已經(jīng)提示過，不過朋友沒有完全理解，必須引用你想要通知的人的言論，才是有效的。引用自己的發(fā)言對方還是不會知道... 2010-05-10 21:06:39

引用回帖:

Originally posted by bohdshen at 2010-05-10 14:53:23:
您好
這個PORI是初始孔隙度
最后POR應(yīng)該在0~1之間
我做的這段深度內(nèi)孔隙度POR實際值已經(jīng)知道在0.09到0.25之間
我開是沒用wasat直接寫問題也一樣是做了1米就不做了
無加深處是FoRWARD測井軟件的固定部分
它 ...

首先感謝你
我試驗了你的方法加了stop
確實運行了但是結(jié)果除了開始的1米沒有問題
下面的結(jié)果都是0
這說明是真的PORI超出了1了吧。。
可是這段范圍在結(jié)論上應(yīng)該在0.03到0.3之間才對
這我該怎么修改呢
敬請指教非常之感謝

贊一下(1人)

回復(fù)此樓

5樓2010-05-10 19:46:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

★
余澤成(金幣+1):謝謝應(yīng)助！ 2010-05-10 20:25:21

你把 WASAT 需要的參數(shù)給出來，你現(xiàn)在給出的程序沒有辦法運行啊……或者你給一個可以運行的，可以反映你的問題的程序……

贊一下(1人)

回復(fù)此樓

6樓2010-05-10 19:59:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bohdshen

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 275.7
帖子: 14
在線: 4.7小時
蟲號: 918750
注冊: 2009-12-02
專業(yè): 應(yīng)用地球物理學(xué)

引用回帖:

Originally posted by snoopyzhao at 2010-05-10 19:59:23:
你把 WASAT 需要的參數(shù)給出來，你現(xiàn)在給出的程序沒有辦法運行啊……或者你給一個可以運行的，可以反映你的問題的程序……

我給你傳數(shù)據(jù)和說明公式，這樣可能方便些
qq798858957,或者傳給你的郵箱
十分感謝

贊一下

回復(fù)此樓

7樓2010-05-10 21:26:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

★
resonant(金幣+1):謝謝。 2010-05-10 23:02:23

我給你的信息中有我的信箱。俺沒有 QQ……

贊一下(1人)

回復(fù)此樓

8樓2010-05-10 21:31:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

★ ★
wangen994(金幣+2):非常熱心，辛苦了。 2010-05-11 08:42:01

看到了你的數(shù)據(jù)了。由于時間的關(guān)系，還沒有來得及仔細(xì)研究。但初步的計算結(jié)果表明，AT值過小，以致于在0.029~1這樣的區(qū)間里，所有的值都不能滿足 AT 這樣的限制要求。

換一種思路看看如何：

如果你確信 POR 在什么樣的范圍之內(nèi)，那么以非常小的步長，遍歷整個區(qū)間，然后以 AT 最小的那個 POR 為你的結(jié)果。

贊一下(1人)

回復(fù)此樓

9樓2010-05-11 07:48:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

bohdshen(金幣+88):感謝你的付出 2010-05-11 10:28:05

下面是修改后的程序，小寫字母部分是我改寫的。主要的目的在于打印出符合要求的 POR, 對于不能獲得的 POR，給出一個提示。

zzz 是一個包括 AC 和 DEN 兩列數(shù)據(jù)的文件。

CODE:

      program main

      open(10,file='zzz',status='old')

      do 5 i=1,509

        read(10,*) ac, den

        pori = 0.029

        AS=92414*DEN/(AC**2)

        CALL WASAT(AS,PORI,POR, ind)

        if(ind .ne. 1) write(*,*) por

5     continue

      END

      SUBROUTINE WASAT(AS,PORI,POR, ind)

      ind = 0

9     AA=(1.0-PORI)**(3.0/(1.0-PORI))

      AM=10.17*((1.0-AA)**2.0)

      AN=AS-47.317*AA

      AT=(AM/AN)-PORI

      IF(ABS(AT) .LE. 0.001) THEN

        POR=PORI*100

      ELSE

        pori=pori+0.000001

        if (pori .ge. .3) then

          write(*,*) 'i can not get the correct value!'

          ind = 1

          return

        endif

        GOTO 9

      ENDIF                                                        

      RETURN

      END

贊一下

回復(fù)此樓

10樓2010-05-11 09:56:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 bohdshen 的主題更新

返回列表

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +6	mikasa的圍巾 2026-03-28	6/300	2026-03-29 01:32 by fmesaito
[考研] 289求調(diào)劑 +13	新時代材料 2026-03-27	13/650	2026-03-29 01:16 by 544594351
[考研] 322求調(diào)劑 +7	宋明欣 2026-03-27	7/350	2026-03-28 21:27 by sanrepian
[考研] 壓國家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7	迷人的哈哈 2026-03-28	7/350	2026-03-28 16:47 by 催化大白
[考研] 312，生物學(xué)求調(diào)劑 +3	小譯同學(xué)abc 2026-03-28	3/150	2026-03-28 15:32 by 落�？伤�
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +3	nhj_ 2026-03-25	3/150	2026-03-28 13:24 by Iveryant
[考研] 085600 286分材料求調(diào)劑 +7	麻辣魷魚 2026-03-27	8/400	2026-03-28 12:17 by zllcz
[考研] 材料與化工考研調(diào)劑 +17	孅華 2026-03-22	17/850	2026-03-28 08:35 by WYUMater
[考研] 311求調(diào)劑 +3	希望上岸阿小楊 2026-03-23	3/150	2026-03-28 07:57 by 熱情沙漠
[考研] 材料求調(diào)劑一志愿哈工大總分298分，前三科223分 +5	dongfang59 2026-03-27	5/250	2026-03-28 04:53 by wxiongid
[考研] 考研調(diào)劑 +4	Sanmu-124 2026-03-26	4/200	2026-03-27 17:49 by kiokin
[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +6	kotoko_ik 2026-03-23	7/350	2026-03-27 12:29 by 惠州彭于晏
[考研] 調(diào)劑推薦 +5	清酒714 2026-03-26	6/300	2026-03-27 11:12 by 不吃魚的貓
[考研] 一志愿吉大071010，316分求調(diào)劑 +3	xgbiknn 2026-03-27	3/150	2026-03-27 10:36 by guoweigw
[考研] 329求調(diào)劑 +7	鈕恩雪 2026-03-25	7/350	2026-03-27 04:28 by wxiongid
[考研] 303求調(diào)劑 +6	藍山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3	min3 2026-03-24	4/200	2026-03-25 18:22 by xcjcqu
[考研] 一志愿武理085500機械專業(yè)總分300求調(diào)劑 +3	an10101 2026-03-24	7/350	2026-03-25 00:00 by 山鬼0-
[考研] 080500求調(diào)劑 +3	zzzzfan 2026-03-24	3/150	2026-03-24 16:38 by barlinike
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518

24小時熱門版塊排行榜

[交流] 【求助】求助 程序我問題 解決送80金幣 已有3人參與

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

[交流] 【求助】求助程序我問題解決送80金幣已有3人參與

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: