ESUM=ESUM+EABS
EMAX=MAX(EMAX,EABS)
EMIN=MIN(EMIN,EABS)
5 CONTINUE
6 CONTINUE
7 CONTINUE
EAVG=ESUM/N
PRINT 600,EMIN,EMAX,EAVG
STOP
END
Subroutine WTRTBL(T,P,RHO,XH,DV)
SUBROUTINE WTRTBL(T,P,RHO,XH,DV)
* THIS SUBROUTINE CALCULATES THE THERMODYNAMIC AND
* TRANSPORT PROPERTIES OF WATER GIVEN THE TEMPERATURE AND
* PRESSURE CONDITIONS.
* INPUT VALUES ARE T(C) AND P(BARS). THE OUTPUTS ARE DENSITY
* RHO (KG/M3), THE ENTHALPY XH(KJ/KG), AND THE DYNAMIC
* VISCOSITY DV IN KG/M-SEC. THE EQUATIONS FOR THE VISCOSITY ARE
* NOT VALID FOR TEMPERATURES GREATER THAN 300 C.
T=T+273.16
* CALL SATLN(T,P)
DATA RI1,B0,B1,B2,B3,B4,B5,B6,B7,B8,B9,Ul ,Wl/22129.,-37444.8692,
*466453.368,-2666876.77,9030271.53,- 19769400.2,28949239.9,
*-28309932.7,17808942.6,-6534676.01,1065198.53,.58620689,
*.41666667/
DATA G,H,RK,RL,RM,RN,Fl,Gl ,Hl,RKl,RLl,RMl,RNl,Ql,Rl,Zl/
*.417,1.139706E-4,9.949927E-5,7.241165E-5,.7676621,1.052358E- 11,
*3.7E8,3.122199E8,199985.,1.72,1.362926E16,1.500705,.6537154,
*62.5,13.10268,1.5108E-5/
DATA A2,A3,A4,A5/.3828209486,.2162830218,.1498693949,.4711880117/
P1=221.287
T1=647.3
T=T/T1
P=P/P1
U=F1-(G1*T*T)-(H1*(T**(-6.)))
W=U+SORT((RK1*U*U)+(RL1*(P-(RM1*T))))
V1=(G/(W**(1./3.4)))-H+(RK*T)
V2=((ABS(RN1-T))**2.)*(RL+(((ABS(RN1-T))**8.)*RM))
V3=(RN*(R1+(R1*P)+(P*P)))/(Z1+(T**11.))
RHO=1./(V1+V2-V3)
H0=B0+(B1*T)+(B2*T*T)+(B3*T**3.)+(B4*T**4.)+(B5*T**5.)+
*(B6*T**6.)+(B7*T**7.)+(B8*T**8.)+(B9*T**9.)
Y1=RL1*RM1/2.
Q=2.*G/RL1
V=(-2.*G1*T*T)+(6.*H1/(T**6.))
H1=((U1 * W)-(W1*((3.4*U)-V)))*W
H1=Q/(W**(1./3.4))*(H1+(Y1*T)-(.72*V*U))
H2=(((RN1-T)*((RL*(RN1+T))+(RM*((ABS(RN1-T))* *8.)*(RN1+(9.*T)))))
85