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

      PROGRAM test_der_rotx

c     test pour le calcul des variables pour le fichier des
c     oscillations non adiabatiques

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

c     CESAM5

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

      USE mod_cesam, ONLY : nbelem
      USE mod_chim, ONLY : nuc
      USE mod_donnees, ONLY : ini_ctes_phys, lit_nl, nom_fich2, print_ctes
      USE mod_kind
      USE mod_numerique, ONLY : pause

      IMPLICIT NONE

      REAL (kind=dp), ALLOCATABLE, DIMENSION(:,:) :: jac   
      REAL (kind=dp), ALLOCATABLE, DIMENSION(:) :: dcomp,
     1 epsilon, ex, xchim
      REAL (kind=dp) :: be7e, b8e,  dkapt, dkapro, dkapx, d2kap_tro,
     1 d2kap_ro2,d2kap_t2, dpdro, dpdt, dpdx, d2pdro2, d2pdrodt,
     2 d2pdt2, dudro, dudt, dudx, d2udro2, d2udrodt, d2udt2,
     3 ero, et, f17e, hhe, kap, n13e, o15e, p, ro, t

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

2000  FORMAT(8es10.3)

      CALL video(lignes=100,colonnes=80,blanc=.TRUE.)

c      CALL pause('entree')

c     lecture de mon_modele.don

      nom_fich2='mon_modele'
      CALL lit_nl ; CALL ini_ctes_phys ; CALL PRINT_ctes(6)

      p=1.9d5 ; t=8.d3 ; ro=1.d-8     !ro : fictif

      CALL nuc(t,ro,xchim,dcomp,jac,.FALSE.,0,
     1 epsilon,et,ero,ex,hhe,be7e,b8e,n13e,o15e,f17e)

	ALLOCATE(xchim(nbelem))
      CALL nuc(t,ro,xchim,dcomp,jac,.FALSE.,1,
     1 epsilon,et,ero,ex,hhe,be7e,b8e,n13e,o15e,f17e)
      WRITE(*,2000)xchim

      CALL pause('premier cas')
      WRITE(*,"('p=',es10.3,', t=',es10.3)")p,t
      CALL der_rotx(p,t,xchim,
     1 dpdro,dpdt,dpdx,d2pdro2,d2pdrodt,d2pdt2,
     2 dudro,dudt,dudx,d2udro2,d2udrodt,d2udt2,
     3 kap,dkapt,dkapro,dkapx,d2kap_tro,d2kap_ro2,d2kap_t2)    
      WRITE(*,*)'dpdro,dpdt,dpdx,d2pdro2,d2pdrodt,d2pdt2' 
      WRITE(*,2000)dpdro,dpdt,dpdx,d2pdro2,d2pdrodt,d2pdt2
      WRITE(*,*)'dudro,dudt,dudx,d2udro2,d2udrodt,d2udt2' 
      WRITE(*,2000)dudro,dudt,dudx,d2udro2,d2udrodt,d2udt2
      WRITE(*,*)'kap,dkapt,dkapro,dkapx,d2kap_tro,d2kap_ro2,d2kap_t2'     
      WRITE(*,2000)kap,dkapt,dkapro,dkapx,d2kap_tro,d2kap_ro2,d2kap_t2

      CALL pause('cas suivant')   
      p=1.843D+15 ; t=4.605D+06
      CALL der_rotx(p,t,xchim,
     1 dpdro,dpdt,dpdx,d2pdro2,d2pdrodt,d2pdt2,
     2 dudro,dudt,dudx,d2udro2,d2udrodt,d2udt2,
     3 kap,dkapt,dkapro,dkapx,d2kap_tro,d2kap_ro2,d2kap_t2)    
      WRITE(*,*)'dpdro,dpdt,dpdx,d2pdro2,d2pdrodt,d2pdt2' 
      WRITE(*,2000)dpdro,dpdt,dpdx,d2pdro2,d2pdrodt,d2pdt2
      WRITE(*,*)'dudro,dudt,dudx,d2udro2,d2udrodt,d2udt2'
      WRITE(*,2000)dudro,dudt,dudx,d2udro2,d2udrodt,d2udt2        
      WRITE(*,*)'kap,dkapt,dkapro,dkapx,d2kap_tro,d2kap_ro2,d2kap_t2'     
      WRITE(*,2000)kap,dkapt,dkapro,dkapx,d2kap_tro,d2kap_ro2,d2kap_t2

      STOP

      CONTAINS
      
       INCLUDE 'der_rotx.f'
       INCLUDE 'video.f'
      
      END PROGRAM test_der_rotx