c*************************************************************************-

      Module MOD_BP_FOR_ALECIAN_NEW

c Module des routines de calcul des accélérations radiatives

c Adaptation à CESAM2k B.Pichon, P.Morel

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

	USE MOD_DONNEES, only: msol, g, nomch => nom_chemin, zi, sigma, nchim,
	1 nom_elem, nucleo, lsol, rsol
	USE mod_kind

	IMPLICIT NONE

	INTEGER, PARAMETER :: nion_ppxi=30, pnchim = 13 , pzi=26, tgrad = 2

	REAL (kind=dp), DIMENSION(30,0:pzi,pnchim) :: niv_e,niv_q
	REAL (kind=dp), DIMENSION(0:nion_ppxi,pnchim) :: acstar, alphstar,
	1 betstar, phistar, psistar, xistar
	REAL (kind=dp), DIMENSION(0:pzi,pnchim) :: el_pot
	REAL (kind=dp), DIMENSION(pnchim) :: C_sol

	INTEGER, DIMENSION(30,0:pzi,pnchim) :: niv_g
	INTEGER, DIMENSION(pnchim,pnchim) :: isotops
	INTEGER, DIMENSION(0:pzi,pnchim) :: niv_flag, niv_nb, niv_z, rar_flag
	INTEGER, DIMENSION(pnchim) :: id_ppxi
        INTEGER  :: igrad  	!0 ou 1 cf. ppxi2_bp, g_rad pour élém. inconnus

	CHARACTER (len=2), DIMENSION(pnchim) :: el_svp

	SAVE

	PRIVATE
	PUBLIC :: acc_rad2_bp

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

      CONTAINS
     	INCLUDE 'acc_rad2_bp.f'
     	INCLUDE 'ft.f'
     	INCLUDE 'gr_gazrare.f'
     	INCLUDE 'gr_parti3.f'
     	INCLUDE 'gr_theta.f'
     	INCLUDE 'gr_popion3.f'
     	INCLUDE 'grad2_bp.f'
     	INCLUDE 'niv.f'
     	INCLUDE 'ppxi2_bp.f'

      End Module MOD_BP_FOR_ALECIAN_NEW