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

	subroutine gethsp

	IMPLICIT NONE

        include 'parameter.par'
	include 'common.com'

	integer j, k

	double precision vtemp

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

	do j=1,n
	   vtemp = temp(j)
           call gethort(vtemp,nsp,hspj)
	   do k=1,nsp
	      hsp(j,k) = hspj(k)*rg(k)*vtemp
	   end do
	end do

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

	return
	end