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

*1.1649D0,1.8393D0,1.4311 D0/
DO 100 I=1,25
100 YI(I)=YC(I)
IF(RON.NE.0D0) GO TO 333
BMM=0D0
DO 3333 I=1,25
3333 BMM=BMM+YI(I)*BMI(I)
333 YS=0D0
DO 55 I=9,25
55 YS=YS+YI(I)
YS1=0D0
DO 67 I=12,21
67 YS1=YS1+YI(I)
YS2=0D0
DO 69 I=22,25
69 YS2=YS2+YI(I)
YI(2)=YI(2)+YI(9)+YI(10)
YI(3)=YI(3)+YI(11)
YI(4)=YI(4)+YS1
YS3=YI(4)+YI(5)
IF(RON.NE.0D0.AND.YI(5).LT.0.01D0.AND.YS3.LT.0.03D0)THEN
YI(4)=YS3
YI(5)=0D0
ENDIF
IF(RON.EQ.0D0.AND.YI(5).LT.0.01D0.AND.YS3.LE.0.03D0)THEN
YI(4)=YS3
YI(5)=0D0
ENDIF
YI(6)=YI(6)+YS2
IF(RON.EQ.0D0) GO TO 555
ROM=0D0
DO 7 I=1,8
7 ROM=ROM+YI(I)*ROI(I)
DO 9 I=1,8
9 GI(I)=YI(I)*ROI(I)/ROM
SUM=0D0
DO 11 I=1,8
11 SUM=SUM+GI(I)/BMI(I)
SUM=1./SUM
DO 13 I=1,8
13 YI(I)=GI(I)*SUM/BMI(I)
555 NC=0
YSUM=0D0
DO 155 I=1,8
IF(YI(I).EQ.0D0) GO TO 155
NC=NC+1
NI(NC)=I
Y(NС)=YI(I)
YSUM=YSUM+Y(NC)
BM(NC)=BMI(I)