应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】求助 程序我问题 解决送80金币
北京石油化工学院2026年研究生招生接收调剂公告
151/2返回列表12下一页
查看: 1299  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

bohdshen

银虫 (初入文坛)

[交流] 【求助】求助 程序我问题 解决送80金币 已有3人参与

这个非常简单是我太笨 请帮忙看下  如我写有问题  请大侠帮我改写
输入DEN AC已知。 计算POR，初值PORI=0.029
公式如下：
问题：只能计算一小段的POR 再往下就不能算 函数关系表达式无误
提示错误：
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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-05-09 20:51:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★
resonant(金币+1):感谢回复：-） 2010-05-10 13:54:40
你的公式与你的子程序表达的似乎不是一样的东西呀……

你是求那个公式的值还是别的？
赞一下(1人) 回复此楼
2楼2010-05-10 12:46:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★
resonant(金币+1):感谢建设性指导：-） 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人) 回复此楼
3楼2010-05-10 14:00:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bohdshen

银虫 (初入文坛)
nono2009:建议通过“引用回复该帖”，以便别人收到你的message. 2010-05-10 15:12:05
您好
这个PORI是初始孔隙度
最后POR应该在0~1之间
我做的这段深度内孔隙度POR实际值已经知道在0.09到0.25之间
我开是没用wasat直接写 问题也一样是做了1米就不做了
无加深处是FoRWARD测井软件的固定部分
它是顺序结构 逐点深度计算 后来别人说我直接写怕有冲突 就让我用wasat
我现在因为下午要去开会
所以晚上来探讨您的回复
谢谢
赞一下(1人) 回复此楼
4楼2010-05-10 14:53:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bohdshen

银虫 (初入文坛)
resonant:nano区长已经提示过，不过朋友没有完全理解，必须引用你想要通知的人的言论，才是有效的。引用自己的发言对方还是不会知道... 2010-05-10 21:06:39
引用回帖:
Originally posted by bohdshen at 2010-05-10 14:53:23:
您好
这个PORI是初始孔隙度
最后POR应该在0~1之间
我做的这段深度内孔隙度POR实际值已经知道在0.09到0.25之间
我开是没用wasat直接写 问题也一样是做了1米就不做了
无加深处是FoRWARD测井软件的固定部分
它 ...

首先感谢你
我试验了你的方法 加了stop
确实运行了 但是结果除了开始的1米没有问题
下面的结果都是0
这说明是真的PORI超出了1了吧。。
可是这段范围在结论上应该在0.03到0.3之间才对
这我该怎么修改呢
敬请指教 非常之感谢
赞一下(1人) 回复此楼
5楼2010-05-10 19:46:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★
余泽成(金币+1):谢谢应助！ 2010-05-10 20:25:21
你把 WASAT 需要的参数给出来，你现在给出的程序没有办法运行啊……或者你给一个可以运行的，可以反映你的问题的程序……
赞一下(1人) 回复此楼
6楼2010-05-10 19:59:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bohdshen

银虫 (初入文坛)
引用回帖:
Originally posted by snoopyzhao at 2010-05-10 19:59:23:
你把 WASAT 需要的参数给出来，你现在给出的程序没有办法运行啊……或者你给一个可以运行的，可以反映你的问题的程序……

我给你传数据和说明公式，这样可能方便些
qq798858957,或者传给你的邮箱
十分感谢
赞一下 回复此楼
7楼2010-05-10 21:26:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★
resonant(金币+1):谢谢。 2010-05-10 23:02:23
我给你的信息中有我的信箱。俺没有 QQ……
赞一下(1人) 回复此楼
8楼2010-05-10 21:31:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★ ★
wangen994(金币+2):非常热心，辛苦了。 2010-05-11 08:42:01
看到了你的数据了。由于时间的关系，还没有来得及仔细研究。但初步的计算结果表明，AT值过小，以致于在0.029~1这样的区间里，所有的值都不能满足 AT 这样的限制要求。

换一种思路看看如何：

如果你确信 POR 在什么样的范围之内，那么以非常小的步长，遍历整个区间，然后以 AT 最小的那个 POR 为你的结果。
赞一下(1人) 回复此楼
9楼2010-05-11 07:48:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
bohdshen(金币+88):感谢你的付出 2010-05-11 10:28:05
下面是修改后的程序，小写字母部分是我改写的。主要的目的在于打印出符合要求的 POR, 对于不能获得的 POR，给出一个提示。

zzz 是一个包括 AC 和 DEN 两列数据的文件。
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
赞一下 回复此楼
10楼2010-05-11 09:56:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bohdshen 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 华南理工0703化学，总分336求调剂 +5 手机用户 2026-04-03 5/250 2026-04-06 09:11 by 无际的草原
[考研] 0703化学调剂325分 +10 15771691647 2026-04-04 10/500 2026-04-06 06:31 by houyaoxu
[考研] 070300化学学硕311分求调剂 +9 梁富贵险中求 2026-04-04 9/450 2026-04-06 06:24 by houyaoxu
[考研] 262求调剂 +7 天下第一文 2026-04-04 8/400 2026-04-05 21:31 by 激流勇渡
[考研] 求调剂 一志愿西南交通大学085701环境工程 282分 +7 多多爱吃汉堡 2026-04-04 7/350 2026-04-05 19:47 by 福农资环_环境基
[考研] 计算机11408，286分求调剂 +7 木子念晞 2026-04-05 7/350 2026-04-05 19:02 by chy09050039
[考研] 一志愿江南大学085501机械工程专硕326分，本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 材料与化工306分找调剂 +12 沧海轻舟e 2026-04-03 13/650 2026-04-04 23:45 by lqwchd
[考研] 349求调剂 +11 zwjjjjjj 2026-03-31 11/550 2026-04-04 19:52 by 蓝云思雨
[考研] 319求调剂 +4 星星不眨眼喽 2026-04-03 4/200 2026-04-04 16:25 by 中飞院空管学院�
[考研] 化工调剂303分，过四级 +28 栖梧待风 2026-04-02 28/1400 2026-04-03 21:40 by qlm5820
[考研] 285求调剂 +5 AZMK 2026-04-03 8/400 2026-04-03 18:17 by AZMK
[考研] 求调剂 +9 akdhjs 2026-03-31 11/550 2026-04-03 13:32 by akdhjs
[考研] 348求调剂 +11 zzzzyk123 2026-04-01 11/550 2026-04-02 16:52 by Wang200018
[考研] 材料化工340求调剂 +5 jhx777 2026-03-30 5/250 2026-04-02 12:45 by smileboy2006
[考研] 322求调剂 +5 熹僖XX 2026-03-31 6/300 2026-04-02 10:08 by 求调剂zz
[考研] 070300化学279求调剂 +15 哈哈哈^_^ 2026-03-31 17/850 2026-04-01 21:37 by 给你你注意休息
[考研] 307分求调剂 +14 (o~o) 2026-03-31 15/750 2026-04-01 20:43 by longlotian
[考研] 材料与化工（0856）304求B区调剂 +8 邱gl 2026-03-30 16/800 2026-04-01 17:58 by 邱gl
[考研] 本科211总分289，08工学真心求调剂 +3 utopiaE 2026-03-30 3/150 2026-03-30 23:42 by ms629
信息提示
关闭
请填处理意见
关闭