25
CALL PTCONT
IF(NCONT.EQ.l) GO TO 134
CALL EA
CALL PHASEA
134 RETURN
END
SUBROUTINE PTCONT
IMPLICIT REAL*8(A-H,O-Z)
COMMON/NCONT/NCONT/Z/Z/P/P/T/T/YA/Y(2)/RON/RON
NCONT=0
IF(RON.LT.0.66D0.OR.RON.GT.1D0) NCONT=1
IF(Y(1).GT.0.2D0.OR.Y(2).GT.0.15D0) NCONT=l
IF(P.LE.0.D0.OR.T.LE.0.D0) NCONT=1
IF(T.LT.250.D0.OR.T.GT.340.D0) NCONT=1
IF(P.GT.12.D0) NCONT=1
IF(NCONT.EQ.1) Z=0D0
RETURN
END
SUBROUTINE EA
IMPLICIT REAL*8(A-H,O-Z)
COMMON/T/T/YA/Y(2)/RON/RON/P/P/PT/PA,TA/BI/B1,B2/T0/T0
PCM=2.9585*(1.608D0-0.05994*RON+Y(2)-.392*Y(1))
TCM=88.25*(0.9915D0+1.759*RON-Y(2)-1.681*Y(1))
PA=0.6714*P/PCM+0.0147
TA=0.71892*T/TCM+0.0007
DTA=TA-1.09D0
F=0D0
IF(PA.GE.0D0.AND.PA.LT.2D0.AND.DTA.GE.0D0.AND.DTA.LT.0.3D0)
F=75D-5*PA**2.3/DEXP(20.*DTA)+
*11D-4*DTA**0.5*(PA*(2.17D0-PA+1.4*DTA**0.5))**2
IF(PA.GE.0D0.AND.PA.LT.1.3D0.AND.DTA.GE.-0.25D0.AND.DTA.LT.0D0)
*F=75D-5*PA**2.3*(2D0-DEXP(20.*DTA))+
*1.317*PA*(1.69D0-PA**2)*DTA**4
IF(PA.GE.1.3D0.AND.PA.LT.2D0.AND.DTA.GE.-0.21D0.AND.DTA.LT.0D0)
*F=75D-5*PA**2.3*(2D0-DEXP(20.*DTA))+
*0.455*(1.3D0-PA)*(1.69*2.D0**1.25-PA**2)*(DTA*(0.03249D0+
*18.028*DTA**2)+DTA**2*(2.0167D0+DTA**2*(42.844D0+200.*DTA**2)))
T1=:TA**5/(TA**2*(6.60756*TA-4.42646D0)+3.22706D0)
T0=(TA**2*(1.77218D0-0.8879*TA)+0.305131D0)*T1/TA**4
B1=2.*T1/3.-TO**2
B0=T0*(T1-T0**2)+0.1*T1*PA*(F-1D0)
B2=(B0+(B0**2+B1**3)**0.5)**(1D0/3D0)
RETURN
END
SUBROUTINE PHASEA
IMPLICIT REAL*8(A-H,O-Z)
COMMON/Z/Z/PT/PA,TA/BI/B1,B2/T0/T0
Z=(1D0+0.00132/TA**3.25)**2*0.1*PA/(B1/B2-B2+T0)
RETURN
END
C *************************************************************
С * *