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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】求助 程序我问题 解决送80金币
151/2返回列表12下一页
查看: 994  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国家资助计划 +4 bearbo 2024-07-01 4/200 2024-07-01 23:55 by 秘蓝123
[基金申请] 希望青基是一个正反馈 +9 tureAurora 2024-06-30 9/450 2024-07-01 23:12 by Yuanl?h
[基金申请] 自科会评 +13 newtime_1 2024-07-01 14/700 2024-07-01 22:48 by Pickfoot
[教师之家] 院系领导要整人、要孤立普通教师很容易，这位教授被整期间没人敢为他发声 +9 zju2000 2024-06-30 12/600 2024-07-01 21:55 by 秋晴果实
[基金申请] 2020面地青已有会评专家接到电话了 +11 dream200606 2024-06-28 13/650 2024-07-01 21:20 by liliwang215
[基金申请] 出结果了 +28 wuyue526 2024-06-28 34/1700 2024-07-01 18:47 by 徐凤年之我在火�
[基金申请] 恳请大家帮我看看如何提高 +12 drwuwenhan 2024-06-28 15/750 2024-07-01 18:43 by 汤姆喵喵喵
[访问学者] 澳大利亚科廷大学招收化学、化工、材料背景访问学者和访问学生 +3 cfwghj 2024-06-27 5/250 2024-07-01 15:56 by ylbt1989
[基金申请] 先中青基、再中博后面上的情况也挺多吧？ +12 WKKWKK 2024-06-29 13/650 2024-07-01 12:34 by miner007
[基金申请] 特助中标的话看不到评审排名吗 +4 outsider1986 2024-06-30 7/350 2024-07-01 08:42 by 小小614
[访问学者] no news is good news +8 wyjecho666 2024-06-29 12/600 2024-06-30 13:49 by chang_jiang
[基金申请] 博后面上如何提高创新分？ +5 yuyiang 2024-06-29 10/500 2024-06-29 18:24 by yuyiang
[论文投稿] 我投稿了Elsevier，后来申请撤稿，你不允许我撤稿，强制发表然后问我要钱 5+4 996247712 2024-06-26 12/600 2024-06-29 15:17 by bobvan
[基金申请] 统计一波计分系数 +5 lyfbangong 2024-06-28 6/300 2024-06-29 13:30 by lddscut
[论文投稿] Scitific Reports +5 智晨君 2024-06-27 5/250 2024-06-29 10:37 by zhulei513
[有机交流] 过柱子，产品，杂质在是 pe:ea=100:1 也一起出来？ +6 w256 2024-06-25 6/300 2024-06-28 18:15 by 我想更懂你29
[基金申请] 博后网上办公系统维护 +14 取名字烦人 2024-06-25 23/1150 2024-06-28 17:22 by cincpac
[基金申请] 第17批博士后特别资助名单 +7 grdsluyu 2024-06-27 7/350 2024-06-28 09:58 by 喵呜呜_c
[有机交流] 大佬们，打的核磁氢谱与chemdraw预测的有偏差 +7 来了个 2024-06-26 7/350 2024-06-27 11:11 by wubiansiya
[有机交流] 酚羟基甲基化 10+4 A好运来啦啦啦 2024-06-25 6/300 2024-06-26 10:08 by 88817753
信息提示
关闭
请填处理意见
关闭