fortran调用refprop焓熵不对

w1_d2

我用fortran调用refprop，用的subroutine TPFLSH (t,p,z,D,Dl,Dv,x,y,q,e,h,s,cv,cp,w,ierr,herr)，计算后得到的焓除以摩尔质量，得到kJ/kg，但与安装的nist软件测得值对比。单工质相同，多组分工质就有大的差异比如差17kJ/kg。求帮助。

w1_d2

程序如下，请帮助下。

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	program ex

02	implicit none

03	integer,parameter          ::ncomp=2

04	character(len=255)         ::hfiles(ncomp),hfmix,herr

05	character(len=3)           ::hrf

06	integer                    ::ierr

07	real(kind=8)               ::t,p,D,Dl,Dv,q,e,h,s,cv,cp,w

08	real(kind=8)               ::z(ncomp),x(ncomp),y(ncomp)

09	real(kind=8)               ::wm

10	real(kind=8),external      ::WMOL

11

12	hfiles(1)="fluids\METHANE.FLD"

13	hfiles(2)="fluids\ETHANE.FLD"

14	hfmix="fluids\hmx.bnc"

15

hrf="DEF"

16

17

t=352.0

18

p=15000.0

19	z=(/0.93,0.07/)

20

21	call SETUP(ncomp,hfiles,hfmix,hrf,ierr,herr)

22	if(ierr .ne. 0)then

23	write(*,*)herr

24

end if

25

26	call TPFLSH(t,p,z,D,Dl,Dv,x,y,q,e,h,s,cv,cp,w,ierr,herr)

27	if(ierr .ne. 0)then

28	write(*,*)herr

29

end if

30	wm=WMOL(z)

31

32

33	write(*,*)"比焓:",h,"J/mol"

34	write(*,*)"比熵:",s,"J/(mol.K)"

35

36	write(*,*)"摩尔质量:",wm,'g/mol'

37

38	end program

​