25
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 *************************************************************
С * *
С * Подпрограмма расчета коэффициента сжимаемости природного *
С * газа по модифицированному уравнению состояния GERG-91. *
С * *
C *************************************************************
$NOTRUNCATE
SUBROUTINE GERG2(ICALC,YA,YY)
IMPLICIT REAL*8(A-H,O-Z)
COMMON/T/T1/P/PRESS/RON/RON/Z/Z
COMMON/XBLOK/X1,X2,X3,X11,X12,X13,X22,X23,X33
COMMON/MBLOK/GM2,GM3,FA,FB,TO,R
DATABMO/.0838137D0/,BM1/-.00851644D0/,WD0/134.2153D0/,
*WD1/1067.943D0/
Z=-1D0
IF(ICALC.EQ.2) GO TO 3
X2=YA