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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】求助 程序我问题 解决送80金币
51/1返回列表
查看: 995  |  回复: 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

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

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的回帖
查看全部 15 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 普通老师、教研室主任、院长副院长都要做很多杂事，学科带头人、专业负责人最舒适 +9 河西夜郎 2024-07-03 11/550 2024-07-04 01:13 by 河西夜郎
[论文投稿] 硕士论文知网可查还能整理成英文论文投稿发表吗 +6 魔芋丝结 2024-07-02 6/300 2024-07-03 21:34 by 投必得科研顾问
[基金申请] 史上最严的基金限项 +10 MathGLA 2024-07-01 12/600 2024-07-03 18:07 by kuaile5420
[基金申请] 孩子，我要告诉你 (金币+5) +26 青涩学子 2024-07-01 29/1450 2024-07-03 17:45 by 超级无敌华子
[基金申请] 科研之友最近访客那出现来自北京的未知用户，什么情况？ +31 lihangshi 2024-06-28 73/3650 2024-07-03 15:04 by 自若远山HH
[基金申请] 国资B档含金量咋样 +8 kkkjgfd 2024-07-02 9/450 2024-07-03 13:17 by 右右0938
[基金申请] 名单中没有，评审结果反馈也没有信息，有没有相同的小伙伴指教 +6 花生微微 2024-06-28 10/500 2024-07-03 09:56 by Nice_bro
[论文投稿] 《化学进展》投稿时3篇代表作 10+3 yicheong 2024-07-01 6/300 2024-07-03 07:31 by kanyechris
[基金申请] 特助没中，请问大家能看到评审意见吗？ +3 Lynn212 2024-07-02 3/150 2024-07-03 06:54 by 勇敢的female
[基金申请] 自科会评 +15 newtime_1 2024-07-01 17/850 2024-07-02 23:17 by kobe0107
[基金申请] 请教大佬，博后面上申请透漏个人信息 +5 Tansund 2024-07-01 7/350 2024-07-02 17:20 by kangshisan
[论文投稿] 投稿JAMT 2+3 sdfuu 2024-06-30 6/300 2024-07-02 13:46 by bobvan
[教师之家] 佛山大学和深职大怎么选？ +11 24601 2024-06-29 22/1100 2024-07-02 12:33 by 24601
[基金申请] 青拔答辩通知出来了 +6 5011fudan 2024-06-28 6/300 2024-07-01 19:01 by 学员qYccD0
[考博] 30岁 想读非全博，求建议 +8 八千里路云雨泥 2024-06-27 22/1100 2024-07-01 09:32 by nanoeye
[基金申请] 特助中标的话看不到评审排名吗 +4 outsider1986 2024-06-30 7/350 2024-07-01 08:42 by 小小614
[基金申请] 中了博后面上系统能查到分数和排名吗？ +5 老黑不黑 2024-06-29 6/300 2024-07-01 06:19 by yqliu130
[教师之家] 唉 +12 小木虫用户100号 2024-06-28 12/600 2024-06-30 23:19 by kanyechris
[访问学者] no news is good news +8 wyjecho666 2024-06-29 12/600 2024-06-30 13:49 by chang_jiang
[基金申请] 今年国自必中，不中回家种田 +11 wutzxt 2024-06-28 12/600 2024-06-30 09:06 by shl2112501
信息提示
关闭
请填处理意见
关闭