*-H)*P
H3=(RN*(Z1+(12.*(T**11.))))/((Z1+(T**11.))**2.)
H3=H3*(Q1+(((R1/2.)+(P/3.))*P))*P
XH=H0+(RI1*(H1+H2-H3))
TC=(T*T1)-273.16
CALL SATLN(TC,PS)
PS=PS/P1
DV=.02414+(10.**(A2/(T-A3)))*(1.+((P-PS)*A4*(T-A5)))
T=T*T1
P=P*P1
DV=DV*1.E-3
T=T-273.16
RETURN
END
Subroutine SATLN(T,P)
SUBROUTINE SATLN(T,P)
* THIS SUBROUTINE CALCULATES THE SATURATION PRESSURE FOR A
* GIVEN WATER TEMPERATURE. THE WATER TEMPERATURE IS IN C
AND
* THE SATURATION PRESSURE IS RETURNED IN BARS.
DATA RK,A,B,C,D,E,F/2.937E5,5.426651,-2005.1,1.3869E-4,1.1965E-11,
*-.0044,-.0057148/
T=T+273.16
T2=T+.01
Y=647.26-T
X=(T2*T2)-RK
T=T/647.3
A1=A+(B/T2)+(C*X/T2*(10.**(D*X*X)-1.))+(E*(10.**(F*(Y* *1.25))))
P=(1.01325*(10.**A1))+((T-.422)*(.577-T)*EXP(-12.+(T**4.))*
*9.80665E-3)
T=T*647.3
T=T-273.16
RETURN
END
Subroutine ROMBRG(FUN,A,B,C,ERR,RES)
SUBROUTINE ROMBRG(FUN,A,B,C,ERR,RES)
C THIS SUBROUTINE COMPUTES INTEGRALS OF A USER SUPPLIED
C FUNCTION USING ROMBERG'S METHOD. THIS SUBROUTINE IS FROM
C "NUMERICAL METHODS FOR ENGINEERING APPLICATION",J.H.
C FERZIGER, JOHN WILEY AND SONS, 1981. THE ARGUMENT ARE:
C FUN = THE FUNCTION TO BE INTEGRATED
C A = LOWER LIMIT
C B = UPPER LIMIT
C C = ARRAY OF FUNCTION DEFINITION PARAMETERS IF REQUIRED
C ERR = THE DESIRED ACCURACY
C RES = THE RESULTING VALUE FOR THE INTEGRAL
86
Previous Page
