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

СОММОN/СРСI/СРС1(20,5),СРС2(20,3)/IDGFD/TOID(8),МСОD(8),МСРD(8)
*/IDGF/CPC(20,8),TOI(8),MCO(8),MCP(8)
COMMON/P/P/Т/Т/Z/Z/TS/RO,PA,W/ETA/ETA
RM=8.31451D0
IF(ICALC.NE.1) GO TO 1
CALL COMPON
IF(Z.EQ.0D0)GO TO 133
DO 11111 J=l,8
DO 11111 I=1,20
IF(J.LE.5) CPC(I,J)=CPC1(I,J)
IF(J.GT.5) CPC(I,J)=CPC2(I,J-5)
11111 CONTINUE
CALL DDIJ(DIJ,LIJ)
DO 75 1=1,NC
TC(I)=TCD(NI(I))
VC(I)=BM(I)/VCD(NI(I))
PII(I)=PIID(NI(I))
MCO(1)=MCOD(NI(I))
MCP(I)=MCPD(NI(I))
TOI(I)=TOID(NI(I))
MP=MCO(I)+MCP(I)+1
DO 23 J=1,MP
23 CPC(J,I)=CPC(J,NI(I))
DO 123 J=1,NC
IF(I.GE.J) GO TO 123
DIJ(I,J)=DIJ(NI(I),NI(J))
LIJ(I,J)=LIJ(NI(I),NI(J))
123 CONTINUE
75 CONTINUE
CALL PARMIX(DIJ,LIJ,TC,VC,P11)
DO 27 1=1,10
DO 27 J=l,8
27 B(I,J)=AIJ(I,J)+BIJ(I,J)*PIM
IF(RON.NE.0D0) THEN
CALL PHASE
RON=0D0
GO TO 133
ENDIF
1 CALL PHASE
133 RETURN
END
SUBROUTINE COMPON
IMPLICIT REAL*8(A-H,O-Z)
DIMENSION BMI(25),ROI(8),GI(8),YI(25)
COMMON/Y/Y(8)/BMM/BMM/BM/BM(8)/YI/YC(25)/NI/NI(8)/NC/NC/RON/RON
DATA BMI/16.043D0,30.07D0,44.097D0,2*58.123D0,28.0135D0,
*44.01D0,34.082D0,26.038D0,28.054D0,42.081D0,3*72.15D0,
*86.177D0,78.114D0,100.204D0,92.141D0,114.231D0,128.259D0,
*142.286D0,4.0026D0,2.0159D0,28.01D0,31.9988D0/
DATA ROI/0.6682D0,1.2601D0,1.8641D0,2.4956D0,2.488D0,