28
END
С Подпрограмма, реализующая схему Кардано для определения
С фактора сжимаемости из уравнения состояния
SUBROUTINE ITER2(P,T,Bm,Cm,Z)
IMPLICIT REAL*8(A-H,O-Z)
B1=1D3*P/2.7715/T
B0=Bl*Bm
C0=Bl**2*Cm
A1=1D0+B0
A0=1D0+1.5*(B0+C0)
A01=A0**2-A1**3
IF(A01.LE.0D0) THEN
Z=0D0
RETURN
ENDIF
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)