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

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

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

	do j=1,n
	   vtemp = temp(j)
           call getgibbs(vtemp+epso2,nsp,gspj_p)
           call getgibbs(vtemp-epso2,nsp,gspj_m)
	   do k=1,nsp
	      gsp(j,k) = ( gspj_p(k) - gspj_m(k) ) / epstemp
	   end do
	end do

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

	return
	end