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

NP=NP-1
ELSE
I=I+1
ENDIF
lF(I.LE.NP) GO TO 9
J=l
15 JS=0
DO 17 I=1,NP
IF(ROP(I,J).EQ.0D0) JS=JS+1
17 CONTINUE
IF(JS.EQ.NP) THEN
IF(J.NE.NT) THEN
DO 19 I=J,NT-1
TI(I)=TI(I+1)
DO 19 K=1,NP
ROP(K,I)=ROP(K,I+1)
PAP(K,I)=PAP(K,I+1)
WP(K,I)=WP(K,I+1)
19 ETAP(K,I)=ETAP(K,I+1)
ENDIF
NT=NT-1
ELSE
J=J+1
ENDIF
IF(J.LE.NT) GO TO 15
CALL PROP(NPROP)
IF(NPROP.EQ.5) GO TO 134
IF(NPROP.EQ.l) CALL TABL(PI,TI,ROP,NP,NT,NPROP)
IF(NPROP.EQ.2) CALL TABL(PI,TI,PAP,NP,NT,NPROP)
IF(NPROP.EQ.3) CALL TABL(PI,TI,WP,NP,NT,NPROP)
IF(NPROP.EQ.4) CALL TABL(PI,TI,ETAP,NP,NT,NPROP)
WRITE(*,'(A\)')
*' Продолжить вывод рассчитанных свойств ? 0 - нет, 1 - да '
READ(*,*)NCONT
IF(NCONT.EQ.l) GO TO 500
ENDIF
134 STOP
END
SUBROUTINE PROP(NPROP)
WRITE(*,1)
1 FORMATS(//
*10X,'─────Рассчитаны следующие физические свойства─────────────'/
*10Х,' '/
*10Х,' 1. Плотность '/
*10Х,' '/
*10Х,' 2. Показатель адиабаты '/
*10Х,' '/
*10Х,' 3. Скорость звука '/
*10Х,' '/
*10Х,' 4. Коэффициент динамической вязкости '/