c******************************************************************
	
	REAL (kind=dp) FUNCTION fesh(h1,h2,he3,he4)
	
c function PUBLIC du module mod_spectroscopie

c calcul de [Fe/H]

c entrées: 
c	h1,h2,he3,he4 : abondances par gramme de H1, H2, He3, He4

c sortie
c	[Fe/H]

c Auteur: P. Morel, Laboratoire Lagange, O.C.A., CESAM2k

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

	USE mod_kind
	
	IMPLICIT NONE
	
	REAL (kind=dp), INTENT(in) :: h1, h2, he3, he4

	REAL (kind=dp), SAVE :: zsx_sol 	
	REAL (kind=dp) :: x, y, z

	LOGICAL, SAVE :: init=.TRUE.
	
c-----------------------------------------------------------------------

	IF(init)THEN
	 init=.FALSE.
	 WRITE(*,1)
1	 FORMAT('saisir Z/X_solaire, Exemple: Z/X=0.0245')
	 READ*,zsx_sol
	 WRITE(*,2)zsx_sol
2	 FORMAT('On utilise Z/X_solaire=',es10.3)
	ENDIF

	x=h1+h2 ; y=he3+he4 ; z=1.d0-x-y 
	IF(MIN(x,z) > 0.d0)THEN
	 fesh=LOG10(z/x/zsx_sol)
	ELSE		
	 fesh=0.d0
	ENDIF
		
	RETURN
	
	END FUNCTION fesh