function CPCO2(T)

      implicit none
      
      real CPCO2,T

c       Computes specific heat of CO2 using a fit to
c       data from www.engineeringtoolbox.com

        CPCO2 = 100.*(4.79 + 1.4416E-2*T - 7.2933E-6*T*T) !MKS
c
      return
      end