28
A=A0-A01**0.5
A2=DABS(A)**(1D0/3D0)
IF(A-LT.0D0) A2=-A2
Z=(1D0+A2+A1/A2)/3.
END
BLOCK DATA BDGRG2
IMPLICIT REAL*8(A-H,O-Z)
COMMON/BBLOK/BR11H0(3),BR11H1(3),BR11H2(3),BR22(3),BR23(3),
*BR33(3)/CBLOK/CR111H0(3),CR111H1(3),CR111H2(3),CR222(3),
*CR223(3),CR233(3),CR333(3)
COMMON/ZETA/Z12,Z13,Y12,Y13,Y123
COMMON/MBLOK/GM2,GM3,FA,FB,TO,R
DATA BR11H0/-.425468D0,.2865D-2,-.462073D-5/,
* BR11H1/.877118D-3,-.556281D-5,.881514D-8/,
* BR11H2/-.824747D-6,.431436D-8,-.608319D-11/,
* BR22/-.1446D0,.74091D-3,-.91195D-6/,
* BR23/-.339693D0,.161176D-2,-.204429D-5/,
* BR33/-.86834D0,.40376D-2,-.51657D-5/
DATA CR111H0/-.302488D0,.195861D-2,-.316302D-5/,
* CR111 H1/.646422D-3,-.422876D-5,.688157D-8/,
* CR111H2/-.332805D-6,.22316D-8,-.367713D-11/,
* CR222/.78498D-2,-.39895D-4,.61187D-7/,
* CR223/.552066D-2,-.168609D-4,.157169D-7/,
* CR233/.358783D-2,.806674D-5,-.325798D-7/,
* CR333/.20513D-2,.34888D-4,-.83703D-7/
DATA Z12/.72D0/,Z13/-.865D0/,Y12/.92D0/,Y13/.92D0/,Y123/1.1D0/
DATA GM2/28.0135D0/,GM3/44.01D0/,
* FA/22.414097D0/,FB/22.710811D0/,
* TO/273.15D0/,R/.0831451D0/
END 46
C **********************************************************
С * *
С * Подпрограмма расчета коэффициента сжимаемости природного *
С * газа по уравнению состояния AGA8-92DC. *
C * *
C **********************************************************
SUBROUTINE AGA8DC(ICALC)
IMPLICIT REAL*8(A-H,O-Z)
REAL*8 KI,KIJ,KD
COMMON/RM/RM/Y1/Y(19)/NC1/NC/NI1/NI(19)/EFI/EI(19),KI(19),
*GI(19),QI(19),FI(19)
*/INTER1/EIJ(19,19),UIJ(19,19),KIJ(19,19),GIJ(19,19)
*/EFD/ED(19),KD(19),GD(19),QD(19),FD(19)/Z/Z
RM=8.31448D0
IF(ICALC.NE.l) GO TO 3
CALL COMPO1
IF(Z.EQ.0D0) GO TO 133
CALL PARIN1
DO 75 I=1,NC
EI(I)=ED(NI(I))
KI(I)=KD(NI(I))
GI(I)=GD(NI(I))