24
IF(J.EQ.NP.AND.I+NL-1.EQ.NT) GO TO 29
WRITE(1,7)
WRITE(*,7)
IF(NYST.EQ.l) PAUSE
*’ Для продолжения вывода вставьте бумагу и нажмите <ВВОД> ’
NW=0
IF(NL.GT.1) WRITE(1,17)LIN2(1),(LIN1(K),K=1,NL-1)
IF(NL.EQ.l) WRITE(1,17)LIN2(1)
WRITE(1,19)AT(NL)
IF(NL.GT.1) WRITE(1,21)LIN4,(LIN2(K),K=1,NL-1)
IF(NL.EQ.l) WRITE(1,21)LIN4
WRITE(1,23)(TI(K),K=I,I+NL-1)
WRITE(1,17)(LIN3(K),K=1,NL)
NW=NW+6
ENDIF
25 CONTINUE
15 CONTINUE
29 CLOSE(l)
WRITE(*,7)
PAUSE ’ Вывод завершен, для продолжения работы нажмите <ВВОД> ’
WRITE(*,66)
66 FORMAT(/’ Назначить другое устройство вывода ?’,
*’, 0 - нет, 1 - да ’\)
READ(*,*)NBOLB
IF(NBOLB.EQ.l) GO TO 22
RETURN
END
C **********************************************************
С * *
С * Подпрограмма расчета коэффициента сжимаемости природного *
С * газа по модифицированному методу NX19. *
C **********************************************************
SUBROUTINE NX19(YA,YY)
IMPLICIT REAL*8(A-H,O-Z)
COMMON/NCONT/NCONT/YA/Y(2)/RON/RON
Y(1)=YA
Y(2)=YY
CALL PTCONT
IF(NCONT.EQ.l) GO TO 134
CALL EA
CALL PHASEA
134 RETURN
END
SUBROUTINE PTCONT
IMPLICIT REAL*8(A-H,O-Z)
COMMON/NCONT/NCONT/Z/Z/P/P/T/T/YA/Y(2)/RON/RON
NCONT=0
IF(RON.LT.0.66D0.OR.RON.GT.1D0) NCONT=1
IF(Y(1).GT.0.2D0.OR.Y(2).GT.0.15D0) NCONT=l
IF(P.LE.0.D0.OR.T.LE.0.D0) NCONT=1
IF(T.LT.250.D0.OR.T.GT.340.D0) NCONT=1
IF(P.GT.12.D0) NCONT=1