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

      SUBROUTINE acc_rad2_bp(lum,ray,t,nel,grav,ychim,ioni,g_rad,dg_rad)

c routine du module mod_bp_for_alecian_new

c     Version 0.0 (2004/09/22)
c     cette subroutine calcule l'accélération radiative et ses dérivées
c     par rapport aux l'abondances. Elle doit être appelée pour chaque
c     couche, et à chaque fois que la diffusion est calculée.
c
c Adaptation a CESAM2k par B. Pichon, OCA
c
c entrées:
c     ray: rayon
c     t : température
c     nel : nombre électrons par volume,
c     ychim : composition chimique/mole
c     ioni : taux d'ionisation
c     grav : gravité
c
c sorties
c     g_rad(i) : accélérations radiatives la gravité est grav+g_rad sur i
c     dg_rad(i,j) : dérivée de g_rad(i) / ion j
c
c--------------------------------------------------------------

	USE mod_donnees, ONLY: lsol, mtot, nchim, nucleo, rsol
	USE mod_kind
	USE mod_numerique, ONLY : entre

	IMPLICIT NONE

      REAL (kind=dp), INTENT(in) ::  grav, lum, nel, ray, t
      REAL (kind=dp), INTENT(in), DIMENSION(:) :: ychim
      REAL (kind=dp), INTENT(in), DIMENSION(0:,:) :: ioni

      REAL (kind=dp), INTENT(out), DIMENSION(:,:) :: dg_rad
      REAL (kind=dp), INTENT(out), DIMENSION(:) :: g_rad

      LOGICAL, SAVE :: init= .TRUE.

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

2000	FORMAT(8es10.3)

	IF(init)THEN
         CALL ppxi2_bp(mtot,nchim)
         init=.FALSE.
	ENDIF

	CALL grad2_bp(grav,lum*lsol,ray*rsol,t,nel,ychim,ioni,g_rad,dg_rad)

	RETURN

	END SUBROUTINE acc_rad2_bp