c*******************************************************************************
c ... find non-dimensional specific heat of the mixture at (j)
c ... assumes cpsp(j,k) have been already computed
c ... uses: cp = SUM(Y_k*cp_k) where cp_k are non-dimensional w.r.t. cpref
c ... if in_n2 is .true. then use a pure N2 gas
c ... otherwise, use the local 1D flame mixture fractions at each point
c ... assumes species(nsp) = N2

	function getcp()

	include 'parameter.par'
#include "precision.com"
	include 'common.com'

c-------------------------------------------------------------------------------

	do j=1,n

	   if(in_n2)then
	      do k=1,nsp-1
		 xsp(k) = 0.0
	      end do
	      xsp(nsp)=1.0
	   else
	      do k=1,nsp
	         xsp(k) = xs(j,k)
	      end do
	   endif

	   call getmsfr(nsp,xsp,ysp)

	   sum = 0.0
	   do k=1,nsp
	      sum = sum + ysp(k)*cpsp(j,k)
	   end do

	   cp(j) = sum

	end do

c-------------------------------------------------------------------------------

	return
	end