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

	REAL (kind=dp) FUNCTION initialise_u(r,rstar)
	
c	function private du module mod_evol
c	valeur initiale de la vitesse angulaire

c entrée :
c	r : rayon en Rsol
c	rstar : rayon total en Rsol

c sortie :
c	u_initial : vitesse de circulaion méridienne

c	Auteur: P.Morel, Département Cassiopée, O.C.A.
c	CESAM2k

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

	USE mod_donnees, ONLY : w_form, w_rot
	USE mod_kind

	IMPLICIT NONE
	
	REAL (kind=dp), INTENT(in) :: r, rstar
	REAL (kind=dp), PARAMETER :: u0=1.d-4
		
c--------------------------------------------------------------------------------

2000	FORMAT(8es10.3)

c valeur initiale de la vitesse de circulation méridienne
c	IF(w_form == 0.d0)THEN
c	 u_initial=0.d0
c	ELSE
c	 u_initial=ABS(w_rot)/(2.d0-(r/rstar)**2)*u0	
c	 u_initial=2.d0*ABS(w_rot)/(1.d0+(r/rstar)**2)*u0
c	ENDIF

	initialise_u=ABS(w_rot)/(2.d0-(r/rstar)**2)*u0

	RETURN
	
	END FUNCTION initialise_u