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

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

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


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

C       РАСЧЕТ ВЛАГОСОДЕРЖАНИЯ R1.FOR

          BLOCK DATA

          COMMON/DAN/T0,TO,N,FN,P

          DATA T0,TO,N,FN,P/-50.,50.,100,100.,101.325/

          END

          COMMON/DAN/T0,TO,N,FN,P

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

  10    CALL UID(I0)

          IF(I0.EQ.0) GOTO 11

          DT=(TO-T0)/N

          WRITE(50,6) P

          NK=N+1

          DO 1 I=1,NK

          T=T0+DT*(I-1)

          CALL DA(P,T,FN,D,PH)

          PRINT 5,T,FN,D,PH

          WRITE(50,5) T,FN,D,PH

  1      CONTINUE

          GOTO 10

  11    CONTINUE

  5      FORMAT(' T=',F5.0,' FN=',F5.0,' D=',F10.5,' PH=',F10.5)

  6      FORMAT(' R1 : P=',F6.2,' КПА.')

          STOP

          END

          SUBROUTINE DA(P,T,FI,D,PH)

          IF(T.GE.0.) PH=0.611*EXP(2.30259*7.5*T/(237.3+T))

          IF(T.LT.0.) PH=0.611*EXP(2.30259*9.02*T/(253.+T))

          PH1=PH*FI/100.

          D=622.*PH1/(P-PH1)

          RETURN

          END

          SUBROUTINE UID(IO)

          COMMON/DAN/T0,TO,N,FN,P

          PRINT 1

PRINT 2,T0,TO,N,FN,P

          PRINT 3

          READ 4,IO

          IF(IO.EQ.0) GOTO 100

          IF(IO.EQ.1) GOTO 100

          PRINT 6,T0,TO,N,FN,P

          READ 8,T0,TO,N,FN,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*)/

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

    *'  T0,C TO,C N FN,\% P,кПа   '/

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

    *  1X,61(1H*))

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

  4      FORMAT(I2) 

  6      FORMAT(

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

   *'  Т0,C TO,C N FN,\% P,кПа   '/

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

  7   FORMAT(1H+,' '$)

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

          RETURN

          END