24
WRITE(*,7)
PAUSE ’ Для продолжения вывода нажмите <ВВОД> ’
WRITE(*,100)
NW=0
WRITE(1,7)
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
IF(NW.EQ.54.AND.NYST.NE.0) THEN
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