c*******************************************************************************
c ... get G for a pure species
c ... linear interpolation from table of g_k(T)

	subroutine cgsp(tc,cgspl)

	IMPLICIT NONE

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


	double precision tc,cgsv,rtab,jtab,jtp

	integer k

	real*8 cgspl(nsmx)

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

	do k=1,n_species
	   cgsv    = cgsptab(k,jtab)
	   cgspl(k) = cgsv + rtab*(cgsptab(k,jtp)-cgsv)     
	end do

	return
	end