      PROGRAM COMBUST


      DIMENSION C(11,4)
      CHARACTER NAME(11,6)

      CP(T,T0,A,B,CC)=A+.5E-3*B*(T+T0)+.333333E-6*CC*(T*T+T*T0+T0*T0)

      OPEN(UNIT=1,FILE='POD')
      OPEN(UNIT=2,FILE='REZ')

C  SPECIFIC HEAT

      T0=273.16

      VM=22.4
      DT=100.

      DO 20 I=1,11
      READ (1,101) (NAME(I,K),K=1,6)
  101 FORMAT (6A1)
      READ (1,100) (C(I,K),K=1,4)
  100 FORMAT (10F8.4)
      IF(MOD((I-1),2).EQ.0) WRITE (2,102)
  102 FORMAT (1H1)
      WRITE (2,103) (NAME(I,K),K=1,6)
  103 FORMAT (/,1X,'COMPONENT  ',6A1)
      WRITE(2,200)
  200 FORMAT ('     TEMP    Cp        cp        cp         h        h')
      WRITE(2,201)
  201 FORMAT ('     st C  kJ/kmolK  kJ/m3K    kJ/kgK     kJ/m3   kJ/kg')

      T=273.16

   30 CP1=CP(T,T0,C(I,1),C(I,2),C(I,3))
      CP2=CP1/VM
      CP3=CP1/C(I,4)

      TC=T-273.16
      ENT1=CP2*TC
      ENT2=CP3*TC

      WRITE (2,110) TC,CP1,CP2,CP3,ENT1,ENT2
  110 FORMAT (F8.0,3F10.4,2F10.2,F10.0,F10.4)

      T=T+DT
      IF(TC.LT.2000.) GOTO 30

   20 CONTINUE

      CLOSE(1)
      CLOSE(2)

      STOP
      END

