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

*' Ввод исходных данных завершен.'
P=.101325D0
T=293.15D0
ICALC=1
CALL EOSVNIC(ICALC)
IF(Z.EQ.0D0) THEN
CALL RANGE(NRANGE)
IF (NRANGE) 134,134,200
ENDIF
ICALC=2
NTS=0
DO 7 I=1,NP
P=PI(I)
DO 7 J=1,NT
T=TI(J)
CALL EOSVNIC(ICALC)
IF(Z.NE.0D0) THEN
NTS=NTS+1
ROP(I,J)=RO
PAP(I,J)=PA
WP(I,J)=W
ETAP(I,J)=ETA
ELSE
ROP(I,J)=0D0
PAP(I,J)=0D0
WP(I,J)=0D0
ETAP(I,J)=0D0
ENDIF
7 CONTINUE
500 WRITE(*,100)
100 FORMAT(25(/))
IF(NTS.EQ.0) THEN
CALL RANGE(NRANGE)
IF (NRANGE) 134,134,200
ELSE
I=1
9 IS=0
DO 11 J=1,NT
IF(ROP(I,J).EQ.0D0) IS=IS+1
11 CONTINUE
IF(IS.EQ.NT) THEN
IF(I.NE.NP) THEN
DO 13 J=I,NP-1
PI(J)=PI(J+1)
DO 13 K=1,NT
ROP(J,K)=ROP(J+1,K)
PAP(J,K)=PAP(J+1,K)
WP(J,K)=WP(J+1,K)
13 ETAP(J,K)=ETAP(J+1 ,K)
ENDIF