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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】有很多错误的程序，找不出来了
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1894  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zyj8119

木虫 (著名写手)

[交流] 【求助】有很多错误的程序，找不出来了 已有2人参与

c     此程序是运用peng-robinson方程计算气体的逸度系数faita(1),faita(2)，是利用牛顿迭代法求解
c     然后通过origin拟合得到这些数据N1,N2,m1,m2,K1,K2
c     求解得到f1,f2,x1的过程,首先利用牛顿迭代法求解x1, 然后利用表达式求解f1,f2
      PROGRAM MAIN
        LOGICAL converged
        REAL EXTERNAL FUNCTION MS MSD  G dG
        integer::S N SQ
        parameter(EPS=1.0E-06,R=83.14,epsilon=1-2**(0.5),sigma=1+2**(0.5)
     *)
        dimension:: faita(7),ns(4),ma(4),K(4),A(1000000),M(7),omiga(7),
     *Tc(7),Pc(7),Z(1000000),b(7),beta(7),alpha(7),q(7),as(7),Z1(7)
      real:: alpha,beta,Tr,omiga,M,rho(7),q,a,b,SV,I(7)
        DATA  (ns=(/2.0,3.0,4.0,5.0/);ma=(/1.0,2.0,3.0,
     *4.0/);K=(/5.0,6.0,7.0,8.0/);M=(/16.043,30.070,44.097,
     *58.123,72.150, 44.01, 28.01/);omiga=(/0.012,0.100,0.152,0.2,
     *0.252,0.224,0.048/);Tc=(/190.6,305.3,369.8,425.1,469.7,304.2,
     *132.9/);Pc=(/45.99,48.72,42.48,37.96,33.70,73.83,34.99/)
        READ(*,*)P1
        P=P1/100
      SV=(1-(T/Tc(N))**(0.5))**2
      alpha(N)=(1+(0.37464+1.54226*omiga(N)-0.26992*omiga(N)**2)*SV)
      as(N)=((0.45724*R**2*Tc(N)**2)/Pc(N))*alpha(N)
      b(N)=0.07779*R*Tc(N)/Pc(N)
      beta(N)=b(N)*P/(R*T)
      q(N)=as(N)*alpha(N)*P/(R**2*T**2)
      MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
     *(q*beta(N)-beta(N)**2-beta(N)**3)
        MSD(Z)=3*Z**2-2*(1-beta(N))*Z+(q(N)-2*beta(N)-3*beta(N)**2)
      IF(N.LT.7.AND.N.GT.1)THEN
        DO WHILE(.not.converged.and.ABS(MSD(Z)).GT.0)
        Z=Z-MS(Z)/MSD(Z)
        converged=(abs(MS(Z)/MSD(Z)).LT.EPS)
      END DO
        ELSE
        END IF
        Z1(N)=Z
     
      DO 4 N=1,7
        I(N)=(1/(sigma-epsilon))*log((Z1(N)+sigma*beta(N))/
     *(Z1(N)+epsilon*beta(N)))
        rho(N)=(P/(Z1(N)*R*T))*M(N)*EPS
      faita(N)=exp(Z1(N)-1-log(Z1(N)-beta)-q(N)*I(N))
4    CONTINUE     
   

        f1=P*y1*faita(1)/X
      f2=P*(1-y1)*faita(2)/(1-X)
        G(X)=ns(1)*LOG(1+K(1)*f1**m(1))+ns(2)*LOG(1+K(2)*f1**ma(2))-
     *ns(3)*LOG(1+K(3)*f2**m(3))+ns(4)*LOG(1+K(4)*f2**ma(4))            
      dG(X)=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+K(1)*f1**ma(1))+
     *(ma(2)*ns(2)*K2*f1*(ma(2)-1))/(1+K(2)*f1**m(2))-
     *(ma(3)*ns(3)*K(3)*f2*(ma(3)-1))/(1+K(3)*f2**ma(3))+
     *(ma(4)*ns(4)*K4*f2*(ma(4)-1))/(1+K(4)*f2**ma(4))
      S=0
      DO WHILE(.not.converged.and.ABS(dG(X)).GT.0)
      X=X-G(X)/dG(X)
      S=S+1
      converged=(abs(G(X)/dG(X)).LT.EPS)
      END DO
        x1=X
20    format(1x,3F9.4)
      write(20,*)f1,f2,x1
        end
      --------------------Configuration: bb - Win32 Debug--------------------
Compiling Fortran...
D:\vc++\MSDev98\MyProjects\22\bb.for
D:\vc++\MSDev98\MyProjects\22\bb.for(14) : Error: Syntax error, found ';' when expecting one of: , )
DATA  (ns=(/2.0,3.0,4.0,5.0/);ma=(/1.0,2.0,3.0,
-------------------------------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(27) : Error: This name has not been declared as an array.   [MS]
      MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
------^
D:\vc++\MSDev98\MyProjects\22\bb.for(27) : Error: The shapes of the array expressions do not conform.
      MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
-----------------------------------------------------------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(27) : Error: The shapes of the array expressions do not conform.   [MS]
      MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
------^
D:\vc++\MSDev98\MyProjects\22\bb.for(29) : Error: This name has not been declared as an array.   [MSD]
MSD(Z)=3*Z**2-2*(1-beta(N))*Z+(q(N)-2*beta(N)-3*beta(N)**2)
--------^
D:\vc++\MSDev98\MyProjects\22\bb.for(29) : Error: The shapes of the array expressions do not conform.   [MSD]
MSD(Z)=3*Z**2-2*(1-beta(N))*Z+(q(N)-2*beta(N)-3*beta(N)**2)
--------^
D:\vc++\MSDev98\MyProjects\22\bb.for(37) : Error: The shapes of the array expressions do not conform.   [Z1]
Z1(N)=Z
--------^
D:\vc++\MSDev98\MyProjects\22\bb.for(43) : Error: The shapes of the array expressions do not conform.   [FAITA]
      faita(N)=exp(Z1(N)-1-log(Z1(N)-beta)-q(N)*I(N))
------^
D:\vc++\MSDev98\MyProjects\22\bb.for(49) : Error: This name has not been declared as an array.   [G]
G(X)=ns(1)*LOG(1+K(1)*f1**m(1))+ns(2)*LOG(1+K(2)*f1**ma(2))-
--------^
D:\vc++\MSDev98\MyProjects\22\bb.for(51) : Error: This name has not been declared as an array.   [DG]
      dG(X)=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+K(1)*f1**ma(1))+
------^
Error executing df.exe.

bb.obj - 10 error(s), 0 warning(s)
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
好好学习，天天向上。
1楼 2010-07-11 11:17:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★
余泽成(金币+1):谢谢参与应助！ 2010-07-11 16:27:52
在贴代码时，最好用 Discuz! 代码那块儿显示的 #。这样别人能够很方便地看或拷贝了……

另外，要自己为所有的变量赋初值，比如 MS(ZM) 中的 ZM 就没有赋初值。另外，MS(ZM) 中用到 beta(N)，这里的 N 又是哪一个值呢？

[ Last edited by snoopyzhao on 2010-7-11 at 16:21 ]
赞一下(1人) 回复此楼
7楼2010-07-11 16:15:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

zyj8119

木虫 (著名写手)
改成这样，还有7个错误   
c     此程序是运用peng-robinson方程计算气体的逸度系数faita(1),faita(2)，是利用牛顿迭代法求解
c     然后通过origin拟合得到这些数据N1,N2,m1,m2,K1,K2
c     求解得到f1,f2,x1的过程,首先利用牛顿迭代法求解x1, 然后利用表达式求解f1,f2
      PROGRAM MAIN
        LOGICAL converged
        REAL EXTERNAL FUNCTION MS MSD  G dG
        integer::S N SQ
        parameter(EPS=1.0E-06,R=83.14,epsilon=1-2**(0.5),sigma=1+2**(0.5)
     *)
        dimension:: faita(7),ns(4),ma(4),K(4),M(7),omiga(7),
     *Tc(7),Pc(7),b(7),beta(7),alpha(7),q(7),as(7),Z1(7)
      real:: alpha,beta,Tr,omiga,M,rho(7),q,a,b,SV,I(7)
        DATA  (ns=(/2.0,3.0,4.0,5.0/);ma=(/1.0,2.0,3.0,4.0/);
     *K=(/5.0,6.0,7.0,8.0/);M=(/16.043,30.070,44.097,
     *58.123,72.150, 44.01, 28.01/);omiga=(/0.012,0.100,0.152,0.2,
     *0.252,0.224,0.048/);Tc=(/190.6,305.3,369.8,425.1,469.7,304.2,
     *132.9/);Pc=(/45.99,48.72,42.48,37.96,33.70,73.83,34.99/)
        READ(*,*)P1
        P=P1/100
      SV=(1-(T/Tc(N))**(0.5))**2
      alpha(N)=(1+(0.37464+1.54226*omiga(N)-0.26992*omiga(N)**2)*SV)
      as(N)=((0.45724*R**2*Tc(N)**2)/Pc(N))*alpha(N)
      b(N)=0.07779*R*Tc(N)/Pc(N)
      beta(N)=b(N)*P/(R*T)
      q(N)=as(N)*alpha(N)*P/(R**2*T**2)
      MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
     *(q*beta(N)-beta(N)**2-beta(N)**3)
        MSD(Z)=3*Z**2-2*(1-beta(N))*Z+(q(N)-2*beta(N)-3*beta(N)**2)
      IF(N.LT.7.AND.N.GT.1)THEN
        DO WHILE(.not.converged.and.ABS(MSD(Z)).GT.0)
        Z=Z-MS(Z)/MSD(Z)
        converged=(abs(MS(Z)/MSD(Z)).LT.EPS)
      END DO
        ELSE
        END IF
        Z1(N)=Z
     
      DO 4 N=1,7
        I(N)=(1/(sigma-epsilon))*log((Z1(N)+sigma*beta(N))/
     *(Z1(N)+epsilon*beta(N)))
        rho(N)=(P/(Z1(N)*R*T))*M(N)*EPS
      faita(N)=exp(Z1(N)-1-log(Z1(N)-beta)-q(N)*I(N))
4    CONTINUE     
   

        f1=P*y1*faita(1)/X
      f2=P*(1-y1)*faita(2)/(1-X)
        G(X)=ns(1)*LOG(1+K(1)*f1**m(1))+ns(2)*LOG(1+K(2)*f1**ma(2))-
     *ns(3)*LOG(1+K(3)*f2**m(3))+ns(4)*LOG(1+K(4)*f2**ma(4))            
      dG(X)=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+K(1)*f1**ma(1))+
     *(ma(2)*ns(2)*K2*f1*(ma(2)-1))/(1+K(2)*f1**m(2))-
     *(ma(3)*ns(3)*K(3)*f2*(ma(3)-1))/(1+K(3)*f2**ma(3))+
     *(ma(4)*ns(4)*K4*f2*(ma(4)-1))/(1+K(4)*f2**ma(4))
      S=0
      DO WHILE(.not.converged.and.ABS(dG(X)).GT.0)
      X=X-G(X)/dG(X)
      S=S+1
      converged=(abs(G(X)/dG(X)).LT.EPS)
      END DO
        x1=X
20    format(1x,3F9.4)
      write(20,*)f1,f2,x1
        end
--------------------Configuration: bb - Win32 Debug--------------------
Compiling Fortran...
D:\vc++\MSDev98\MyProjects\22\bb.for
D:\vc++\MSDev98\MyProjects\22\bb.for(14) : Error: Syntax error, found ';' when expecting one of: , )
DATA  (ns=(/2.0,3.0,4.0,5.0/);ma=(/1.0,2.0,3.0,4.0/);
-------------------------------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(27) : Error: This name has not been declared as an array.   [MS]
      MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
------^
D:\vc++\MSDev98\MyProjects\22\bb.for(27) : Error: The shapes of the array expressions do not conform.   [MS]
      MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
------^
D:\vc++\MSDev98\MyProjects\22\bb.for(29) : Error: This name has not been declared as an array.   [MSD]
MSD(Z)=3*Z**2-2*(1-beta(N))*Z+(q(N)-2*beta(N)-3*beta(N)**2)
--------^
D:\vc++\MSDev98\MyProjects\22\bb.for(43) : Error: The shapes of the array expressions do not conform.   [FAITA]
      faita(N)=exp(Z1(N)-1-log(Z1(N)-beta)-q(N)*I(N))
------^
D:\vc++\MSDev98\MyProjects\22\bb.for(49) : Error: This name has not been declared as an array.   [G]
G(X)=ns(1)*LOG(1+K(1)*f1**m(1))+ns(2)*LOG(1+K(2)*f1**ma(2))-
--------^
D:\vc++\MSDev98\MyProjects\22\bb.for(51) : Error: This name has not been declared as an array.   [DG]
      dG(X)=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+K(1)*f1**ma(1))+
------^
Error executing df.exe.
赞一下 回复此楼
好好学习，天天向上。
2楼2010-07-11 11:20:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★ ★ ★
余泽成(金币+3):辛苦了！ 2010-07-11 16:27:27
zyj8119(金币+5):DO循环不是给了范围了吗？那么多公式？ 2010-07-11 17:41:32
看不懂你要做什么，不过有几点建议：

你的程序似乎是 Fortran77 的固定格式，而且里面的循环什么的也是按 Fortran 77 样式写的。但里面有些写法似乎又是 Fortran 90 的样式，这样很不好，建议按一种风格来写你的程序，要么就是纯 Fortran 77 的，要么就是 Fortran90 的。

你的程序中有很多变量没有初值，比如
SV=(1-(T/Tc(N))**(0.5))**2
中的 N；
再如：
MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-(q*beta(N)-beta(N)**2-beta(N)**3)
中的 Z。如果你这里的 MS 是一个函数（就像你前面申明的那样），那么最好用别的符号来代替 Z，否则会引起歧义。

这个程序中最大的问题也就在这些个函数的定义上面了……
赞一下(1人) 回复此楼
3楼2010-07-11 13:45:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)
引用回帖:
Originally posted by snoopyzhao at 2010-07-11 13:45:13:
看不懂你要做什么，不过有几点建议：

你的程序似乎是 Fortran77 的固定格式，而且里面的循环什么的也是按 Fortran 77 样式写的。但里面有些写法似乎又是 Fortran 90 的样式，这样很不好，建议按一种风格来写你 ...

我想定义一个函数，但是错误里面好像非要我搞成数组啊，MS(Z)是一个Z的函数，MSD(Z)是MS(Z)的导数，这样可以利用牛顿迭代法求解，下面的也是，G(X)是函数本身，  dG(X)是G(X)的导数。
赞一下 回复此楼
好好学习，天天向上。
4楼2010-07-11 14:30:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工272求调剂 +16 阿斯蒂芬2004 2026-03-28 16/800 2026-03-29 22:57 by tantai88888
[考研] 085601材料工程找调剂 +14 oatmealR 2026-03-29 15/750 2026-03-29 21:48 by plmuchong
[考研] 一志愿北京工业大学，324分求调剂 +6 零八# 2026-03-28 6/300 2026-03-29 21:20 by nanaliuyun
[考研] 一志愿华东师范大学有机化学专业，初试351分，复试被刷求调剂! +5 真名有冰 2026-03-29 6/300 2026-03-29 20:53 by 唐沐儿
[考研] 一志愿北京理工大学本科211材料工程294求调剂 +8 mikasa的围巾 2026-03-28 8/400 2026-03-29 12:48 by 无际的草原
[考研] 318一志愿吉林大学生物与医药 求调剂 +4 笃行致远. 2026-03-28 4/200 2026-03-29 11:55 by longlotian
[考研] 本科新能源科学与工程，一志愿华理能动285求调剂 +3 AZMK 2026-03-27 5/250 2026-03-28 16:19 by xxxsssccc
[考研] 081200-314 +3 LILIQQ 2026-03-27 4/200 2026-03-28 09:41 by 保护地球你我做�
[考研] 材料求调剂一志愿哈工大324 +7 闫旭东 2026-03-28 9/450 2026-03-28 08:51 by Xu de nuo
[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +9 丹青奶盖 2026-03-26 10/500 2026-03-28 07:45 by barnett0632
[考研] 328求调剂 +7 嗯滴的基本都 2026-03-27 7/350 2026-03-28 04:19 by fmesaito
[考研] 315分求调剂 +7 26考研上岸版26 2026-03-26 7/350 2026-03-28 04:05 by fmesaito
[考研] 07化学280分求调剂 +10 722865 2026-03-23 10/500 2026-03-27 15:51 by Plutoqq
[考研] 305求调剂 +5 哇卢卡库 2026-03-26 5/250 2026-03-27 14:01 by laoshidan
[考研] 343求调剂 +4 赠我一本书 2026-03-23 4/200 2026-03-27 00:40 by wxiongid
[考研] 一志愿天津大学339材料与化工求调剂 +3 江往卖鱼 2026-03-26 3/150 2026-03-26 09:42 by 王小欠i
[考研] 303求调剂 +6 蓝山月 2026-03-25 6/300 2026-03-25 22:47 by 418490947
[考研] 334分 一志愿武理-080500 材料求调剂 +4 李李不服输 2026-03-25 4/200 2026-03-25 21:26 by 星空星月
[考研] 300分，材料，求调剂，英一数二 +5 超赞的 2026-03-24 5/250 2026-03-24 21:07 by 星空星月
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
信息提示
关闭
请填处理意见
关闭