c********************************************************************
	
	SUBROUTINE t_bvf(v,dv,b,db,fesh,dfesh,teff,dteff)
		
c calcul de la Teff en utilisant les Mv, Mb et [Fe/H]
c suivant Magain 1987 A&A 181, 323

c routine du module mod_photo

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

c entrées
c	v, dv, b, db: magnitudes absolues V et B et erreur
c	fesh, dfesh : [Fe/H] et erreur

c sorties
c	teff, dteff : Teff et erreur

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

	USE mod_kind
	
	IMPLICIT NONE
	
	REAL (kind=dp), intent(in) :: v, dv, b, db, fesh, dfesh
	REAL (kind=dp), intent(out) :: teff, dteff
	
	REAL (kind=dp), SAVE :: ln10
	REAL (kind=dp) :: bmv, dteff_bmv, dteff_fesh

	LOGICAL, SAVE :: init=.TRUE.

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

	IF(init)THEN	!initialisations
	 init=.FALSE.
	 WRITE(*,10)
10	 FORMAT('calcul de la Teff a partir de B&V et [Fe/H]',/,
	1 '!!!!Attention!!!! pour naines de type solaire',
	2 'suivant Magain 1987 A&A 181, 323',/)
	
	 ln10=LOG(10.d0)	
	ENDIF
	
	bmv=b-v ; teff=7950.d0-4230.d0*bmv*(1.d0-0.204d0*10.d0**fesh)	
	dteff_bmv=4230.d0*(1.d0-0.204d0*10.d0**fesh)
	dteff_fesh=4230.d0*bmv*0.204d0*10.d0**fesh*ln10
	dteff=sqrt(55.d0**2+(dteff_bmv*db)**2+(dteff_bmv*dv)**2
	1 +sqrt(dteff_fesh*dfesh)**2)
				
	RETURN
	
	END SUBROUTINE t_bvf