c***********************************************************

	REAL (kind=dp) FUNCTION calc_entropie(p,t,xchim)

c routine public du module mod_etat

c calcul de l'entropie spécifique par etat_irwin en utilisant la structure EOS1
c voir les commentaires dans la routine etat_irwin lignes 91 à 103

c entrées
c	p, t, xchim : pression, température, comp.chim (par gramme)

c Auteur: P.Morel, laboratoire Lagrange, O.C.A.

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

	USE mod_kind

	IMPLICIT NONE

	REAL (kind=dp), INTENT(in), DIMENSION(:) :: xchim
	REAL (kind=dp), INTENT(in) :: p, t

	REAL (kind=dp) :: alfa, beta, cp, dcpp, dcpt, dcpx, delta,
	1 deltap, deltat, dgradadp, dgradadt, dgradadx, deltax, drop,
	2 drot, drox, dup, dut, dux, entropie, gamma1, gradad, ro, u

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

2000	FORMAT(8es10.3)

	CALL etat_irwin(p,t,xchim,.FALSE.,
	1 ro,drop,drot,drox,u,dup,dut,dux,
	1 delta,deltap,deltat,deltax,cp,dcpp,dcpt,dcpx,entropie,
	2 gradad,dgradadp,dgradadt,dgradadx,alfa,beta,gamma1)
c	WRITE(*,2000)entropie ; PAUSE'entropie'

c sortie
	calc_entropie=entropie

	RETURN

	END FUNCTION calc_entropie