c*******************************************************************************
c ... get hspl(k) the enthalpy of species k at (i,j)  
c ... linear interpolation from table of h_k(T)

	subroutine hsp(tc,hspl)

	IMPLICIT NONE

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

	integer jtab,jtp,k

	double precision tc,rtab,hsv


	real*8 hspl(nsmx)

	jtab = int((tc-temp1tab)*odtab) + 1
	rtab = (tc-temptab(jtab))*odtab
	jtp  = jtab + 1
	
	do k=1,n_species
	   hsv     = hsptab(k,jtab)
	   hspl(k) = hsv + rtab*(hsptab(k,jtp)-hsv) 
	end do

	return
	end