c*******************************************************************************
c ... get cp for each of the pure species at given (i,j)
c ... dimensional

	subroutine getcpsp

	IMPLICIT NONE

        include 'parameter.par'
	include 'common.com'

	integer j, k

	double precision vtemp

c-------------------------------------------------------------------------------
c ... get cp/R of species 

	do j=1,n
	   vtemp = temp(j)
           call getcpor(vtemp,nsp,cpspj)
	   do k=1,nsp
	      cpsp(j,k) = cpspj(k)*rg(k)
	   end do
	end do

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

	return
	end