版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

[交流] 【求助】有很多错误的程序，找不出来了已有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)

回复此楼

» 猜你喜欢

一志愿南航 335分 | 0856 | GPA 4.07 | 有科研经历已经有8人回复
求调剂已经有5人回复
求化学调剂已经有5人回复
085600，材料与化工321分求调剂已经有10人回复
0703 化学求调剂，一志愿山东大学 342 分已经有5人回复
一志愿南开大学0710生物学359求调剂已经有3人回复
085600，专业课化工原理，320分求调剂已经有4人回复
材料与化工272求调剂已经有16人回复
290求调剂已经有3人回复
327求调剂已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

Perl 生物信息学编程求助，运行出现问题，烦请帮忙查找一下错误，谢谢！已经有3人回复
求高人帮忙看一下我这个Fortran程序的错误出在哪里？已经有12人回复
出现应用程序错误已经有8人回复
求助！！！！安装MS5.0程序过程中出现的错误提示已经有12人回复
帮我看看这里错误提示是什么意思？已经有3人回复
【求助】此程序没有错误，怎么运行不出来？已经有6人回复
【求助】程序运行后出现的错误 access violation 已经有15人回复

好好学习，天天向上。

1楼 2010-07-11 11:17:43

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

改成这样，还有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

至尊木虫 (职业作家)

程序强帖: 16
应助: 157 (高中生)
贵宾: 0.02
金币: 18844.7
红花: 29
帖子: 3803
在线: 1422.4小时
虫号: 183750
注册: 2006-02-13
专业: 污染生态化学

★ ★ ★
余泽成(金币+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

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

引用回帖:

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的回帖

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

$    此程序是运用peng-robinson方程计算气体的逸度系数faita(1),faita(2)，是利用牛顿迭代法求解
$    然后通过origin拟合得到这些数据N1,N2,m1,m2,K1,K2
$    求解得到f1,f2,x1的过程,首先利用牛顿迭代法求解IAST理论得到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),ZM
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
$    计算peng-robinson方程中需要的一些参量
DO 3 N=1,7
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)
3    CONTINUE
$    运用牛顿迭代法开始求解压缩因子的方程
MS(ZM)=ZM**3-(1-beta(N))*ZM**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
   *(q*beta(N)-beta(N)**2-beta(N)**3)
MSD(ZM)=3*ZM**2-2*(1-beta(N))*ZM+(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(ZM)).GT.0)
ZM=ZM-MS(ZM)/MSD(ZM)
converged=(abs(MS(ZM)/MSD(ZM)).LT.EPS)
   END DO
ELSE
END IF
Z1(N)=ZM
$    求解密度和逸度
   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

$    运用IAST理论以及上面得到的逸度系数，开始求解IAST的方程
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

赞一下

回复此楼

好好学习，天天向上。

5楼2010-07-11 14:53:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

程序强帖: 16
应助: 157 (高中生)
贵宾: 0.02
金币: 18844.7
红花: 29
帖子: 3803
在线: 1422.4小时
虫号: 183750
注册: 2006-02-13
专业: 污染生态化学

★
余泽成(金币+1):谢谢参与应助！ 2010-07-11 16:27:44

下面这两个方程，是不能定义成这个样子的。
      MS(ZM)=ZM**3-(1-beta(N))*ZM**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
   *(q*beta(N)-beta(N)**2-beta(N)**3)
      MSD(ZM)=3*ZM**2-2*(1-beta(N))*ZM+(q(N)-2*beta(N)-3*beta(N)**2)
从你的代码上看，这两个方程中的 beta(N), q(N) 啥的都是变量，呵呵……

你这么搞搞看：

MS(ZM,betaN,qN)=ZM**3-(1-betaN)*ZM**2+(qN-2*betaN-3*betaN**2)*ZM-(qN*betaN-betaN**2-betaN**3)
.....
然后调用的时候，记得给 ZM, betaN 和 qN 赋值。其它的函数应该类似。

赞一下(1人)

回复此楼

6楼2010-07-11 15:46:58

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

程序强帖: 16
应助: 157 (高中生)
贵宾: 0.02
金币: 18844.7
红花: 29
帖子: 3803
在线: 1422.4小时
虫号: 183750
注册: 2006-02-13
专业: 污染生态化学

★
余泽成(金币+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的回帖

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

$    此程序是运用peng-robinson方程计算气体的逸度系数faita(1),faita(2)，是利用牛顿迭代法求解
$    然后通过origin拟合得到这些数据N1,N2,m1,m2,K1,K2
$    求解得到f1,f2,x1的过程,首先利用牛顿迭代法求解IAST理论得到x1, 然后利用表达式求解f1,f2
   PROGRAM MAIN
LOGICAL converged
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),ZM,SSA,SSB,SSC,
   *SSD
DATA  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
$    计算peng-robinson方程中需要的一些参量
DO 3 N=1,7
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)
3    CONTINUE
$    运用牛顿迭代法开始求解压缩因子的方程

   IF(N.LT.7.AND.N.GT.1)THEN
ZM=1
DO WHILE(.not.converged.and.ABS(MSD(ZM)).GT.0)
SSA=MS(ZM,beta(N),q(N))
SSB=MSD(ZM,beta(N),q(N))
ZM=ZM-SSA/SSB
converged=(abs(SSA)/(SSB))).LT.EPS)
   END DO
ELSE
END IF
do 4 N=1,7
Z1(N)=ZM
4    CONTINUE
$    求解密度和逸度
   DO 5 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)=SVF(Z1(N),beta,q(N),I(N))
5 CONTINUE

$    运用IAST理论以及上面得到的逸度系数，开始求解IAST的方程
S=0
X=1
SSC=iast(X,P,y1,faita(1),faita(2))
SSD=diast(X,P,y1,faita(1),faita(2))
   DO WHILE(.not.converged.and.ABS(SSD).GT.0)
   X=X-SSC/SSD
   S=S+1
   converged=(abs(SSC/SSD).LT.EPS)
   END DO
x1=X
20 format(1x,3F9.4)
   write(20,*)f1,f2,x1
end


function MS(SA,SB,SC)
   MS(SA,SB,SC)=SA**3-(1-SB)*SA**2+(SC-2*SB-3*SB
   ***2)*SA-(SC*SB-SB**2-SB**3)
end

function MSD(SD,SE,SF)
   MSD(SD,SE,SF)=3*SD**2-2*(1-SE)*SD+(SF-2*SE-
   *3*SE**2)
end

function iast(MADD,P,YY,SG1,SG2)
   dimension  ns(4),ma(4),K(4)
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/)
f1=P*YY*SG1/MADD
   f2=P*(1-YY)*SG2/(1-MADD)
iast(MADD,P,YY,SG1,SG2)=ns(1)*LOG(1+K(1)*f1**ma(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))
   end

function  diast(MAD,S,Y,SG3,SG4)
   dimension  ns(4),ma(4),K(4)
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/)
f1=S*Y*SG3/MAD
   f2=S*(1-Y)*SG4/(1-MAD)
diast(MAD,S,Y,SG3,SG4)=(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))
END

FUNCTION SVF(SG5,SG6,SG7,SG8)
   SVF(SG5,SG6,SG7,SG8)=exp(SG5-1-log(SG5-SG6)-SG7*SG8)
END
改成这样，还是有9个错误啊。。。

赞一下

回复此楼

好好学习，天天向上。

8楼2010-07-11 17:34:13

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

-----------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(76) : Error: This global name is invalid in this context. [MSD]
   MSD(SD,SE,SF)=3*SD**2-2*(1-SE)*SD+(SF-2*SE-
------^
D:\vc++\MSDev98\MyProjects\22\bb.for(75) : Warning: The return value of this FUNCTION has not been defined. [MSD]
function MSD(SD,SE,SF)
-----------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(83) : 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(87) : Error: This global name is invalid in this context. [IAST]
iast(MADD,P,YY,SG1,SG2)=ns(1)*LOG(1+K(1)*f1**ma(1))+ns(2)*LOG
--------^
D:\vc++\MSDev98\MyProjects\22\bb.for(81) : Warning: The return value of this FUNCTION has not been defined. [IAST]
function  iast(MADD,P,YY,SG1,SG2)
------------------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(95) : 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(99) : Error: This global name is invalid in this context. [DIAST]
diast(MAD,S,Y,SG3,SG4)=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+
--------^
D:\vc++\MSDev98\MyProjects\22\bb.for(93) : Warning: The return value of this FUNCTION has not been defined. [DIAST]
function  diast(MAD,S,Y,SG3,SG4)
------------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(106) : Error: This global name is invalid in this context. [SVF]
   SVF(SG5,SG6,SG7,SG8)=exp(SG5-1-log(SG5-SG6)-SG7*SG8)
------^
D:\vc++\MSDev98\MyProjects\22\bb.for(105) : Warning: The return value of this FUNCTION has not been defined. [SVF]
FUNCTION SVF(SG5,SG6,SG7,SG8)
-----------------^
Error executing df.exe.

bb.obj - 9 error(s), 5 warning(s)

赞一下

回复此楼

好好学习，天天向上。

9楼2010-07-11 17:39:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

beta(N)和q(N)本来就是计算得到的数组数值啊，怎么给初值呢，前面已经算了啊。

赞一下

回复此楼

好好学习，天天向上。

10楼2010-07-11 17:40:51

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 zyj8119 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 南京大学化学调剂 +9	景随风 2026-03-29	14/700	2026-03-29 21:22 by chixmc
[考研] 277跪求调剂 +6	1915668 2026-03-27	10/500	2026-03-29 16:03 by 王亮_大连医科大
[考研] 298求调剂 +3	种圣赐 2026-03-29	3/150	2026-03-29 12:06 by longlotian
[考研] 086000生物与医药调剂 +5	Feisty。 2026-03-28	9/450	2026-03-29 12:02 by longlotian
[考研] 279求调剂 +4	蝶舞轻绕 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 调剂求院校招收 +6	鹤鲸鸽 2026-03-28	6/300	2026-03-29 08:15 by fmesaito
[考研] 一志愿华北电力大学能动专硕，293，求调剂 +3	15537177284 2026-03-23	5/250	2026-03-28 16:11 by xxxsssccc
[考研] 081200-314 +3	LILIQQ 2026-03-27	4/200	2026-03-28 09:41 by 保护地球你我做�
[考研] 085405 考的11408求各位老师带走 +3	Qiu学ing 2026-03-28	3/150	2026-03-28 09:19 by 乐呵呵的追梦人
[考研] 张芳铭-中国农业大学-环境工程专硕-298 +4	手机用户 2026-03-26	4/200	2026-03-28 07:17 by mmm just
[考研] 考研调剂 +4	Sanmu-124 2026-03-26	4/200	2026-03-27 17:49 by kiokin
[考研] 308求调剂 +7	墨墨漠 2026-03-25	7/350	2026-03-27 14:47 by 狂炫麦当当
[考研] 一志愿郑大085600，310分求调剂 +5	李潇可 2026-03-26	5/250	2026-03-27 11:14 by 不吃魚的貓
[论文投稿] Journal of Mechanical Science and Technology +3	Russ_ss 2026-03-25	5/250	2026-03-27 10:49 by 陆小果画大饼
[考研] 336材料求调剂 +7	陈滢莹 2026-03-26	9/450	2026-03-27 00:20 by wxiongid
[考研] 生物学 296 求调剂 +4	朵朵- 2026-03-26	6/300	2026-03-26 19:01 by 不吃魚的貓
[考研] 334分一志愿武理-080500 材料求调剂 +4	李李不服输 2026-03-25	4/200	2026-03-25 21:26 by 星空星月
[考研] 302求调剂 +4	锦衣卫藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫疯狂
[有机交流] 有机合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 292求调剂 +4	鹅鹅鹅额额额额� 2026-03-24	4/200	2026-03-24 16:41 by peike

24小时热门版块排行榜

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

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