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

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

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


5.2.3. вывод результатов расчета

Вывод результатов расчета производится подпрограммой PR, текст которой приведен в п. 2 и используется во всех последующих работах. В процессе ее работы непосредственно на экран выводится справочная информация, сообщающая о проценте заполнения файла результатов R5U.REZ в виде: FILE R5U.REZ: 100 \%, где значения цифр быстро изменяются от нуля до N.

Результаты расчета помещаются в специальный файл результатов R5U.REZ, который создается автоматически. После окончания работы с программой в этом файле можно будет просмотреть все просчитанные варианты в хронологическом порядке. Выходные данные представлены в удобной для чтения табличной форме, формат которой задается в подпрограмме SUBROUTINE PR.

5.3. Текст программы (файл R5U.FOR)

 

C       НАГРЕВ ВЛАЖНОГО ВОЗДУХА R5U.FOR

          BLOCK DATA

          COMMON/DAN/T0,TO,N,D0,P,IRT

          DATA T0,TO,N,D0,P,IRT/-50.,50.,100,20.,101.325,1/

          END

          COMMON/DAN/T0,TO,N,D0,P1,IRT

          COMMON/REZ/B(7,100)

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

  11    CALL UID(IO)

          IF(IO.EQ.0) GOTO 10

          TSH=(TO-T0)/N

          CALL OX1(TSH)

          CALL PR

          GOTO 11

  10    CONTINUE

          STOP

          END

          SUBROUTINE OX1(TSH)

          COMMON/DAN/T0,TO,N,D0,P1,IRT

          COMMON/REZ/B(7,100)

          T1=T0

          P2=P1

          DO 1 I=1,N

          T2=T1+TSH

          CALL FCPV(T1,P1,T2,P2,D0,IRT,CPX)

          CPV=CPB(T1,P1,D0,IRT)

          Q=(T2-T1)*CPV

          TK=T1+Q/CPX

          B(1,I)=I*1.

          B(2,I)=T0+TSH*I

          B(3,I)=CPV

          B(4,I)=CPX

          B(5,I)=Q

          IF(IRT.EQ.2) GOTO 11

          IF(T1.GE.0.) GOTO 11

          IF(T2.LE.(-0.1)) GOTO 11

          CALL UTK(I,T1,T2,TK)

          PRINT 6

  11    CONTINUE

          B(6,I)=TK

          B(7,I)=B(2,I)-TK

          T1=TK

          PRINT 5,I

  1      CONTINUE

  5      FORMAT(1H+,' GO R5U :',I4,' \%'$)

  6      FORMAT(' GO R5U :')

          RETURN

          END