Название: Компьютерное моделирование теплофизических процессов - Методические указания (А.В. Чичиндаев)

Жанр: Информатика

Просмотров: 982


2.3. текст программы (файл r2.for)

C       ТЕМПЕРАТУРА ТОЧКИ РОСЫ R2.FOR

          BLOCK DATA

          COMMON/DAN/D0,DO,N,P,TR

          DATA D0,DO,N,P/0.,20.,100,101.325/

          END

          COMMON/DAN/D0,DO,N,P,TR

          OPEN(50,FILE='R2.REZ')

  11    CALL UID(IO)

          IF(IO.EQ.0) GOTO 10

          WRITE(50,3) P

          DSH=(DO-D0)/N

          DO 1 I=1,N

          D=D0+DSH*I

          CALL TTR(D,P)

  1      CONTINUE

          GOTO 11

  10    CONTINUE

  3      FORMAT(' P=',F8.2)

          STOP

          END

          SUBROUTINE TTR(D,P)

          PHU=P/(622./D+1.)

          IF(T.GE.0.) TR=237.3/(7.5/A-1.)

          IF(T.LT.0.) TR=253./(9.02/A-1.)

          PRINT 5,D,PHU,TR

          WRITE(50,5) D,PHU,TR

  5      FORMAT(' D=',F6.2,' PHU=', F6.2,' TR=',F7.2)

          RETURN

          END

          SUBROUTINE UID(IO)

          COMMON/DAN/D0,DO,N,P,TR

          PRINT 1

          DO 99 I=1,10

          PRINT 2,D0,DO,N,P

          PRINT 3

          READ 4,IO

          IF(IO.EQ.0) GOTO 100

          IF(IO.EQ.1) GOTO 100

          PRINT 6,D0,DO,N,P

          READ 8,D0,DO,N,P

  99    CONTINUE

 100 CONTINUE

C     FORMATNYJ BLOK

  1   FORMAT(1X,60(1H-),1H./

   *'  ПРОГРАММА УТОЧНЕНИЯ ИСХОДНЫХ ДАННЫХ UID     :'/

   *'    ПРОГРАММА ЗАПРАШИВАЕТ УПРАВЛЯЮЩИЙ СЧЕТЧИК IO, ОБОЗНАЧАЮЩИЙ ::'/

   *'   0 - ВЫХОД ИЗ ПРОГРАММЫ       :'/

   *'   1 - ПАРАМАТРЫ НЕ МЕНЯТЬ ;    :'/

   *'   2 - СМЕНИТЬ ПАРАМЕТРЫ ;       :'/

   *     1X,60(1H-),1H:)

  2      FORMAT(1X,61(1H*)/

   *'  ИСХОДНЫЕ ПАРАМЕТРЫ ТЕПЛОНОСИТЕЛЯ '/

   *'  D0,г DO,г N P,кПа   '/

   *  1X,F5.1,F5.1,I4,F6.1/

   *  1X,61(1H*))

  3      FORMAT(/' ВВЕДИТЕ ЖЕЛАЕМОЕ ЗНАЧЕНИЕ IO = '$)

  4      FORMAT(I2) 

  6      FORMAT(

   *'  ВВЕДИТЕ НОВЫЕ ЗНАЧЕНИЯ ПАРАМЕТРОВ '//

   *'  D0,г DO,г N P,кПа   '/

   *  1X,F5.1,F5.1,I4,F6.1/1H+$)

  7      FORMAT(1H+,' '$)

  8      FORMAT(F5.1,F5.1,I4,F6.1) 

          RETURN

          END

 

2.4. Варианты заданий

Номер

варианта

DO,

г/кг

D0,

г/кг

N

P,

кПа

1

2

3

4

5

6

7

8

9

10

0

0

0

0

0

0

0

0

0

0

20

20

20

20

20

20

20

20

20

20

100

100

100

100

100

100

100

100

100

100

25

50

75

100

125

150

175

200

225

250

Номер

варианта

DO,

г/кг

D0,

г/кг

N

P,

кПа

11

12

13

14

15

16

17

18

19

20

0

0

0

0

0

0

0

0

0

0

20

20

20

20

20

20

20

20

20

20

100

100

100

100

100

100

100

100

100

100

275

300

325

350

375

400

425

450

475

500