c	INCLUDE '../SOURCE/mod_donnees.f'
	INCLUDE '../SOURCE_T/mod_atm.f'
c	INCLUDE '../SOURCE/mod_cesam.f'	
c	INCLUDE '../SOURCE/mod_nuc.f'	
	
c*******************************************************************

	PROGRAM test_atm

c	programme de test pour les routines d'atmosphère
c	CESAM2k

c	Auteur: P.Morel, Département J.D. Cassini, O.C.A.
c	CESAM2k

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

	USE mod_atm, ONLY : atm	
	USE mod_donnees, ONLY : ihe4, ini_ctes, lit_nl, mtot,
	1 m_qs, nchim, nom_elem, nom_fich2, nucleo, n_atm, precix,
	2 rep_atm, x0, y0
	USE mod_etat, ONLY : etat
	USE mod_kind
	USE mod_nuc, ONLY : nuc
	USE mod_numerique, ONLY : pause
	USE mod_variables, ONLY : mstar

	IMPLICIT NONE
	
	REAL (kind=dp), DIMENSION(:,:), ALLOCATABLE :: jac	 
	REAL (kind=dp), DIMENSION(:), ALLOCATABLE :: comp, dcomp, epsilo,
	1 ex, xchim

	REAL (kind=dp) :: be7, b8, dmsdl, dmsdr, dpsdl, dpsdr, dptsdl,
	1 dptsdr, dtsdl, dtsdr, et, ero, f17, hh, l_rac, m_rac, n13,
	2 o15, pt_rac, p_rac, r_rac, t_eff, t_rac, wrot
	
	INTEGER :: i

	LOGICAL :: list=.TRUE.

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

2000	FORMAT(8es10.3)

	PRINT*,'rep_atm0',rep_atm
	
c	nom_fich2='mon_modele'
c	nom_fich2='ksi_hya'		
 	nom_fich2='bid'
	  
	CALL lit_nl(wrot) ; CALL ini_ctes

	l_rac=6.545E-01 ; r_rac=8.6630E-01 ; precix=1.d-4
	n_atm=50 ; m_qs=1 ; mstar=mtot
	
c	l_rac=1.330E+01 ; r_rac=1.842E+00
	l_rac=6.448D+01 ; r_rac=2.558D+00
			
c	appel d'initialisation pour tabulation des réactions nucléaires
c	allocations fictives

	ALLOCATE(xchim(0),dcomp(0),jac(0,0),epsilo(0),ex(0))	
	CALL nuc(1.5d+07,1.5d+02,xchim,dcomp,jac,.FALSE.,0,
	1 epsilo,et,ero,ex,hh,be7,b8,n13,o15,f17)
	
c	détermination des abondances initiales

	DEALLOCATE(xchim) ; ALLOCATE(xchim(nchim))
	CALL nuc(1.5d+07,1.5d+02,xchim,dcomp,jac,.FALSE.,1,
	1 epsilo,et,ero,ex,hh,be7,b8,n13,o15,f17)
	xchim=xchim*nucleo
	WRITE(*,2000)xchim
		
	CALL atm(list,l_rac,r_rac,xchim,pt_rac,dptsdl,dptsdr,
	1 t_rac,dtsdl,dtsdr,m_rac,dmsdl,dmsdr,p_rac,dpsdl,dpsdr,t_eff)
	
	WRITE(*,2000)p_rac,t_rac

c	PAUSE'sortie'

	END PROGRAM test_atm