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

	subroutine getcpsp

        include 'parameter.par'
#include "precision.com"
	include 'common.com'
	
c-------------------------------------------------------------------------------
c ... get cp/R of species 

	do j=1,n
	   vtemp = temp(j)
           call getcpor(vtemp+epso2,nsp,cpspj_p)
           call getcpor(vtemp-epso2,nsp,cpspj_m)
	   do k=1,nsp
	      cpsp(j,k) = rg(k) * ( cpspj_p(k) - cpspj_m(k) ) / epstemp
	   end do
	end do

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

	return
	end