Межгосударственный стандарт ГОСТ 30319.3-96 "Газ природный. Методы расчета физических свойств. Определение физических свойств по уравнению состояния" (введен в действие постановлением Государственного комитета РФ по стандартизации, метрологии и сертификац стр. 16

155 CONTINUE
CALL MOLDOL(YI,YS)
DO 551 I=1,NC
551 Y(I)=Y(I)/YSUM
RETURN
END
SUBROUTINE MOLDOL(YI,YS)
IMPLICIT REAL*8(A-H,O-Z)
DIMENSION YI(25)
COMMON/Z/Z
Z=-1D0
IF(YI(l).LT.0.5D0.OR.YI(2).GT.0.2D0.OR.YI(3).GT.0.05D0.OR.
*YI(4).GT.0.03D0.OR.YI(5).GT.0.03D0.OR.YS.GT.0.01D0)Z=0D0
IF(YI(6).GT.0.3D0.OR.YI(7).GT.0.3D0.OR.YI(8).GT.0.3D0)Z=0D0
RETURN
END
SUBROUTINE DDIJ(DIJ,LIJ)
IMPLICIT REAL*8(A-H,0-Z)
REAL*8 LIJ(8,8)
DIMENSION DIJ(8,8)
DO 1 I=1,8
DO 1 J=l,8
LIJ(I,J)=0.D0
1 DIJ(I,J)=0.D0
DIJ(l,2)=0.036D0
DIJ(l,3)=0.076D0
DIJ(l,4)=0.121D0
DIJ(l,5)=0.129D0
DIJ(l,6)=0.06D0
DIJ(l,7)=0.074D0
DIJ(2,6)=0.106D0
DIJ(2,7)=0.093D0
DIJ(6,7)=0.022D0
DIJ(l,8)=0.089D0
DIJ(2,8)=0.079D0
DIJ(6,8)=0.211D0
DIJ(7,8)=0.089D0
LIJ(l,2)=-0.074D0
LIJ(l,3)=-0.146D0
LIJ(l,4)=-0.258D0
LIJ(l,5)=-0.222D0
LIJ(l,6)=-0.023D0
LIJ(l,7)=-0.086D0
LIJ(6,7)=-0.064D0
LIJ(7,8)=-0.062D0
RETURN
END
SUBROUTINE PARMIX(DIJ,LIJ,TC,VC,PII)
IMPLICIT REAL*8(A-H,O-Z)
REAL*8 LIJ(8,8)