19
11 CONTINUE
IF(IC.EQ.NT) THEN
IF(I.NE.NP) THEN
DO 13 J=I,NP-1
PI(J)=PI(J+1)
DO 13 K=1,NT
13 ZP(J,K)=ZP(J+1,K)
ENDIF
NP=NP-1
ELSE
I=I+1
ENDIF
IF(I.LE.NP) GO TO 9
J=l
15 JS=0
DO 17 I=1,NP
IF(ZP(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
ТI(I)=ТI(I+1)
DO 19 K=1,NP
19 ZP(K,I)=ZP(K,I+1)
ENDIF
NT=NT-1
ELSE
J=J+1
ENDIF
IF(J.LE.NT) GO TO 15
CALL TABL(YA,YY,PI,TI,ZP,NP,NT,NVAR,AR)
ENDIF
GO TO 200
134 STOP
END
SUBROUTINE VAR(NVAR)
WRITE(*,1)
1 FORMAT(//
*10X,’ Расчет коэффициента сжимаемости природного газа’//
*10Х,’ ----------------Метод расчета----------------- ’/
*10Х,’ ’/
*10Х,’ 1. Модифицированный метод NX 19 ’/
*10Х,’ ’/
*10Х,’ 2. Уравнение состояния GERG-91 ’/
*10Х,’ ’/
*10Х,’ 3. Уравнение состояния AGA8-92DC ’/
*10Х,’ ’/
*10Х,’ 4. Уравнение состояния ВНИЦ СМВ ’/
*10Х,’ ’/
*10Х,’---------------------------------------------------’/)
WRITE(*,5)
5 FORMAT(/,3X,