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

	PROGRAM test_degene

c test pour la dégénérescence
c on suppose le milieu totalement ionisé

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

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

	USE mod_donnees, ONLY : lit_nl, nchim, nom_fich2, nucleo, w_rot
	USE mod_kind
	USE mod_nuc, ONLY : nuc, t_nuc
		
	REAL (kind=dp), ALLOCATABLE, DIMENSION(:,:) :: jac		
	REAL (kind=dp), ALLOCATABLE, DIMENSION(:) :: comp, dcomp, ex
	REAL (kind=dp), DIMENSION(5) :: epsilon			
	REAL (kind=dp) :: be7e, b8e, et, ero,f17e, hhe, n13e, o15e, p, psi,
	1 t, ro

	INTEGER :: fait
	
	LOGICAL :: deriv=.FALSE.
	
c---------------------------------------------------------------

2000	FORMAT(8es10.3)

c lecture du fichier de données pour initialisations
	nom_fich2='1.0' ; CALL lit_nl(w_rot)

	ALLOCATE(jac(0,0),comp(0),dcomp(0),ex(0))

	fait=0
	CALL nuc(t,ro,comp,dcomp,jac,deriv,fait,
	1 epsilon,et,ero,ex,hhe,be7e,b8e,n13e,o15e,f17e)
	
	DEALLOCATE(comp,dcomp) ; ALLOCATE(comp(nchim),dcomp(nchim))

c calcul des abondances initiales
	fait=1
	CALL nuc(t,ro,comp,dcomp,jac,deriv,fait,
	1 epsilon,et,ero,ex,hhe,be7e,b8e,n13e,o15e,f17e)
	
	DEALLOCATE(jac,dcomp,ex)

	t=4.117d+07 ; ro=2.440d+05 ; p=3.019d+21	
        comp=(/ 1.294d-27,3.078d-60,9.797d-01,4.915d-05,1.348d-05,8.691d-03,
	1 3.690d-07,3.955d-03,1.722d-03,1.969d-03,2.225d-03,4.782d-07,
	2 7.470d-04,8.163d-04,6.524d-05 /)
	comp=comp/nucleo
	
	WRITE(*,2000)t,t_nuc(comp,p,ro,t)
	
	p=1.419d+17 ; t=1.352d+07 ; ro=7.717d+01 	
        comp=(/ 7.000d-01, 8.827d-05, 2.799d-01, 3.432d-03, 4.137d-05,
	1 1.061d-03, 4.177d-06, 9.642d-03, 3.903d-06, 1.959d-03, 2.214d-03,
	2 2.175d-05, 7.433d-04, 8.123d-04, 6.491d-05 /)
 	
	comp=comp/nucleo
	
	WRITE(*,2000)t,t_nuc(comp,p,ro,t)
		

	STOP
	
	END PROGRAM test_degene