c*******************************************************************************
c ... get cpspl(k) the specific heat of species k at (i,j)  
c ... linear interpolation from table of cpsp_k(T)

	subroutine cpsp(tc,cpspl)

	IMPLICIT NONE

        include 'parameter.par'
	include 'therm.com'

	integer jtab,jtp,k

	double precision tc,rtab,cpv

	real*8 cpspl(nsmx)

	jtab = int((tc-temp1tab)*odtab) + 1
	rtab = (tc-temptab(jtab))*odtab
	jtp  = jtab + 1

	do k=1,n_species
	   cpv      = cpsptab(k,jtab)
	   cpspl(k) = cpv + rtab*(cpsptab(k,jtp)-cpv) 
	end do

	return
	end