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

	REAL (kind=dp) FUNCTION initialise_w(r,rstar)
	
c routine 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	w_initial : vitesse angulaire

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

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

	USE mod_donnees, ONLY : langue, pi, w_form, w_rot
	USE mod_kind

	IMPLICIT NONE
	
	REAL (kind=dp), INTENT(in) :: r, rstar
	
	LOGICAL, SAVE :: init=.TRUE.
		
c--------------------------------------------------------------------------------

2000	FORMAT(8es10.3)

	IF(init)THEN
	 init=.FALSE.

	 SELECT CASE(langue)
	 CASE('english')
	  WRITE(*,1000)w_form ; WRITE(2,1000)w_form
1000	  FORMAT('The angular velocity is initialized with q=',es10.3)
	 CASE DEFAULT	 
	  WRITE(*,1)w_form ; WRITE(2,1)w_form
1	  FORMAT('Initialisation de la vitesse angulaire, q=',es10.3)
	 END SELECT
	ENDIF

c valeur initiale de la vitesse angulaire
	
	initialise_w=ABS(w_rot)*(1.d0+w_form*(1.d0+COS(pi*r/rstar)))
	
	RETURN
	
	END FUNCTION initialise_w