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

	SUBROUTINE dnu0_int(r,vson,n,dnu0)
	
c routine du module mod_sismologie
	
c calcul de delta nu0 par l'intégrale de l'inverse de la vitesse du son

c entrées: 
c	r : rayons
c	vson : vitesse du son
c	n : nombre de couches
c sortie:
c	dnu0: delta nu0
	
c Auteur: P. Morel, Departement Lagrange, O.C.A., CESAM2k

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

	USE mod_kind
	USE mod_numerique, ONLY : bsp1dn, sum_n
	
	IMPLICIT NONE
	
	REAL (kind=dp), INTENT(in), DIMENSION(:) :: r, vson 
	INTEGER, INTENT(in) :: n
	REAL (kind=dp), INTENT(inout) :: dnu0
	
	INTEGER, PARAMETER :: m=4
	REAL (kind=dp), DIMENSION(1,n) :: unsvson
	REAL (kind=dp), DIMENSION(n+m) :: rt
	REAL (kind=dp), DIMENSION(1) :: ff, dff
	
	INTEGER :: knot, l=m
	
c--------------------------------------------------------------------	

	unsvson(1,:)=1.d0/vson
	
	CALL bsp1dn(1,unsvson,r,rt,n,m,knot,.FALSE.,r(1),l,ff,dff)
	CALL sum_n(1,unsvson,rt,m,knot,.FALSE.,r(1),r(n),ff)
	dnu0=0.5d0/ff(1)
	
	RETURN
	
	END SUBROUTINE dnu0_int