11
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
SUBROUTINE TABL(PI,TI,ZP,NP,NT,NPROP)
IMPLICIT REAL*8(A-H,O-Z)
CHARACTER*26 AR,FNAME
CHARACTER PROP(4)*58,A*6,LIN1(5)*9,LIN2(5)*9,LIN3(6)*9,LIN4*9,
*AT(6)*28,RAZM(4)*39
CHARACTER*70 F,FZ(11,2),FW(11,2)
DIMENSION PI(100),TI(100),ZP(100,100),ZPP(6)