11
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 FORMAT(//
*10X,’???? Рассчитаны следующие физические свойства ???’/
*10Х,’ ’/
*10Х,’ 1. Плотность ’/
*10Х,’ ’/
*10Х,’ 2. Показатель адиабаты ’/
*10Х,’ ’/
*10Х,’ 3. Скорость звука ’/
*10Х,’ ’/
*10Х,’ 4. Коэффициент динамической вязкости ’/
*10Х,’ ’/
*10Х,’???????????????????????????’/)
WRITE(*,5)
5 FORMAT(/,3X,
*’Введите порядковый номер свойства для вывода результатов расче’,
*’та’/
*’ или 5 для выхода в ДОС ’\)
READ(*,*)NPROP
RETURN
END
SUBROUTINE RANGE(NRANGE)
IMPLICIT REAL*8(A-H,O-Z)
COMMON/Z/Z
WRITE(*,1)
1 FORMAT(//
*’ Метод расчета при заданных параметрах “не работает”’/
*’ Продолжить работу программы ? 0 - нет, 1 - да ’\)
READ(*,*)NRANGE
RETURN
END