26
С * Подпрограмма расчета коэффициента сжимаемости природного *
С * газа по модифицированному уравнению состояния 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
X3=YY
IF(RON.LT.0.66D0.OR.RON.GT.1D0)Z=0D0
IF(X2.LT.0D0.OR.X2.GT.0.2D0)Z=0D0
IF(X3.LT.0D0.OR.X3.GT.0.15D0) Z=0D0
IF(Z.EQ.0D0) GO TO 133
X1=1D0-X2-X3
Х11=Х1*Х1
X12=X1*X2
X13=X1*X3
X22=X2*X2
X23=X2*X3
X33=X3*X3
Z=1D0-(.0741*RON-.006D0-.063*YA-.0575*YY)**2
BMNG=24.05525*Z*RON
Y1=1D0-YA-YY
BMY=(BMNG-28.0135*YA-44.01*YY)/Y1
С Расчет теплоты сгорания эквивалентного углеводорода (Н)
H=47.479*BMY+128.64D0
RETURN
3 Т=Т1
ТС=Т1-Т0
P=PRESS
IF(PRESS.LE.0D0.OR.PRESS.GT.12D0)Z=0D0
IF(T1.LT.250D0.OR.T1.GT.340D0)Z=0D0
IF(Z.EQ.0D0) GO TO 133
CALL B11BER(T,H,B11)
CALL BBER(T,B11,B,Z)
IF(Z.EQ.0D0) GO TO 133
CALL CBER(T,H,C,Z)
IF(Z.EQ.0D0) GO TO 133
CALL ITER2(P,T,B,C,Z)
133 RETURN
END
SUBROUTINE B11BER(T,H,B11)
IMPLICIT REAL*8(A-H,O-Z)
COMMON/BBLOK/BR11H0(3),BR11H1(3),BR11H2(3),BR22(3),BR23(3),BR33(3)
T2=T*T