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

	subroutine gethsp

        include 'parameter.par'
#include "precision.com"
	include 'common.com'

c-------------------------------------------------------------------------------
c ... get h/RT of species

	do j=1,n
	   vtemp = temp(j)
           call gethort(vtemp+epso2,nsp,hspj_p)
           call gethort(vtemp-epso2,nsp,hspj_m)
	   do k=1,nsp
	      hspj_p(k) = hspj_p(k) * rg(k) * (vtemp+epso2)
	      hspj_m(k) = hspj_m(k) * rg(k) * (vtemp-epso2)
	   end do
	   do k=1,nsp
	      hsp(j,k) = ( hspj_p(k) - hspj_m(k) ) / epstemp
	   end do
	end do

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

	return
	end