\section{Test de fonctionnement}\index{test de fonctionnement}
\label{sec:btst}
Se placer dans le sous-directory SOURCE de \cesam2k.
Créer un module exécutable du programme {\tt cesam2k}
en utilisant la
procédure {\tt exe2k}. Taper 0 (zéro) "pour arrêter". On obtient:

\begin{verbatim}
morel@pcmorel% exe2k cesam2k
Encountered 0 errors, 0 warnings in file CESAM2k_noplot.f.

 Pour arreter : taper 0 puis RETURN

 Pour poursuivre une evolution : taper 1 puis RETURN

 Pour initialiser un modele de ZAMS : taper 2 puis RETURN

 Pour initialiser un modele de PMS : taper 3 puis RETURN
0
  ARRET

\end{verbatim}
Ci-dessus, on a entré 0 (zéro) ``pour arrêter''.
Se placer ensuite dans le directory {\sc tests} pour y lancer
l'exécution. On suppose, pour ce faire, que le fichier
{\tt .cshrc} indique le
chemin du directory {\tt CESAM2k/SOURCE}, ainsi qu'il est suggéré au
\paragr{sec:impl}. Si tel n'est pas le cas, copier l'exécutable
{\tt CESAM2k.out} dans le directory {\tt tests}.
Vérifier aussi que le fichier {\tt test.don} est bien l'exemple
décrit au \paragr{sec:fdon}. On obtient:
\begin{verbatim}
morel@pcmorel% CESAM2k_noplot.out

 Pour arreter : taper 0 puis RETURN

 Pour poursuivre une evolution : taper 1 puis RETURN

 Pour initialiser un modele de ZAMS : taper 2 puis RETURN

 Pour initialiser un modele de PMS : taper 3 puis RETURN
 \end{verbatim}
 entrer {\bf\tt 2}, on obtient:
 \begin{verbatim}
2
 le modele initial de ZAMS est-il donne en binaire ? o/n
 \end{verbatim}
 entrer {\bf\tt n}, on obtient:
 \begin{verbatim}
 entrer le nom du fichier ASCII du modele initial
 Exemples: m010.zams, m020.zams, m050.zams
 \end{verbatim}
 entrer {\bf\tt m010.zams}, on obtient:
 \begin{verbatim}
 CESAM utilise le modele m010.zams

 entrer l'identificateur du modele
 Exemple: mon_modele, modele_euler
 \end{verbatim}
 entrer {\bf\tt test}, on obtient:
 \begin{verbatim}

 identificateur des fichiers du modele : test

             ***************************************************************

              MODELE DE STRUCTURE INTERNE calcule par CESAM version 5.0.0.0

             ***************************************************************

                        Debut du calcul le: 08 Janvier 2003 a 17h07

         NAMELISTS du fichier: test.don

 &NL_NOMS NOM_CTES_PHYS=ctes_phys85,
 NOM_DES=no_des,
 NOM_PERTE=perte_ext,
 NOM_ATM=lim_atm,
 NOM_TDETAU=edding,
 NOM_ABON=solaire_gn,
 NOM_DIFFM=diffm_mp,
 NOM_DIFFT=difft_nu,
 NOM_CONV=conv_jmj,
 NOM_ETAT=etat_eff,
 NOM_OPA=opa_int_zsx,
 NOM_NUC=ppcno9 ,
 NOM_NUC_CPL=NACRE,
 NOM_OSC=osc_adia,
 NOM_CHEMIN=/home/morel/SUN_STAR_DATA/
 /
 &NL_ESPACE MTOT=1.000000000000000,
 MDOT=0.000000000000000E+00,
 DER_NUM=F,
 PRECISION=sp,
 N_MAX=1500
 /
 &NL_ATMOSPHERE TAU_MAX=10.00000000000000,
 LIM_RO=T
 /
 &NL_TEMPS AGEMAX=50.00000000000000,
 DTLIST=1000000000000.000,
 X_STOP=-0.1000000000000000,
 LOG_TEFF=10.00000000000000,
 T_STOP=50000000.00000000,
 ARRET=else
 /
 &NL_CHIM X0=0.7000000000000000,
 Y0=0.2800000000000000,
 ZSX0=0.000000000000000E+00,
 D_TURB=0.000000000000000E+00,
 RE_NU=1.000000000000000,
 DIFFUSION=F,
 F_RAD=F
 /
 &NL_ROT W_ROT=0.000000000000000E+00,
 DW=0.000000000000000E+00,
 ROT_SOLID=T
 /
 &NL_CONV ALPHA=1.800000000000000,
 CPTURB=0.000000000000000E+00,
 OVSHTS=0.000000000000000E+00,
 OVSHTI=0.000000000000000E+00,
 JPZ=F,LEDOUX=F
 /
 &NL_ETAT F_EOS=eos_opal_190.bin
 &NL_OPA F_OPA=opa_yveline_etendu.bin v9/OPINTPATH_AX
 /
 &NL_NUC MITLER=F
 /

    Constantes fondamentales selon Handbook ctes. phys. + GONG + NACRE

 masse solaire= 1.98910000E+33
 luminosite solaire= 3.84600000E+33
 rayon solaire= 6.95990000E+10
 charge de l'electron= 4.80320680E-10
 constante de Boltzmann= 1.38065800E-16
 constante de Planck= 6.62607550E-27
 constante gravitation= 6.67259000E-08
 constante radiation= 7.56591220E-15
 celerite de la lumiere= 2.99792458E+10
 electron volt= 1.60217733E-12
 masse atomique unite= 1.66054020E-24
 masse de l'electron= 9.10938970E-28
 zero des Mbol L= 3.05500000E+33
 constante de Newton G Msol= 1.32724488E+26

 masses des noyaux, donnees de NACRE 1999

    H= 1.00782500E+00     H2= 2.01401800E+00      He3= 3.01602930E+00
  He4= 4.00260330E+00    Li6= 7.01600400E+00      Li7= 6.01512100E+00
  Be7= 7.01692920E+00    Be9= 9.01218210E+00      C12= 1.20000000E+01
  C13= 1.30033548E+01    N13= 1.30057386E+01      N14= 1.40030740E+01
  N15= 1.50010890E+01    O16= 1.59949146E+01      O17= 1.69991315E+01
  O18= 1.79991640E+01    F18= 1.80009377E+01      F19= 1.89984032E+01
 Ne20= 1.99924402E+01   Ne22= 2.19913855E+01     Na23= 2.29897697E+01
 Mg24= 2.39850419E+01   Mg25= 2.49858370E+01     Mg26= 2.59825930E+01
 Fe56= 5.58470000E+01   neut= 1.00866500E+00     prot= 1.00727647E+00

             Physique utilisee

EOS: etat_eff                Opa.: opa_int_zsx             conv.: conv_jmj
reac. nuc.: ppcno9           compil. nuc.: NACRE           abon_ini.: solaire_gn
atmos.: lim_atm              T(tau): edding                pertes: perte_ext
diff. micro.: diffm_mp       diff. turb.: difft_nu         ctes.: ctes_phys85

----------parametres numeriques utilises-----------

modele calcule en super precision

m_qs= 2, m_ch= 2, ordre= 2, precix= 1.0E-04, loc_zc= 1.0E-02
precit= 2.0E-01, psi0= 1.0E-01, mvt_dis= T

utilisation des variables lagrangiennes

 modele sans diffusion
 modele sans perte de masse

modele d'initialisation en ASCII: m010.zams
nombre de couches: 197
 Reactions nucleaires tabulees pour: ppcno9
 PP+CNO, 9 elements, H2, Li7, Be7 a l'equilibre

 reactions de la compilation NACRE

 elements chimiques utilises:
  H1
 He3
 He4
 C12
 C13
 N14
 N15
 O16
 O17

 reactions nucleaires utilisees:
 H1(p,e+ nu)H2
 H2(p,g)He3
 He3(He3,2p)He4
 He3(a,g)Be7
 Li7(p,a)He4
 Be7(e-,nu g)Li7
 Be7(p,a)He4
 C12(p,g)C13
 C13(p,g)N14
 N14(p,g)N15
 N15(p,g)O16
 N15(p,a)C12
 O16(p,g)O17
 O17(p,a)N14

ecrantage faible
Temperature d'amorcage des reactions nucleaires: 5.000E+05
Temperature maximale pour les reactions nucleaires: 1.000E+08

Mixture solaire de Grevesse & Noels 93

Rapports isotopiques suivant Anders & Grevesse
sauf He3/He4 de Gautier & Morel 97

Ex, element fictif # CNO : Si28
masse de l'element fictif :  28
charge de l'element fictif :  14

 Reactions thermonucleaires des cycles PP, CNO

 nombre de reactions :  14
 nombre d'elements chimiques :  10

abondances initiales deduites de X0= 7.000E-01, Y0= 2.800E-01, Z0= 2.000E-02
Z0/X0= 2.857E-02
H1=X0, H2+He3+He4=Y0, avec H2 dans He3
Z0 = 1-X0-Y0 = C12+C13+N14+N15+O16+O17+Ex

H1: 7.000E-01, He3: 8.827E-05, He4: 2.799E-01, C12: 3.425E-03, C13: 4.128E-05
N14: 1.059E-03, N15: 4.168E-06, O16: 9.642E-03, O17: 3.903E-06, Ex: 5.824E-03

rapports en nombre:
C/Z: 2.455E-01, N/Z: 6.458E-02, O/Z: 5.130E-01, Ex/Z: 1.769E-01

 abondances negligeables:
H1: 1.000E-03, He3: 5.000E-07, He4: 1.000E-03, C12: 5.000E-06, C13: 1.000E-07
N14: 1.000E-06, N15: 5.000E-09, O16: 1.000E-05, O17: 5.000E-09, Ex: 1.000E-06

 pour l'evolution temporelle, test de precision sur H et He4


rotation solide, vitesse angulaire initiale = 0.000E+00 rd/sec
modele sans rotation

------ Integration du modele quasi statique ---------

Calcul du Jacobien par derivees analytiques
saut de la fonction de repartition par couche = 1.000E-01
modif. nb. couches si variation cte.rep. > 2%
constantes de repartition utilisees:
ctep=-1.00E+00, ctem= 1.50E+01

modif du nb. couches pour mod. quasi-stat. de  300 a  437
|fct.rep.-psi0|/psi0= 4.578E-02 > dpsi= 2.000E-02

ecart des limites ZR/ZC au plus proche noeud: 1.000E-02

 modele sans overshooting

 ------------- La Convection (debut) -----------------

dans Z, masse elem. chim. hors network : 6.072E-18


equation d'etat EFF

opacites a Z variable des tables OPAL prolongees Kurucz:
klgx..z..., X=0.70, 0.35, 0.00; Z= 0.010, 0.019, 0.030
extrapolation lineaire en X= 0.80
compilation de Y. Lebreton et J. Fernandez
extrapolation et approximation par B-splines 4D de P. Morel
ordre des splines: mr, mt, mx, mz:  2,   2,   2,   2
 Z0= 2.000E-02
debut de ZC localisee a 52% entre les couches #    1 et #    2, limite ZR/ZC #  1
fin de ZC localisee a 24% entre les couches #   36 et #   37, limite ZR/ZC #  2
debut de ZC localisee a 11% entre les couches #  225 et #  226, limite ZR/ZC #  3
fin de ZC localisee a 96% entre les couches #  431 et #  432, limite ZR/ZC #  4
 suppression d'une limite ZR/ZC trop centrale # 1

 Limites ZR/ZC retenues apres affinement de la localisation:
  pour limite ZR/ZC #   1 a  23% entre couches #   36 et #   37, fin de ZC
   R_zc/Rtot= 1.381E-01, M_zc/Mstar= 9.248E-02, 1-M_zc/Mstar= 9.075E-01
  pour limite ZR/ZC #   2 a  12% entre couches  225 et  226, debut de ZC
   R_zc/Rtot= 7.364E-01, M_zc/Mstar= 9.776E-01, 1-M_zc/Mstar= 2.244E-02
  pour limite ZR/ZC #   3 a  95% entre couches #  431 et #  432, fin de ZC
   R_zc/Rtot= 9.999E-01, M_zc/Mstar= 1.000E+00, 1-M_zc/Mstar= 2.784E-11

 ------------- La Convection (fin) -----------------
Utilisation de la methode de collocation
coefficient pour Newton-Raphson modifie, nrm= 1.500E-01

Calcul exact du TdS

 ---------THERMO----------

 convection: critere de SCHWARZSCHILD

 sans pression turbulente

 gradient convectif calcule par MLT :
 formulation tenant compte de l'ep. op. de la bulle
 V/Al=  2.222E-01, ksi0= 6.173E-03, phi0= 2.250E+00

------------------ATMOSPHERE-----------------------------

loi t(tau,teff,grav) Eddington, loi purement radiative
tau_min=1.d-4, ro_ext=3.55d-9

atmosphere restituee avec raccord du gradient
tau au fond de l'atmosphere tau_max= 1.000D+01
rayon et masse pris a tau* tel que T(tau*,Teff)=Teff
tau exterieur 1.000D-04, nombre de couches:  75, R* a la couche 56
 limite externe portant sur la densite

 convection dans atmosphere: critere de SCHWARZSCHILD

 sans press. turb. dans atmosphere

 atmosphere iter.  1 err. max. 6.5E-02, variable 1, couche 75, corr 1.0E+00

 atmosphere iter.  2 err. max. 1.4E-02, variable 1, couche 74, corr 1.0E+00

 atmosphere iter.  3 err. max. 9.7E-04, variable 2, couche  1, corr 1.0E+00

 atmosphere iter.  4 err. max. 9.0E-06, variable 2, couche  1, corr 1.0E+00

 atmosphere iter.  5 err. max. 4.1E-10, variable 2, couche  1, corr 1.0E+00

 derivees partielles numeriques / r et l

 atmosphere iter.  1 err. max. 1.0E-05, variable 3, couche 74, corr 1.0E+00

 atmosphere iter.  2 err. max. 4.1E-12, variable 1, couche 65, corr 1.0E+00

 atmosphere iter.  1 err. max. 1.0E-05, variable 3, couche 75, corr 1.0E+00

 atmosphere iter.  2 err. max. 7.2E-12, variable 1, couche 34, corr 1.0E+00

 ------- L'atmosphere (fin) ------

iter. globale:  1, err. max: 1.3E-01, variable: 2, couche:   429, corr: 4.0E+00
variables Pgaz, T, R, L, M : 1.969E+05 1.041E+04 8.777E-01 6.927E-01 1.000E+00

modif du nb. couches pour mod. quasi-stat. de  437 a  631
|fct.rep.-psi0|/psi0= 4.458E-02 > dpsi= 2.000E-02

 ------------- La Convection (debut) -----------------

fin de ZC localisee a 25% entre les couches #   52 et #   53, limite ZR/ZC #  1
debut de ZC localisee a  2% entre les couches #  327 et #  328, limite ZR/ZC #  2
fin de ZC localisee a 46% entre les couches #  627 et #  628, limite ZR/ZC #  3
 suppression d'une limite ZR/ZC trop externe # 3

 Limites ZR/ZC retenues apres affinement de la localisation:
  pour limite ZR/ZC #   1 a  25% entre couches #   52 et #   53, fin de ZC
   R_zc/Rtot= 1.385E-01, M_zc/Mstar= 9.230E-02, 1-M_zc/Mstar= 9.077E-01
  pour limite ZR/ZC #   2 a   2% entre couches  327 et  328, debut de ZC
   R_zc/Rtot= 7.367E-01, M_zc/Mstar= 9.773E-01, 1-M_zc/Mstar= 2.266E-02

 ------------- La Convection (fin) -----------------

 atmosphere iter.  1 err. max. 2.9E-03, variable 3, couche 57, corr 1.0E+00

 atmosphere iter.  2 err. max. 7.8E-05, variable 1, couche 63, corr 1.0E+00

 atmosphere iter.  3 err. max. 2.3E-08, variable 1, couche 62, corr 1.0E+00

 derivees partielles numeriques / r et l

 atmosphere iter.  1 err. max. 1.0E-05, variable 3, couche 74, corr 1.0E+00

 atmosphere iter.  2 err. max. 4.3E-12, variable 1, couche 65, corr 1.0E+00

 atmosphere iter.  1 err. max. 1.0E-05, variable 3, couche 75, corr 1.0E+00

 atmosphere iter.  2 err. max. 7.8E-12, variable 1, couche 63, corr 1.0E+00

 ------- L'atmosphere (fin) ------

iter. globale:  2, err. max: 1.7E-01, variable: 3, couche:   630, corr: 2.0E+00
variables Pgaz, T, R, L, M : 1.617E+05 7.181E+03 9.262E-01 6.597E-01 1.000E+00

no modif. nb. couches pour var. quasi-stat. n=  631

..................................

..................................

 atmosphere iter.  1 err. max. 1.0E-05, variable 3, couche 74, corr 1.0E+00

 atmosphere iter.  2 err. max. 7.8E-12, variable 1, couche 62, corr 1.0E+00

 ------- L'atmosphere (fin) ------

iter. globale:  7, err. max: 5.2E-05, variable: 2, couche:   325, corr: 1.0E+00
variables Pgaz, T, R, L, M : 1.322E+14 2.504E+06 6.204E-01 6.498E-01 9.676E-01

 modele de ZAMS converge

 pas de dessin on line

 atmosphere iter.  1 err. max. 2.9E-05, variable 3, couche 74, corr 1.0E+00

 atmosphere iter.  2 err. max. 5.8E-11, variable 1, couche 62, corr 1.0E+00

 ------- L'atmosphere (fin) ------
*********
age= 0.000E+00, log Teff= 3.747E+00, log L/Lsol=-1.872E-01, log R/Rsol=-6.508E-02
Log10 g= 4.568E+00, Pc= 1.437E+17, Tc= 1.360E+07, Roc= 7.767E+01, Xc= 7.000E-01
en. PP= 94%, en. CNO=  6%, en. 3 alpha=  0%, en. grav=  0%, Yc= 2.799E-01
Var. rel. de masse : 0.000E+00, M*= 1.000E+00Msol, modele de la serie principale
*********
Modele suivant, pas temporel dt= 1.000E+01

no modif. nb. couches pour var. quasi-stat. n=  631


 ------------- La Convection (debut) -----------------

fin de ZC localisee a  1% entre les couches #   54 et #   55, limite ZR/ZC #  1
debut de ZC localisee a 54% entre les couches #  325 et #  326, limite ZR/ZC #  2

 Limites ZR/ZC retenues apres affinement de la localisation:
  pour limite ZR/ZC #   1 a   1% entre couches #   54 et #   55, fin de ZC
   R_zc/Rtot= 1.467E-01, M_zc/Mstar= 9.448E-02, 1-M_zc/Mstar= 9.055E-01
  pour limite ZR/ZC #   2 a  57% entre couches  325 et  326, debut de ZC
   R_zc/Rtot= 7.240E-01, M_zc/Mstar= 9.685E-01, 1-M_zc/Mstar= 3.152E-02

 ------------- La Convection (fin) -----------------

limites ZR/ZC utilisees pour le melange convectif

fin de ZC couche :   54
masse= 9.448D-02, 1-m/Mstar= 9.055D-01, rayon= 1.263D-01, rayon/Rstar= 1.467D-01

debut de ZC couche : 326
masse= 9.685D-01, 1-m/Mstar= 3.152D-02, rayon= 6.232D-01, rayon/Rstar= 7.240D-01

 nombre de couches pour la comp. chim. 467

 integration par Lobatto IIIc d'ordre 2

EVOL couche   1, var. max. de H ou He4: 1.409E-03
dt= 1.000E+01, dt optimal= 1.200E+01 masse / abondance / variation/scale / couche:
 0.0E+00 1.8E-01 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.0E+00 1.5E-01
 7.0E-01 2.1E-04 2.8E-01 2.8E-03 4.6E-04 1.4E-03 9.3E-08 9.6E-03 3.9E-06 5.8E-03
 6.6E-04 1.4E+00 1.4E-03 1.8E-01 1.0E+01 2.8E-01 9.8E-01 3.4E-06 1.0E-02 1.3E-17
     1      86       1       1       1       1       1       1       1      77

 atmosphere iter.  1 err. max. 6.2E-12, variable 1, couche 75, corr 1.0E+00

 derivees partielles numeriques / r et l

 atmosphere iter.  1 err. max. 1.0E-05, variable 3, couche 74, corr 1.0E+00

.......................................

.......................................

iter. globale:  7, err. max: 5.2E-03, variable: 4, couche:     3, corr: 1.0E+00
variables Pgaz, T, R, L, M : 1.461E+17 1.335E+07 2.299E-02 5.110E-03 6.973E-04
abon. -iers elem. : 6.967E-01 5.079E-05 2.827E-01 3.055E-04 1.094E-04 4.630E-03

 atmosphere iter.  1 err. max. 5.5E-05, variable 3, couche 74, corr 1.0E+00

 atmosphere iter.  2 err. max. 3.5E-09, variable 1, couche 75, corr 1.0E+00

 derivees partielles numeriques / r et l

 atmosphere iter.  1 err. max. 1.0E-05, variable 3, couche 74, corr 1.0E+00

 atmosphere iter.  2 err. max. 4.3E-12, variable 1, couche 65, corr 1.0E+00

 atmosphere iter.  1 err. max. 1.0E-05, variable 3, couche 75, corr 1.0E+00

 atmosphere iter.  2 err. max. 7.6E-12, variable 1, couche 34, corr 1.0E+00

 ------- L'atmosphere (fin) ------

iter. globale:  8, err. max: 3.5E-05, variable: 4, couche:     3, corr: 1.0E+00
variables Pgaz, T, R, L, M : 1.461E+17 1.335E+07 2.298E-02 5.106E-03 6.972E-04
abon. -iers elem. : 6.967E-01 5.079E-05 2.827E-01 3.055E-04 1.094E-04 4.630E-03

 modele converge


 atmosphere iter.  1 err. max. 7.6E-07, variable 3, couche 57, corr 1.0E+00

 atmosphere iter.  2 err. max. 7.7E-13, variable 1, couche 62, corr 1.0E+00

 ------- L'atmosphere (fin) ------

 creation du fichier pour le calcul des osc. adia.
 test.osc_adia
*********
age= 5.000E+01, log Teff= 3.749E+00, log L/Lsol=-1.754E-01, log R/Rsol=-6.125E-02
Log10 g= 4.560E+00, Pc= 1.484E+17, Tc= 1.344E+07, Roc= 8.134E+01, Xc= 6.967E-01
en. PP= 96%, en. CNO=  4%, en. 3 alpha=  0%, en. grav=  0%, Yc= 2.827E-01
Var. rel. de masse : 0.000E+00, M*= 1.000E+00Msol, modele de la serie principale
*********

Fin d'evolution avec CESAM version: 5.0.0.0
Sortie car agemax atteint
Fin du calcul le 08 Janvier 2003 a 17h11

nom du modele:
test

nom du fichier des NAMELISTs des donnees:
test.don

nom du fichier du listing des resultats en ASCII:
test.lis

nom du fichier du modele d'age 0 PMS en binaire:
test_B.pms

nom du fichier du modele d'age 0 homogene en binaire:
test_B.hom

nom du fichier du modele de ZAMS en binaire:
test_B.zams

nom du fichier du modele de POST en binaire:
test_B.post

nom du fichier du modele de COHE en binaire:
test_B.cohe

nom du fichier du modele de COCA en binaire:
test_B.coca

nom du fichier du modele final en binaire:
test_B.dat

nom du fichier du modele intermediaire en binaire:
test_B.rep

nom du fichier pour dessins du diag. HR et ZC :
test.HR

nom du fichier pour le calcul des osc. adia. :
test.osc_adia

nom du fichier pour le calcul des osc. non adia. :
test.osc_nadia

nom du fichier du modele d'atmosphere en binaire:
test_B.atm

Le code d'evolution stellaire CESAM a ete elabore
dans le cadre du Groupement de Recherche Structure Interne
des Etoiles et des Planetes Geantes. Si son utilisation
vous a donne satisfaction, le but poursuivi par tous
ceux qui y ont contribue aura ete atteint.

P. Morel, ON. Decembre 1989, CESAM1
P. Morel, OCA. Octobre 1991, CESAM2
P. Morel, OCA. Avril 1993, CESAM3
P. Morel, OCA. Decembre 1997, CESAM4
P. Morel, OCA. Decembre 2002, CESAM5

***********************************

morel@pcmorel%
\end{verbatim}
