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

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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)

回复此楼

» 猜你喜欢

好好学习，天天向上。

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

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

snoopyzhao

至尊木虫 (职业作家)

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

★
resonant(金币+1):指导费：－） 2010-07-11 22:39:13

另外，

CODE:

      SV=(1-(T/Tc(N))**(0.5))**2

这里的 T 未定义。程序中 converged 也未赋初始值。

Warning 还是有一些的，比如一些未使用的变量……

赞一下(1人)

回复此楼

16楼2010-07-11 21:25:41

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

查看全部 19 个回答

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

查看全部 19 个回答

24小时热门版块排行榜

zyj8119

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

» 猜你喜欢

snoopyzhao

zyj8119

snoopyzhao

zyj8119

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