c*******************************************************************************
c ... get gsp for each of the pure species at given (i,j)
c ... gsp() is the standard state gibbs free energy of each species
c ... in SI units (J/kg)

	subroutine getgsp

	IMPLICIT NONE

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

	integer j, k

	double precision vtemp

c-------------------------------------------------------------------------------
c ... get gsp of each species

	do j=1,n
	   vtemp = temp(j)
           call getgibbs(vtemp,nsp,gspj)
	   do k=1,nsp
	      gsp(j,k) = gspj(k)
	   end do
	end do

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

	return
	end