\section{Conseils pour l'exploitation}\index{conseils}
Un code numérique, aussi bien fait soit-il, ne peut être "clefs en mains".
Pour entamer une exploitation, on doit toujours être guidé par la
remarque triviale~:
\begin{verse}\begin{em}
Un programme qui ne fonctionne pas, ou mal, pour des cas simples \\
ne fonctionnera pas, ou encore plus mal, pour des cas compliqués.
\end{em}\end{verse}
Un code de structure interne n'échappe, hélas pas, à cette règle.
A cet effet, bien que \cesam2k ait prévu
différents types de ``précision'' standards, il est presque toujours
nécessaire,
à l'aide de ``réglages'', d'optimiser les paramètres numériques en
fonction des conditions rencontrées au cours du calcul.
Aborder un problème en introduisant {\em ad initio} toute la
complexité que l'on désire atteindre est une perte de temps.
Afin d'acquérir une expérience sur la façon dont l'évolution se déroule,
il faudra très rapidement reprendre des cas simples. Il est conseillé de
``dégrossir'' le problème en utilisant des options numériques et physiques
robustes et efficaces, afin d'apprécier rapidement ``ce qui se passe'',
puis d'introduire progressivement la complexité souhaitée.

\medskip
\PB{Les descriptions des routines données plus avant ne sont que succinctes,
se reporter aux commentaires figurants au début des routines pour plus d'informations.}

\section{Création d'un espace de travail}\label{sec:espt}
Pour un ensemble d'applications, il est pratique de créer un espace de
travail spécifique. \cesam2k offre la possibilité de simplifier
l'initialisation de l'environnement du
calcul.\index{create\_evol\_space}\index{espace de travail}
Par exemple, pour calculer un ensemble de modèles et effectuer des exploitations, de
w-virginis, on créera et on se placera dans un directory nommé {\tt W\_VIR}:

\begin{enumerate}
\item Saisir: {\tt create\_evol\_space W\_VIR}. Le directory {\tt W\_VIR} est créé.
Le {\tt README\_usage\_de\_evol\_space} décrit l'utilisation des exécutables qui y ont
été introduits. Les données générales (modèles
initiaux etc..) y sont recopiées de {\tt CESAM2k/EXPLOIT}. Se placer dans
{\tt W\_VIR} qui devient l'environnement des calculs.
\item Renommer le fichier de données {\tt mon\_modele.don} qu'on y trouvera:
{\tt cp mon\_modele.don w\_vir.don} et
adapter les paramètres. A cette fin, on trouvera dans
le fichier {\tt arguments\_du\_point\_don} les arguments
reconnus par la présente version \cesam2k. Ils seront détaillés ci-après.

\item Pour lancer un calcul saisir: {\tt cesam2k.out}
\end{enumerate}
Suivant la définition du moniteur, le graphique peut être tronqué ou mal
positionné. Il convient alors adapter la hauteur et la largeur des cadres en
bricolant dans le fichier {\tt device}, \cfparagr{sec:device}. Une fois des
réglages satisfaisants obtenus, les
reproduire dans le fichier {\tt device} du directory {\tt CESAM2k/EXPLOIT} afin
de les récupérer à chaque fois qu'un nouvel espace de travail sera créé par
{\tt create\_evol\_space}.

Le fichier {\tt mon\_modele.zoom} permet d'adapter, si besoin, les coordonnées
du dessin {\em on line} du diagramme HR, \cfparagr{sec:zoom}.


\section{Le fichier de données~: {\tt mon\_modele.don}}\label{sec:fichdon}
Les options physiques à utiliser et les données sont transmises
au programme principal par des {\tt NAMELIST} regroupées dans un fichier en
ASCII dont le nom a {\bf nécessairement} l'extension {\tt ".don},
par exemple~: {\tt mon\_modele.don}.
On a disposé un exemple de ce fichier dans le sous-directory {\tt EXPLOIT}.
\index{mon\_modele.don} \\

\begin{itemize}
 \item {\tt \& NL\_CESAM}~: regroupe les conditions générales du calcul.
 \index{NL\_CESAM}
 \item {\tt \& NL\_MASS}~: regroupe les paramètres physiques relatifs à
  la variable spatiale.\index{NL\_MASS}
 \item {\tt \& NL\_EVOL}~: regroupe les paramètres physiques relatifs à
  la variable temporelle.\index{NL\_TEMPS}
 \item {\tt \& NL\_CHIM}~: regroupe les paramètres physiques relatifs à
  la composition chi\-mi\-que.\index{NL\_CHIM}
  \item {\tt \& NL\_CONV}~: regroupe les paramètres physiques relatifs à
  la convection.\index{NL\_CONV}
 \item {\tt \& NL\_DIFF}~: regroupe les paramètres physiques relatifs à
  la diffusion.\index{NL\_DIFF}
 \item {\tt \& NL\_ROT}~: regroupe les paramètres physiques relatifs à
  la rotation.\index{NL\_ROT}
 \item {\tt \& NL\_ETAT}~: regroupe les noms des fichiers,
 {\em d'accès externe},
  éventuellement nécessaires au calcul de l'\EOS.\index{NL\_ETAT}
 \item {\tt \& NL\_OPA}~: regroupe les noms des fichiers,
 {\em d'accès externe}, éventuellement nécessaires au calcul de
  l'opacité.\index{NL\_OPA}
 \item {\tt \& NL\_NUC}~: définit des options relatives au calcul
  des taux de réactions thermonucléaires.\index{NL\_NUC}
 \item {\tt \& NL\_ATM}~: regroupe les paramètres physiques
  relatifs à la restitution de l'atmosphère.\index{NL\_ATM}
\end{itemize}

\section{Exemple de fichier de données}\label{sec:fdon}
Cet exemple est volontairement différent de celui utilisé au
\paragr{sec:tfonc}.
\index{nom\_ctes}\index{nom\_des}\index{nom\_perte}
\index{nom\_atm}\index{nom\_tdetau}\index{nom\_abon}
\index{nom\_diffm}\index{nom\_difft}\index{nom\_conv}
\index{nom\_etat}\index{nom\_opa}\index{nom\_nuc}\index{nom\_chemin}
\index{fichier de données}
Les types, real, integer... des paramètres numériques et logiques est facilement identifiable par les valeurs données ci dessous à
titre indicatif.
\begin{tt}\begin{verbatim}
 &NL_CESAM
 NOM_CHEMIN='~/CESAM2k/SUN_STAR_DATA/',
 NOM_CTES='ctes_94',
 NOM_DES='des_m',
 NOM_OUTPUT='no_outputR',
 N_MAX=6000,
 PRECISION='AV'
 /
 &NL_MASS
 MTOT=2.d0,
 NOM_PERTM='pertm_ext$',
 MDOT=0.d0
 /
 &NL_EVOL
 AGEMAX=12.d3,
 ARRET='else',
 DTLIST=1.d10,
 LOG_TEFF=10.d0,
 NB_MAX_MODELES=40,
 HE_CORE=-0.1d0,
 R_STOP=3.0d0,
 T_STOP=2.0d9,
 X_STOP=-0.1d0
 /
 &NL_CHIM
 GRILLE_FIXE=.FALSE.,
 NOM_ABON='solaire_gs',
 MODIF_CHIM=.FALSE.,
 GARDE_XISH=.FALSE.,
 X0=0.70,
 Y0=0.28,
 ZSX0=0.0d0
 /
 &NL_CONV
 NOM_CONV='conv_jmj',
 ALPHA=1.8d0,
 OVSHTS=0.d0,
 OVSHTI=0.d0,
 JPZ=F,
 CPTURB=0.d0,
 LEDOUX=F
 /
 &NL_DIFF
 DIFFUSION=.TRUE.,
 NOM_DIFFM='diffm_mp',
 NOM_DIFFT='difft_nu',
 D_TURB=10.d0,
 RE_NU=1.d0,
 NOM_FRAD='no_frad'
 /
 &NL_ROT
 W_ROT=0.d0,
 UNIT='kms/s',
 NOM_DIFFW='diffw_0',
 NOM_THW='rot_0',
 NOM_PERTW='pertw_0',
 NOM_DES_ROT='end_mod'
 /
 &NL_ETAT
 NOM_ETAT='etat_eff',
 F_EOS='eos_opal_250.bin',7*' '
 /
 &NL_OPA
 NOM_OPA='opa_houdek12',
 F_OPA='opa_yveline.bin','H12/OPINTPATH',6*' ',
 NOM_OPA_COND='opa_cond_mestel'
 /
 &NL_NUC
 NOM_NUC='ppcno3aco',
 NOM_NUC_CPL='NACRE',
 MITLER=F
 /
 &NL_ATM
 NOM_ATM='lim_atm',
 NOM_TDETAU='hopf',
 TAU_MAX=30.d0,
 LIM_RO=.TRUE.
 /
\end{verbatim}
\end{tt}
Certains arguments des NAMELISTs acceptent des\index{suffixes}
"SUFFIXES"\footnote{Toujours en "upper case".} permettant
différentes actions "secondaires", Exemple: avec le suffixe {\tt M} de l'item
{\tt PRECISION}
de la NAMELIST {\tt NL\_CESAM}, le calcul sera effectué avec le nombre de couches
fixé à la valeur maximale {\tt N\_MAX} indiquée dans le
fichier de données. Exemple: {\tt PRECISION='AVM'}

\noi S'il y a reprise ou poursuite d'une évolution, \cesam2k accepte que
des suffixes apparaissent/disparaissent des {\tt NAMELISTs}, leurs actions
seront alors prises en compte/ignorées dans l'évolution ultérieure.

\noi Ces SUFFIXES sont indiqués pour mémoire dans le fichier {\tt arguments\_du\_point\_don}
et leurs actions détaillées dans les paragraphes suivants,

\medskip

\PB{Les noms permettant des affectations sont {\embf case sensitive},
par exemple, il faut coder {\tt "NACRE"} et non {\tt "nacre"}; il
en est de même pour les chemins.}

\section{Arguments de la NAMELIST {\tt NL\_CESAM}}\label{sec:nlcesam}
\index{NL\_CESAM}\index{n\_max}\index{nom\_chemin}\index{routines disponibles}
\PB{Les arguments utilisables dans les NAMELISTs du fichier de données sont
indiqués dans le fichier {\tt argument\_du\_point\_don} des sous-directory
{\tt SOURCE} et {\tt EXPLOIT}.}
\medskip
\begin{itemize}
 \item{\tt NOM\_CHEMIN='$\sim$/SUN\_STAR\_DATA/'}~: chemin et nom du
 répertoire où sont
 disposées les données physiques, les tables d'opacité par exemple. Il est
parfois nécessaire de coder le chemin complet, Ex: {\tt
'/home/morel/CESAM2k/SUN\_STAR\_DATA/'}. Parfois {\tt '\$HOME/SUN\_STAR\_DATA/'}
n'est pas accepté.
 \index{nom\_ctes}\index{ctes\_85.f}\index{ctes\_94.f}\index{NL\_NOMS}
 \index{ctes\_94m.f}\index{ctes\_31}\index{ctes\_ba}\index{ctes\_gaia}
 \item {\tt NOM\_CTES}~: indique le nom de la routine d'initialisation des
 principales cons\-tan\-tes de physique à utiliser\footnote{Dans chaque cas,
 sont seules implémentées les différences connues avec les constantes
 par défaut de \cesam2k, \ie{\tt ctes\_94}.}~:
 \begin{itemize}
  \item {\tt ctes\_31}~: constantes du LATT, \citet{cvr07}.
  \item {\tt ctes\_85}~: constantes de {\sc gong}, \citet{c88}.
  \item {\tt ctes\_94}~: constantes de \cesam2k (défaut).
  \item {\tt ctes\_94m}~: identique à ctes\_94 avec des valeurs des masses
  des noyaux posées égales aux valeurs entières les plus proches.
  \item {\tt ctes\_ba}~: constantes de \citet{ba08}.
  \item {\tt ctes\_gaia}~: constantes de GAIA.
 \end{itemize}
 \item{\tt NOM\_DES}~: nom de la routine de dessin à utiliser~:
 \begin{itemize}\index{des\_m.f}\index{des\_r.f}\index{no\_des}\index{HHe}
 \index{routine de dessin}\index{des\_mC}
 \index{des\_r.f}\index{no\_des}\index{routine de dessin}
  \item {\tt des\_m/des\_r}~: dessin en fonction de la masse/rayon,
  utilise {\sc pgplot}. Dans le cadre des abondances, les valeurs indiquées sont
  normalisées par rapport à leur maximum sauf pour l'hydrogène et l'hélium
  indiquées en valeurs réelles.
   \item {\tt zoom}~: forme personnalisable de {\tt des\_m} décrite au
   \paragr{sec:zoom}, pratique lorsqu'on désire atteindre des conditions
   particulières, par exemple une localisation dans le diagramme HR. \\
  \noi{\bf Suffixes C, HHe, S}:\index{suffixes de dessin C, HHe \& S}
  \begin{description}
   \item[C]: {\tt des\_mC/des\_rC} permet d'indiquer, dans le cadre
   des abondances, les valeurs centrales
   en lieu et place des maximales (utile pour les stades avancés).
   \item[S]: {\tt des\_mS/des\_rS} permet d'indiquer les valeurs de surface
   des abondances en lieu et place des maximales (utile avec diffusion).
   \item[HHe]: \cesam2k formera un fichier, de nom
   {\tt mon\_modele.mHHe}, contenant les limites des shells source H et He,
   \cfparagr{sec:HHe}. Exemples: {\tt des\_mHHe/des\_rHHe}.
   \end{description}
  \item {\tt no\_des}~: sans dessin {\em on line}, à utiliser nécessairement
  si on ne dispose pas du logiciel {\sc pgplot}, voir \paragr{sec:impl}.
   \par
  \noi {\bf Suffixes 0, 1, 2} : {\tt no\_des0}~: sans dessin {\em on line}.
  {\tt no\_des1}~: graphes des variables quasi-statiques à chaque itération
  Newton-Raphson. {\tt no\_des2}~: graphes des variables quasi-statiques à l'issue de
  chaque pas temporel.
 \end{itemize}
 \label{sec:output1}
\item{\tt NOM\_OUTPUT}~: type de fichier ASCII à générer
 à l'issue du calcul~:\index{sorties ASCII}
 \begin{itemize}\index{routine de sorties ASCII}
  \index{nom\_output}\index{osc\_adia.f}\index{no\_output}
  \index{osc\_nadia.f}\index{mon\_modele-ad.osc}
  \index{mon\_modele-ascii}
  \index{mon\_modele-nad.osc}\index{mon\_modele-inv.osc}\index{osc\_invers.f}
  \index{fichier de sorties ASCII}\index{numérotation des modèles}
  \item {\tt osc\_adia}~: génération du fichier ASCII de nom
  {\tt mon\_modele-ad.osc} pour le calcul des oscillations adiabatiques et
  la réalisation de certains dessins,
  \item {\tt all\_adia}~: génération des fichiers ASCII,  \index{all\_adia}
  pour le calcul des oscillations adiabatiques et des dessins,
  pour {\embf TOUS} les modèles. Ces fichiers ont pour dénomination
  {\tt {\em nnnn}-mon\_modele-ad.osc}, {\em nnnn} étant le numéro du modèle,
  \item {\tt osc\_invers}~: génération des fichier ASCII de nom
  {\tt mon\_modele-inv.osc} pour les inversions,
  \item {\tt all\_invers}~: génération de {\embf TOUS} les fichiers ASCII
  de noms {\tt {\em nnnn}-mon\_modele-inv.osc} pour les inversions,
  \index{all\_invers} {\em nnnn} étant le numéro du modèle,
  \item {\tt osc\_nadia}~: génération du fichier ASCII de nom
  {\tt mon\_modele-nad.osc} pour le calcul des oscillations non adiabatiques,
  \item {\tt all\_nadia}~: génération de {\embf TOUS} les fichiers ASCII de noms
  {\tt {\em nnnn}-mon\_modele-nad.osc} pour le calcul des oscillations non
  adiabatiques, {\em nnnn} étant le numéro du modèle,\index{all\_nadia}
  \item {\tt no\_output}~: il n'y a pas de génération de fichier ASCII,
  \item {\tt ascii}~: forme à personnaliser du fichier ASCII
  décrite au \paragr{sec:ascii}, avec la dénomination {\tt mon\_modele-ascii}
  \item {\tt all\_ascii}~: génération de {\embf TOUS} les fichiers ASCII
  de forme à personnaliser de nom {\tt {\em nnnn}-mon\_modele-ascii},
  {\em nnnn} étant le numéro du modèle.\index{all\_ascii}
% \end{itemize}
% \item{\bf Suffixes:}
 \begin{description}
  \item{\bf Suffixes A, T, L, R :}
  \item [A] Le suffixe A (Ex: osc\_adiaA) permet d'utiliser
  l'ancienne formulation de la fréquence de Brunt-Väissälä.
  \item [T, L] Les suffixes T \& L (Ex: osc\_adiaT) permettent de générer des
  fichiers de sortie ASCII exploitables directement par IDL. Ils ont été créés
  à l'initiative de T.Guillot selon les desiderata respectivement de
  T.Guillot (TR) et L.Bigot (LB); \cfparagr{sec:tr},  \cfparagr{sec:lb}
{\embf Les suffixes T \& L sont incompatibles}.\index{suffixes de output A, T \& L}
\item [R] Le suffixe R (Ex: no\_outputR) permet de créer un fichier condensé pour
le tracé du diagramme HR \cfparagr{sec:fhr}; il ne contient que l'âge, $L$ et
$T_{\rm eff}$.
Ce fichier a pour extension {\tt .LR}, Ex: {\tt mon\_modele.LR}, \cfparagr{sec:flr} .
\end{description}
 \end{itemize}

\item{\tt N\_MAX=1000}~: le nombre maximum de couches est défini par
{\tt MAX(N\_MAX,nnnn)}, {\tt nnnn} est un nombre maximal de couches fixé dans
la routine {\tt cesam.f} suivant le niveau de précision requis;
le nombre minimal de couches est fixé à 500. {\bf A condition} que la
formation d'un
fichier ASCII {\bf soit requise} ({\tt NOM\_OUTPUT$\neq$no\_output}), une valeur
 négative de
{\tt N\_MAX} permet de calculer le dernier modèle de l'évolution avec un nombre
de couches égal à {\tt |N\_MAX|}.

\item{\tt PRECISION='pr'}~: niveau de précision requis. On trouve dans la
 \tab{tab:opt1}, les valeurs des paramètres numériques utilisées suivant
 les différentes options définies dans la routine {\tt cesam.f}.\index{precision}
Si le type de précision est codé en minuscules, Ex:{\tt 'sa'}, les variables
quasi-statiques utilisées sont: $\ln P,\ \ln T,\ (R/R_{\odot})^2,
\ (L/L_{\odot})^{2/3},\ (M/M_{\odot})^{2/3}$. S'il est codé
en majuscules, Ex:{\tt 'SA'}, elles seront $\ln P,\ \ln T,
\ R/R_{\odot},\ L/L_{\odot},\ (M/M_{\odot})^{1/3}$.\\
\begin{description}
\item[Suffixes A, T, L, M:]\index{suffixes de précision A, L \& M}
Tous les types de précision peuvent être affectés des suffixes {\tt A}
et/ou {\tt L} et/ou {\tt M}, Exemples: {\tt prA}, {\tt saLA}, {\tt SALA}
{\tt ceM}. L'ordre des suffixes est indifférent, Exemple: {\tt PRAL}.
%\cesam2k accepte qu'à l'occasion d'une reprise les suffixes
%apparaissent/disparaissent des {\tt NAMELISTS} l'évolution ultérieure se
%poursuivra avec/sans leur action.
%\begin{description}
\item [A]: ("ajuste") indique que le pas temporel sera ajusté
de façon à atteindre précisément les objectifs d'arrêt de la
NAMELIST {\tt NL\_EVOL} \cf\paragr{sec:nlevol}.
\item[L]: ("lisse") indique qu'un lissage local des discontinuités
de composition chimique sera effectué à la suite de chaque ajustement
temporel.
\item[M]: ("maximum") indique que le nombre de couches sera inchangé au long de
l'évolution. Il sera fixé au nombre maximal {\tt N\_MAX} indiqué dans la NAMELIST
 {\tt NL\_CESAM} \cf\paragr{sec:nlcesam}.
\end{description}
De plus des suffixes spécifiques peuvent affecter certains des types de précision
implantés dans \cesam2k \footnote{Les "/" indiquent
une alternative.}~:
\index{np}\index{pr}\index{lm}\index{sa}\index{sp}\index{rg}
\index{NP}\index{PR}\index{LM}\index{SA}\index{SP}\index{AV}\index{CE}\index{He}
\begin{itemize}
\item {\tt 'np/NP'} (précision normale) utile pour se donner l'idée générale
d'une évolution,
\item {\tt 'pr/PR'} (précision réaliste) pour une évolution sans une recherche
particulière de précision c'est le {\tt\bf "DEFAUT"},
\item {\tt 'sp/SP'} (super précision) et {\tt 'sa/SA'} (précision solaire) si on
recherche une précision extrême. Avec {\tt 'sa/SA'}, à condition que la formation
d'un fichier ASCII soit requise, les derniers modèles d'une
évolution sont calculés avec le nombre maximum de couches {\tt N\_MAX},
\item {\tt 'av/AV, ce/CE, He'} précisions pour des évolutions vers les stades avancés.
\index{cepheide}\index{mon\_modele.cephe}
{\tt 'ce/CE'} (Céphéïdes) identique à 'AV' avec un contrôle du pas temporel
dans la zone d'instabilité de façon à la traverser en 10 pas temporel environ,
et création d'un fichier {\tt mon\_modele.cephe}, \cfparagr{sec:cephe}.\\
\noi {\bf Suffixe O}: \index{suffixe de precision O}
En codant 'ceO/CEO' le fichier ASCII de chaque modèle de la bande
d'instabilité sera formé.
\item {\tt 'He'} (flash de l'hélium) identique à 'AV', le pas temporel
minimal pouvant descendre au niveau de l'heure. La luminosité pouvant devenir
négative dans le noyau la précision {\tt 'he'} n'est pas prévue.
\item {\tt 'rg'} (réglages) avec ce type de précision le choix des paramètres
numériques est confié à l'utilisateur/trice {\embf averti/e}. Ils sont
introduits
par l'intermédiaire d'un fichier de nom {\tt 'reglages'} ou, mieux
{\tt 'mon\_modele.rg'}, placé dans l'environnement du calcul. On trouvera
au \paragr{sec:reg} un exemple ainsi que la description de son contenu.
La précision {\tt 'rg'} accepte les extensions de precision.
\end{itemize}

\end{itemize}
Les constantes de répartition sont fixées suivant le type de
précision. Les valeurs les plus utilisées sont :
{\tt ctel=0, ctep=-1, ctem=15, cter=0, ctet=-1}. Pour
utiliser des paramètres différents \cf\paragr{sec:reg}.
\index{A} \index{L} \index{precision} \index{majuscules/minuscules}

\begin{table}\index{rg}\index{pr}\index{sp}\index{lm}\index{av}\index{co}\index{
sa}
\caption{\small Valeurs des paramètres numériques de la {\tt precision} par défaut
({\tt pr/PR}) et valeurs adaptées aux autres précisions.
%Pour ces paramètres, en codant {\tt precision='rg'},
%l'utilisateur peut utiliser ses propres valeurs en les codant
%dans le fichier de nom {\tt reglages} placé dans l'environnement du calcul,
%\cf \paragr{sec:reg}.
Les dénominations s'interprètent de la façon suivante~: {\tt pr/PR}~:
précision
réaliste ({\bf "le défaut"}), {\tt sp/SP}~: super précision,
{\tt lm/LM}~: pour petites masses (low mass), {\tt av/AV}~: pour modèles
pouvant
atteindre les stades avancés, {\tt np/NP}~: pour modèles
de précision normale, {\tt co/CO}~: pour modèles destinés à l'exploitation des
résultats du satellite CoRoT, {\tt sa/SA}~: pour modèles solaires.
Pour les Céphéides utiliser le type de précision {\tt ce/CE}, analogue à
{\tt av/AV}, et produisant des fichiers d'exploitation spécifiques.
La signification des différents paramètres numériques est explicitée au
\paragr{sec:reg} et lignes 805-1237 dans la routine
{\tt cesam.f} du sous-directory {\tt SOURCE}.
La valeur initiale du pas temporel et ses limites de variations ultérieures
dépendent de la masse et du type de précision requis; les valeurs portées
dans ce tableau n'ont qu'une valeur indicative. L'utilisation des
{\tt minucules/MAJUSCULES} (Exemple: {\tt pr/PR}) du type de précision définit
le type de variables quasi-statiques à utiliser, \cfparagr{sec:nlcesam}.
Une astérique* signale ces paramètres peuvent être modifiés lorsque
la masse initiale différe de l'unité ou que
la diffusion microscopique est prise en compte.}
\label{tab:opt1}
\begin{tt}
\begin{tabular}{cccccccc} \\ \hline \\
         &pr/PR    &sp/SP    &lm/LM   &av/AV    &np/NP   &co/CO    &sa/SA \\ \\
\hline \\
m\_qs    &2        &         &        &         &1 \\
m\_ch*   &2        &3        &        &         &        &3         &3\\
m\_rot   &2        &         &        &         &1       &3\\
m\_tds	 &2        \\
m\_ptm 	 &2        \\
ordre    &4	   &         &        &         &1  \\
precic   &$10^{-4}$\\
precix*  &$10^{-3}$&$10^{-4}$&	      &         &0.005  &$10^{-4}$&$10^{-5}$\\
precit*  &0.1      &0.02     &        &         &       &0.05     &0.02\\
psi0	 &0.03	   &0.02     &        &         &0.1    &         &0.02\\
loc\_zc  &$10^{-3}$&$10^{-4}$&        &         &0.005  &$10^{-4}$&$10^{-5}$\\
dt0*     &10       \\
dtmin*   &1an      &         &        &1sec \\
dtmax*   &100      &         &300     &         &300    &50       &50\\
d\_lum   &0.1      \\
d\_press &0.1      \\
d\_ray   &0.1      \\
d\_temp  &0.1      \\
age\_car*&1000     \\
ini0	 &3        \\
n\_atm   &75       &100      &        &         &50     &100      &100\\
kipp     &.TRUE.   &.FALSE.  &        &         &       &.FALSE.  &.FALSE. \\
ctel     &0        \\
ctep     &-1       \\
ctem     &15       \\
cter     &0.       \\
ctet     &-1.      \\
\\ \hline \\
\end{tabular}
\end{tt}
\end{table}

\begin{table}
\caption{\small Valeurs des paramètres numériques (suite).}
\begin{tt}
\begin{tabular}{cccccccc} \\ \hline \\
         &pr/PR    &sp/SP    &lm/LM   &av/AV    &np/NP  &co/CO    &sa/SA \\ \\
\hline \\
dn\_fixe &0.05	   \\
dpsi     &0.05	   \\
mu\_saha &.TRUE.   &         &.FALSE. &         &.FALSE.\\
ajuste*  &.FALSE.  \\
lisse*   &.FALSE.  \\
npt\_lisse*&10     \\
q0       &0.05     &0.01     &0.      &0.       &0.     &0.01     &0.01\\
l0       &0        &         &        &         &       &5        &5\\
new\_bv  &.TRUE.   &         &        &         &       &         &.FALSE.   \\
fcv      &.FALSE.  &         &.TRUE.  &         &.TRUE. \\
no\_discon &.FALSE.&         &        &.TRUE.   \\
n\_min\_ZC &5      \\
d\_zc    &10       &5        &100     &         &100     &2        &2\\
l\_fac   &.FALSE.  & .TRUE.  &        &         &        &.TRUE.   &.TRUE.\\
\\ \hline \\
\end{tabular}
\end{tt}
\end{table}

\section{Arguments de la NAMELIST {\tt NL\_MASS}}\label{sec:nlm}
\index{NL\_MASS}\index{mtot}\index{nom\_pertm}\index{mdot}
\begin{itemize}
 \item {\tt MTOT=1.d0}~: masse totale initiale, en $M_\odot$,
 \item {\tt NOM\_PERTM}~: nom du formalisme de perte de masse externe~:
 \index{nom\_pertm}\index{pertm\_ext.f}\index{pertm\_msol.f}\index{pertm\_tot.f}
 \index{pertm\_waldron.f}\index{pertm\_dominik.f}  \index{pertm\_reimer.f}
 \index{pertm\_reim-domi.f}
 \begin{itemize}\index{pertm\_waldron.f}
  \item {\tt pertm\_dominik}~: perte de masse selon l'article
  "stationary dust-driven winds", de Dominik \& al. 1990 A\&A, 240, 365
  \item {\tt pertm\_ext}~: perte/gain de masse externe
  linéaire en fonction du temps.
  \item {\tt pertm\_msol}~: perte/gain de masse externe linéaire
  en fonction du temps, la masse de l'étoile restant au moins
  égale à $1\,M_\odot$,
  \item {\tt no\_pertm}~: pas de perte/gain de masse externe,
  \item {\tt pertm\_reimers}~: perte de masse solaire empirique de Reimers
  cité dans l'article Rasio \& al.ApJ 470, 1187, 1996 (taux=0.6d-13 Msol/y)
  \item {\tt pertm\_reim-domi}~: association de {\tt pertm\_reimers} et
  {\tt pertm\_dominik},
  \item {\tt pertm\_waldron}~: perte de masse empirique de \citet{w85}.
 \end{itemize} \index{pertm\$}
 Si le nom du formalisme de perte de masse externe est affecté d'un '\$',
 Exemple: {\tt pertm\_msol\$}, il sera tenu compte de la perte de masse 'interne'
 due aux réactions thermonucléaires ($E=mc^2$), dont on tiendra compte avec
 {\tt no\_pertm\$}\footnote{Cette perte est de l'ordre de $10^{-14}M_\odot$/an
 pour le soleil.}.

\item {\tt MDOT=1.d-14}~: taux de perte/gain de masse, en $M_\odot$/an,
 ce taux est positif pour un gain de masse, négatif pour une perte de
 masse. Dans le cas standard, la composition chi\-mi\-que de l'apport/retrait
 de masse est celle de la couche la plus externe du modèle. Voir
 \paragr{sec:vent} comment personnaliser, si besoin, cette composition
 chi\-mi\-que.
 \end{itemize}

\section{Arguments de la NAMELIST {\tt NL\_EVOL}}\label{sec:nlevol}
\index{NL\_EVOL}\index{agemax}\index{arret}
\index{dtlist}\index{log\_teff}\index{R\_stop}
\index{nb\_max\_modeles}\index{t\_stop}\index{x\_stop}\index{he\_core}
\begin{itemize}
 \item {\tt AGEMAX=4.5d4}~: âge maximal à atteindre en million d'années,
  \item {\tt ARRET='post'}~: arrêt sur la ZAMS ({\tt'zams'}), à la fin de
  la ZAMS ({\tt'post'}\footnote{{\tt'tams'} est accepté.}),
  au début de la combustion de l'hélium ({\tt'cohe'}) ou du carbone
  ({\tt'coca'})
  ou du néon ({\tt'cone'}) ou de l'oxygène ({\tt'coox'}) ou autre
({\tt'else'}).
 \item {\tt DTLIST=1.d3}~: intervalle de temps minimum, en million d'années,
 séparant deux listes détaillées du modèle (fichier {\tt mon\_modele.lis}).
 Avec DTLIST$< 0$, Ex: {\tt DTLIST=-50.d0}, possibilité de construire des
 fichiers de sortie ASCII ou de reprise à intervalles réguliers ici
 tous les 50 modèles. Ces fichiers seront dénommés
 {\tt {\em nnnn}-mon\_modele-ad.osc}, et/ou
 {\tt {\em nnnn}-mon\_modele\_B.rep}; {\em nnnn} étant le numéro du modèle.
 Avec {\tt DTLIST=-1.d0} ces fichiers intermédiaires seront formés pour
 tous les modèles.
 \item {\tt LOG\_TEFF=-3.7d0}~: arrêt si cette valeur de $\log_{10}T_{\rm eff}$
 a été traversée, dans le sens croissant si {\tt log\_teff} est positif,
 dans le sens décroissant si {\tt log\_teff} est négatif,
 \item {\tt NB\_MAX\_MODELES=5000}~: arrêt après avoir calculé le nombre de
 modèles indiqué, éventuellement 0. A l'issue de chaque pas temporel le fichier
 binaire {\tt mon\_modele\_B.rep} remplace celui obtenu au pas temporel
 précédent. Le nombre maximal est 99999.
 Il est possible d'écrire et de conserver ces fichiers binaires
 à intervalles réguliers. Pour ce faire coder, par exemple,
 {\tt NB\_MAX\_MODELES=-5000} et {\tt DTLIST=-10.d0}, alors tous les
 10 modèles, le fichier binaire de reprise sera écrit et conservé dans
 l'environnement avec le nom {\tt{\em nnnn}-mon\_modele\_B.rep}
 Exemple: {\tt 0145-sun\_B.rep}.
 Le numéro 0000 est affecté aux modèles d'initialisation de \ZAMS\ ou de \PMS;
 en cas de reprise, les numéros reprennent à partir de celui du modèle repris.

 \index{numérotation des modèles}
 \item {\tt HE\_CORE=0.1d0}~: arrêt, dès que la masse du c\oe ur d'hélium
 atteint la valeur {\tt he\_core} ($M_\odot$),
 \item {\tt R\_STOP=0.9d0}~: arrêt, lorsque la valeur du rayon
 bolométrique traverse {\tt R\_STOP} en $R_\odot$,
 \item {\tt T\_STOP=1.d7}~: arrêt, si au centre, la température dépasse cette
 valeur,
 \item {\tt X\_STOP}=0.1d0~: arrêt si cette valeur de X, au centre, est
 traversée.
\end{itemize}
Si le suffixe {\tt A} affecte la variable de précision,
({\tt PRECISION='saA'})\index{suffixe de précision A}
\cfparagr{sec:nlcesam}, le pas temporel est ajusté de façon à obtenir
les valeurs requises.

\subsection{Caractérisation des types de modèles}
\index{zams}\index{post}\index{cohe}\index{coca}\index{coox}\index{cone}
\index{hom}\index{rep}\index{pms}\index{tams}
Au cours d'une évolution, à l'issue de chaque initialisation et de chaque pas
temporel, un fichier binaire contenant l'ensemble des variables du modèle est créé.
Ces fichiers
servent d'une part à poursuivre le calcul de l'évolution et d'autre part à
exploiter les modèles obtenus, calcul des fréquences d'oscillation par exemple.
A l'issue de chaque pas temporel le
fichier du nouveau modèle remplace celui du précédent. Les modèles qui
correspondent à des époques caractéristiques d'une évolution sont conservés, par
exemple à l'amorçage du 3$\alpha$. Ces fichiers dits {\em permanents}
permettent, en particulier, de reprendre l'évolution à ces époques caractéristiques.
Sous réserve d'acceptation par \cesam2k, il est alors possible de
modifier un ou plusieurs paramètres ou d'introduire des modifications de la physique,
par exemple la longueur de mélange ou
l'opacité; modifications qui seront prises en compte dans la suite de
l'évolution. Ces modèles binaires se différencient par leurs extensions qui
sont les suivantes:
\begin{itemize}
\item {\tt 'hom'} : Extension du nom d'un fichier de modèle de ZAMS homogène d'âge 0,
Exemple: {\tt mon\_modele.hom}, ces fichiers sont permanents.
\item {\tt 'pms'} : Extension du nom d'un fichier de modèle de PMS homogène d'âge 0,
Exemple: {\tt mon\_modele.pms}, ces fichiers sont permanents.
\item {\tt 'rep'} : Extension du nom d'un fichier d'un modèle en cours d'évolution,
Exemple: {\tt mon\_modele.rep}, Au fur et à mesure de l'évolution ces fichiers,
dits de "reprise", sont effacés et remplacés par celui du  suivant, ils correspondent  donc au dernier modèle calculé.
\item {\tt 'zams'} : Un modèle décrivant la pré-séquence principale est
identifié comme modèle de ZAMS dès que la génération d'énergie d'origine
nucléaire surpasse celle d'origine graviphique. A l'époque de la transition le fichier de reprise est copié dans un fichier
permanent d'extension
{\tt zams}, Exemple: {\tt mon\_modele.zams}.
\item {\tt 'post'} ou {\tt 'tams'} : Un modèle décrivant la séquence principale est
identifié comme modèle de post-main séquence dès que l'abondance centrale
d'hydrogène devient inférieure à 0.01,
A l'époque de la transition le fichier de reprise est copié dans un fichier permanent
d'extension {\tt 'post'} ou {\tt 'tams'}.
Exemple: {\tt mon\_modele.post} ou {\tt mon\_modele.tams}.
\item {\tt 'cohe'} : Un modèle décrivant la post-séquence principale est
identifié comme modèle avec combustion d'hélium dès que la température
centrale dépasse $1.\,10^8$K. A l'époque de la transition le fichier de
reprise est copié dans un fichier permanent d'extension {\tt 'cohe'},
Exemple: {\tt mon\_modele.cohe}
\item {\tt 'coca'} : Un modèle décrivant la séquence de combustion de
l'hélium est identifié comme modèle avec combustion du carbone dès que
la température centrale dépasse $7.\,10^8$K. A l'époque de la transition le fichier de reprise est copié dans un fichier permanent
d'extension {\tt 'coca'}.
\item {\tt 'cone'} : Un modèle décrivant la séquence de combustion de
carbone est identifié comme modèle avec combustion du néon dès que
la température centrale dépasse $1.\,10^9$K. A l'époque de la transition le fichier de reprise est copié dans un fichier permanent
d'extension {\tt 'cone'}.
\item {\tt 'coox'} : Un modèle décrivant la séquence de combustion du
néon est identifié comme modèle avec combustion d'oxygène dès que
la température centrale dépasse $1.9\,10^9$K. A l'époque de la transition le fichier de reprise est copié dans un fichier permanent
d'extension {\tt 'coox'}.
\end{itemize}

\section{Arguments de la NAMELIST {\tt NL\_CHIM}}\label{sec:nlchim}
\index{zsx0}\index{y0}\index{x0}\index{NL\_CHIM}\index{grille\_fixe}
\index{garde\_xish}\index{modif\_mix}\index{modif\_chim}
\begin{itemize}
 \item {\tt GRILLE\_FIXE=.TRUE.}~: pour utiliser une grille fixe pour la
composition chi\-mi\-que et la diffusion du moment cinétique\footnote{Inactivée
dans la présente version de \cesam2k.},
 \item {\tt NOM\_ABON}~: nom de la mixture
  initiale\footnote{\cf{\EOS\ OPAL: {\sc web} site
 \tt http://www-phys.llnl.gov/V\_Div/OPAL/opal.html}.}~:
 \index{nom\_abon}\index{abondances initiales}
 \index{enhan\_al}\index{enhan\_cha}\index{enhan\_w}\index{météorites}
 \index{solaire\_gn}\index{solaire\_gs}\index{solaire\_ags\_03}
 \index{solaire\_ags\_05}
 \begin{itemize}
  \item{\tt enhan\_al}~: abondances $\alpha$-enhanced de Allard,
  \item{\tt enhan\_cha}~: abondances $\alpha$-enhanced de Chaboyer,
  \item{\tt enhan\_w}~:  abondances $\alpha$-enhanced de Weiss,
  \item{\tt meteorites\_ag}~: abondances météoritiques de \citet{ag89},
  \item{\tt meteorites\_gs}~: abondances météoritiques de \citet{gs98},
  \item{\tt meteorites\_ags05}~: abondances météoritiques de \citet{ags05},
  \item{\tt meteorites\_agss}~: abondances météoritiques de \citet{agss}
  \item{\tt mixture}~: forme à personnaliser de la compilation d'abondances
  décrite au \paragr{sec:mix},
  \item{\tt solaire\_gn}~: mixture solaire de \citet{gn93},
  \item{\tt solaire\_gs}~: mixture solaire de \citet{gs98}.
  \item{\tt solaire\_ags\_03}~: mixture solaire de \citet{a03}
  \item{\tt solaire\_ags\_05}~: mixture solaire de \citet{ags05}
  \item{\tt solaire\_agss}~: abondances météoritiques de \citet{agss}.
 \end{itemize}
 Ces abondances sont initialisées dans la routine {\tt abon\_ini.f}.
 \item {\tt MODIF\_CHIM=.FALSE.}~: signifie que \cesam2k devra utiliser, sans
 modification, les données du sous programme {\tt abon\_ini.f}. Si tel n'est pas
 le cas, coder {\tt MODIF\_CHIM=.TRUE.}, \cesam2k recherchera alors dans
 l'environnement les fichiers permettant les modifications souhaitées.
 Ces fichiers de modification sont:
 \begin{itemize}
 \item{modif\_mix} : composition chi\-mi\-que initiale modifiée, \cfparagr{sec:rapab},
 \item{rap\_iso} : rapports isotopiques modifiés, \cfparagr{sec:rapiso},
 \item{planet} : composition chimique des planétoïdes différente de celle du modèle,
 \cfparagr{sec:cplanet},
 \item{vent} : composition chimique du vent différente de celle du modèle,
 \cfparagr{sec:vent}.
 \end{itemize}
 \cesam2k invitera à supprimer ces fichiers avant d'effectuer des
 calculs qui n'en tiendraient pas compte. Cette disposition
 {\embf à but sécuritaire} a pour fin
 d'éviter de prendre en compte des fichiers qui, par mégarde,
 n'auraient pas été soustraits de l'environnement.
 Si aucun de ces fichiers n'est détecté, \cesam2k le signalera et l'instruction
 {\tt MODIF\_CHIM=.TRUE.} sera ignorée.

\item {\tt GARDE\_XISH=.FALSE.}~: ce paramètre fait l'objet du paragraphe
 suivant,
 \item {\tt X0=0.7}~: abondance initiale, par unité de masse, de H,
 \item {\tt Y0=0.28}~: abondance initiale, par unité de masse, de He,
 \item {\tt ZSX0=0.d0}~: valeur initiale de Z/X,
  \begin{itemize}
  \item Les abondances du modèle initial sont déterminées à partir
  des quantités {\tt X0}, {\tt Y0} et {\tt ZSX0}; en raison de la relation
  $1=X+Y+Z$, seulement 2 de ces 3 quantités sont à prendre en
  considération. La valeur initiale de Y utilisée dans le calcul
  est {\embf toujours} {\tt Y0}; elle est fractionnée en les divers isotopes
  retenus pour décrire l'évolution de l'hélium,
  \item si {\tt ZSX0 $\leq$ 0},
   l'abondance initiale d'hydrogène est initialisée à {\tt X0} et $Z=1-X-Y$
   est déduit de {\tt X0} et {\tt Y0},
  \item sinon, \ie\ {\tt ZSX0 $>$ 0}, l'abondance initiale d'hydrogène
  est déduite de {\tt Y0} et de {\tt ZSX0},
   la valeur indiquée pour {\tt X0} étant ignorée.
   \end{itemize}
 \end{itemize}

 \noi Les valeurs globales initiales {\tt X0}, {\tt Y0} et {\tt Z0} sont fractionnées
 entre les divers éléments et chacun d'eux entre leurs isotopes.

\noi Avec diffusion et/ou après un premier dredge-up, la composition
chi\-mi\-que
des couches extérieures, \ie\ celle observée, diffère de la composition
chi\-mi\-que initiale. Il y a lieu d'ajuster la valeur initiale de
la composition chi\-mi\-que, afin d'obtenir celle observée à l'issue d'une
évolution. Les observations portant sur les rapports métal/H, ou encore
[Fe/H], à l'initiative de F.Thévenin, \cesam2k offre divers outils
pour réaliser cet ajustement~:
\begin{itemize}
 \item Conserver le rapport $Z/X$ de la mixture initiale, au lieu de
 {\tt ZSX0}, objet du \paragr{sec:msh},
 \item utiliser une mixture originale, \cf\paragr{sec:mix},
 \item modifier les rapports isotopiques, \cf\paragr{sec:rapiso},
 \item modifier les rapports d'abondances de la mixture initiale,
  \cf\paragr{sec:rapab}.
 \end{itemize}

\subsection{Alternative~: conservation métal/H ou métal/Z}
\index{métal/H}\index{métal/Z}\index{calibration solaire}
\index{garde\_xish}\label{sec:msh}
\noi Dans \cesam2k, l'utilisation de tabulations est incontournable pour le calcul de l'opacité et de l'équation l'équation d'état.
Dans les tables les valeurs des paramètres
globaux  de composition chimique $X$, $Y$, $Z$ sont fixés.
Il y a des incohérences chroniques entre les valeurs {\em fixes} utilisées pour les
tabulations et celles {\em variables} obtenus au cours de l'évolution. Le rôle de
{\tt GARDE\_XISH} est de préciser la nature de ces incohérences que l'on doit
accepter, faute de mieux.\medskip

\noi On désigne par {\em mixture de base} la mixture proposée par \cesam2k,
Exemple {\tt solaire\_gn}, ou celle fournie par l'utilisateur \cfparagr{sec:mix}.
Elle inclut éventuellement les aménagements décrits aux
\paragr{sec:rapab}, ou encore \paragr{sec:rapiso}.\par
\noindent La {\em mixture initiale} est celle utilisée pour
initialiser une évolution. Elle
est déduite de la mixture de base en tenant compte des indications codées
dans le fichier de données.\medskip

\noi La valeur initiale de $Y$ est {\embf TOUJOURS} la valeur lue
dans le fichier de données, \ie\ {\tt Y0}, par exemple la valeur primordiale
{\tt Y0=0.27}.

\medskip

\noi {\tt GARDE\_XISH=.TRUE.} \respt{{\tt GARDE\_XISH=.FALSE.}}
permet de fixer la valeur initiale de $Z$ \respt{$X$} de façon à conserver
le rapport {\tt métal/H} \respt{{\tt métal/Z}}
\footnote{Alternative réalisée dans la routine {\tt abon\_ini.f}.}.
\par

\noi Ainsi, en codant\footnote{De meilleures dénominations seraient, par
exemple~:
{\tt GARDE\_XISH\_MIX}, {\tt GARDE\_XISH\_MIX\_INI}, {\tt GARDE\_ZSX}
 ou encore {\tt GARDE\_ZSX\_MIX\_INI},\etc}~:
\begin{itemize}
\item  {\tt GARDE\_XISH=.TRUE.},
la valeur de $Z/X$ de la mixture initiale sera déduite
des rapports d'abondances métal/H de la mixture de base et non pas de celle,
{\tt ZSX0}, lue dans le fichier de données.
La métallicité utilisée $Z/X$, {\bf sera différente} de celle ZSX0 indiquée dans
le fichier de données, sauf dans le cas particulier
où les valeurs de Y0 et de ZSX0 correspondent
{\bf exactement} à celles de la mixture de base.
Codage qu'il convient d'utiliser quand on désire que la mixture initiale et la
mixture de base soient "{\em physiquement}" identiques.
\item {\tt GARDE\_XISH=.FALSE.}, la valeur {\tt ZSX0}
lue dans le fichier de données, sera utilisée. Les rapports métal/Z des éléments
chimiques seront ceux de la mixture de base. Codage qu'il convient d'utiliser
pour une {\embf calibration solaire}, pour laquelle on vise à atteindre une
valeur fixée de $Z/X$, en utilisant, pour l'initialisation,le Y0 du fichier de données
et les rapports
métal/Z de la mixture de base qui peuvent éventuellement différer du $Z/X$ à atteindre.
\end{itemize}

\medskip

\PB{L'utilisation de {\tt GARDE\_XISH=.TRUE.} {\embf EST DELICATE}, il est
conseillé de bien vérifier, dans le fichier {\tt mon\_modele.lis},
que les valeurs initiales des abondances utilisées par \cesam2k
correspondent bien à ce que l'on désire.}

\medskip
\PB{Il faut être conscient que les modifications de la mixture ne peuvent
être prises en compte, ni dans les opacités, ni dans l'équation d'état.}

\section{Arguments de la NAMELIST {\tt NL\_CONV}}\label{sec:nlconv}
\index{alpha}\index{cpturb}\index{ovshts}\index{ovshti}\index{jpz}
\index{Ledoux}\index{NL\_CONV}
\begin{itemize}
 \item {\tt NOM\_CONV}~: nom de la routine calculant le gradient de
 température dans les zones convectives.



 \begin{itemize}
  \item{\tt conv\_a0}~: convection MLT, avec longueur de mélange
  $l\rightarrow 0$ aux limites ZR/ZC,
  \item{\tt conv\_cgm\_reza}~: convection suivant \citet{cm91}, avec
  $l=\alpha H_{\rm p}$, tient compte de la prescription de Bernkopf,
  \item{\tt conv\_cm}~: convection suivant \citet{cm91}, avec
  $l=\alpha H_{\rm p}$,
  \item{\tt conv\_cml}~: convection suivant Canuto \& Mazitelli avec
  avec longueur de mélange égale à la plus courte distance des
  bords de la ZC,
  \item{\tt conv\_cm\_reza}~: convection suivant \citet{cm91}, avec
  $l=\alpha H_{\rm p}$, tient compte de $\delta$,
  \item{\tt conv\_jmj}~: convection MLT avec $l=\alpha H_{\rm p}$ \citep{b58}.
  \item {\bf Suffixe S}: Pour les modèles utilisant les types de précision "corot"
  ({\tt co, CO}), "super précision" ({\tt sp,SP}), "précision solaire"
  ({\tt sa, SA}), \cesam2k place chaque limite ZR/ZC à une distance d'un
  point de grille au plus égale à une valeur fixée. Pour éviter cette disposition
  qui ralentit le déroulement du calcul, accoler le suffixe {\bf S} au nom de la
  routine de convection, Exemple: {\tt conv\_jmjS}.\index{suffixe de convection S}
 \end{itemize}

 \item {\tt alpha=1.7d0}~: longueur de mélange,
 \item {\tt ovshts=0.05d0}~: coefficient d'overshoot supérieur, avec une valeur
 positive il y aura une {\bf discontinuité} du gradient pris égal au gradient
 adiabatique
 dans la partie overshootée, avec une valeur négative il y a {\bf continuité} du
 gradient pris égal égal au gradient radiatif dans la partie overshootée,
 \item {\tt ovshti=0.1d0}~: coefficient d'overshoot inférieur, comme pour
 l'overshoot supérieur, dans la partie overshootée une valeur
 positive affecte le gradient adiabatique,
 une valeur négative affecte le gradient radiatif,
 \item {\tt jpz=.TRUE.}~: utilisation des prescriptions de JpZ,
  \item {\tt cpturb=0.d0}~: coefficient de pression turbulente,
  {\bf indisponible dans cette version} de \cesam2k,
 \item {\tt ledoux=.TRUE.}~: utilisation du critère de
 Ledoux. Avec semi-convection on doit coder {\tt ledoux=.FALSE.}
 le formalisme de Ledoux faisant partie intégrante de la description de la
 semi-convection, \cfparagr{sec:nldiff}.
\end{itemize}

\medskip

%\PB{L'utilisation d'un paramètre de pression turbulente non nul est délicate.}

\section{Arguments de la NAMELIST {\tt NL\_DIFF}}\label{sec:nldiff}
 \index{diffusion mi\-cros\-co\-pi\-que}\index{diffm\_br.f}\index{diffm\_mp.f}
 \index{diffm\_0}\index{nom\_diffm}\index{diffusion}\index{difft\_gab.f}
 \index{difft\_nu.f}\index{d\_turb}\index{re\_nu}\index{nom\_frad}
 \index{alecian1.f}\index{alecian2.f}\index{no\_frad}\index{difft\_smc.f}
\index{diffusion turbulente}\index{difft\_nut.f}\index{nom\_difft}
\index{semi\_convection}\index{suffix V}
 \begin{itemize}
 \item {\tt DIFFUSION=.TRUE.}~: on tiendra compte de la diffusion
 mi\-cros\-co\-pi\-que des éléments chi\-mi\-ques.
 %, éventuellement de la diffusion du moment cinétique,
 \item {\tt NOM\_DIFFM}~: nom de la routine calculant les coefficients
 de diffusion mi\-cros\-co\-pi\-que~:
 \begin{itemize}
  \item{\tt diffm\_br}~: calcul des coefficients de diffusion
  mi\-cros\-co\-pi\-que suivant le formalisme de Burgers en tenant compte,
  éventuellement, des accélérations radiatives,
  \item{\tt diffm\_mp}~: calcul des coefficients de diffusion
  mi\-cros\-co\-pi\-que
  suivant le formalisme simplifié de Michaud \& Proffit,
  \item{\tt diffm\_0}~: coefficients de diffusion mi\-cros\-co\-pi\-que nuls.
 \end{itemize}
 \item {\tt NOM\_DIFFT}~: nom de la routine calculant les coefficients
 de diffusion turbulente~:
 \begin{itemize}
  \item{\tt difft\_nu}~: calcul des coefficients de diffusion turbulente
  incluant la diffusivité radiative.
   \item{\tt difft\_nut}~: calcul des coefficients de diffusion turbulente
  incluant la diffusivité radiative et la diffusivité turbulente dans la
  tachocline de la \ZC\ externe suivant \citet{cvr07}.
 \item{\tt difft\_gab}~: calcul des coefficients de diffusion turbulente
  selon une pres\-crip\-tion simple de M.Gabriel.
 \item{\tt difft\_smc}~: calcul des coefficients de diffusion turbulente
  incluant la semi-convection et la diffusivité radiative.
\item{\tt difft\_sun}~: coefficients de diffusion turbulente, sous la \ZC
 solaire suivant \cite{g97}.
 \end{itemize}

 \noi{\bf Suffixe V}:\index{suffixe de diffusion V}
 Quand le suffixe "V" est accolé à la fin du nom de la routine de diffusion
 turbulente à utiliser, Exemple : difft\_nuV, le coefficient de diffusion
 permettant le mélange
 convectif est pris égal à $1/3 V_{\rm conv} L_{\rm mel}$ .

\item {\tt D\_TURB=1.d1}~: coefficient de diffusion turbulente isotrope,
 \item {\tt RE\_NU=1.d0}~: coefficient de diffusivité radiative,
 \item {\tt NOM\_FRAD}~: nom de la routine de calcul des accélérations
 radiatives,
 \begin{itemize}\index{accélérations radiatives} \label{sec:accrad}
  \item{\tt frad\_test}~: utilisation d'accélérations radiatives pour tests,
  \item{\tt alecian\_bp}~: calcul des accélérations radiatives suivant le
  second for\-ma\-lisme de G.Alécian\footnote{Les données peuvent être
  obtenues sur sur simple demande auprès de Bernard.Pichon@oca.eu},
  \item{\tt no\_frad}~: les accélérations radiatives seront ignorées.
 \end{itemize}

% \PB{Les accélérations radiatives ne sont pas disponibles dans la présente
% version de \cesam2k.}
\end{itemize}

\section{Arguments de la NAMELIST {\tt NL\_ROT}}\label{sec:rot}
\index{w\_rot}\index{dws}\index{NL\_ROT}
\PB{\bf La diffusion du moment cinétique n'est pas disponible dans
 la présente version de \cesam2k.}
\begin{itemize}
 \item{\tt W\_ROT=50.d0}~: vitesse angulaire initiale.
 \item{\tt UNIT}~: unité utilisée pour la vitesse angulaire initiale.
 différents choix sont proposés de façon à
 s'adapter facilement aux observations~:
 \begin{itemize}\index{unit}\index{kms/s}\index{jours}\index{rad/s}
  \item{\tt jours}~: période initiale en jours, correspond à la vitesse
  angulaire initiale de rotation de la couche externe,
  \item{\tt kms/s}~: vitesse de la couche externe en kilomètres/seconde,
  le rayon du modèle d'initialisation étant utilisé pour déterminer la vitesse
  angulaire initiale, il convient d'ajuster ce rayon en reprenant plusieurs fois
  le calcul du modèle d'âge 0,
  \item{\tt rad/s}~: la vitesse angulaire initiale est en radians/seconde.
 \end{itemize}
 \index{diffw\_p03.f}\index{diffw\_mpz.f}\index{diffw\_0.f}
  \item{\tt NOM\_DIFFW}~: nom de la routine du calcul des coefficients de
  diffusion du moment cinétique~:
  \begin{itemize}
   \item{\tt diffw\_p03}~: coefficients de
  diffusion du moment cinétique suivant \citet{p03}
   \item{\tt diffw\_mpz}~: coefficients de
  diffusion du moment cinétique suivant \cite{mpz04}
   \item{\tt diffw\_0}~: coefficients de diffusion du moment cinétique nuls.
   \item{\tt diffw\_cte}~: coefficients de diffusion du moment cinétique
   constants~: \\ ($D_{\rm eff}=300,\ D_{\rm h}=1\,000\,000,\ D_{\rm v}=250$).
   \item{\tt diffw\_toul}~: coefficients de diffusion  du moment cinétique
    suivant \cite{cvr07} \index{diffw\_toul.f}
  \end{itemize}
  \item{\tt NOM\_THW}~: Nom de la théorie utilisée pour l'évolution de
  la vitesse angulaire~:
  \begin{itemize}\index{rot\_0}\index{nom\_thw}\index{cons\_glob\_mnt\_cin}
   \item{\tt rot\_0}~: Vitesse angulaire nulle, le modèle ne tient pas compte
   de la rotation.
   \item{\tt rot\_cte}~: Vitesse angulaire constante au cours de l'évolution,
   \ie rotation rigide.\index{rot\_cte}
   \item{\tt cons\_glob\_mnt\_cin}~: Evolution avec rotation rigide, et
   conservation globale du moment cinétique. A l'issue de chaque pas temporel,
   la vitesse angulaire, constante dans tout le modèle, est réajustée de façon à
   conserver le moment cinétique total.\index{diff\_tz97}\index{diff\_mz04}
   \item{\tt diff\_tz97}~: Evolution avec diffusion du moment cinétique, suivant
   la théorie de \citet{tzmm97}, {\bf indisponible dans cette version} de
\cesam2k.
   \item{\tt diff\_mz04}~: Evolution avec diffusion du moment cinétique, suivant
   la théorie de \citet{mz04}, {\bf indisponible dans cette version} de
\cesam2k.
   \item{\tt cons\_loc\_mnt\_cin}~: évolution avec conservation locale du moment
   cinétique avec rotation rigide des \ZCs.
  \end{itemize} \index{cons\_loc\_mnt\_cin}
 \item{\tt NOM\_PERTW}~: Nom de la routine de calcul de perte de moment
 cinétique,
  \begin{itemize}
   \item{\tt pertw\_sch}~: variation de vitesse angulaire proportionnelle à
   $\Omega^3$ (Schumanish), {\bf indisponible dans cette version} de \cesam2k,
   \item{\tt pertw\_loc}~: perte de moment cinétique proportionnelle à
   l'énergie cinétique de rotation locale, {\bf indisponible dans cette version}
    de \cesam2k,
   \item{\tt pertw\_ptm}~: perte de moment cinétique, conséquence d'une
   perte de masse, {\bf indisponible dans cette version} de \cesam2k,
   \item{\tt pertw\_0}~: pas de perte de moment cinétique\footnote{Seule routine
  de perte de moment cinétique disponible dans la présente version
  de \cesam2k.}.
  \end{itemize}
 \item{\tt P\_PERTW=-1.d-9}~: paramètre de perte de moment cinétique attaché
 respectivement à chacune des routines précédentes
 \item{\tt  NOM\_DES\_ROT='end\_mod'}~: Contrôle de la génération de fichiers
  ASCII des variables de la diffusion du moment cinétique, \cfparagr{sec:frot}.
  Si le dessin on-line n'est pas requis ({\tt NOM\_DES=no\_des}),
  le fichier postscript des dessins des variables de la rotation est
  néansmoins automatiquement créé, {\bf indisponible dans cette version} de
   \cesam2k.
  \index{nom\_des\_rot}
 \begin{itemize}
 \item{\tt no\_des}~: Pas de formation de fichier.
 \item{\tt end\_evol}~: Formation du fichier à l'issue de l'évolution, le nom
du
 fichier est {\tt mon\_modele\_coeff\_rota.dat}.
 \item{\tt all\_mod}~: Formation du fichier à l'issue de chaque pas temporel,
  le nom du fichier est {\tt {\em nnnn}-mon\_modele\_coeff\_rota.dat},
 {\em nnnn} étant le numéro du modèle. Tous les fichiers ainsi créés sont
 permanents.
 \item{\tt end\_mod}~: Formation du fichier à l'issue de chaque pas temporel,
  le nom du fichier est {\tt mon\_modele\_coeff\_rota.dat}, seul le dernier
  fichier créé étant disponible.
 \item{\tt all\_iter}~: Formation du fichier à l'issue de chaque itération de
la
 résolution du système d'équations de la diffusion du moment cinétique;
 fichiers destinés aux mises au  point.\index{end\_evol}
 \index{all\_iter}\index{no\_des}\index{end\_mod}\index{all\_mod}
 \end{itemize}
\end{itemize}
\index{p\_pertw}\index{pertw\_0}\index{pertw\_ptm.f}
\index{pertw\_loc.f}\index{pertw\_sch.f}\index{nom\_pertw}
\index{nom\_diffw}

\section{Arguments de la NAMELIST {\tt NL\_ETAT}}\label{sec:nletat}
\index{f\_eos}\index{NL\_ETAT}
\begin{itemize}
\item {\tt NOM\_ETAT}~: nom de la routine d'équation d'état~:
 \begin{itemize}
  \item{\tt etat\_ceff}~: équation d'état de \cite{eff73}, avec
  corrections Coulombiennes ({\em appel à {\tt etat\_eff} en cas de
  difficulté}).
  \item{\tt etat\_eff}~: équation d'état de \cite{eff73},({\em appel
  à {\tt etat\_gong2} en cas de difficulté}).
  \item{\tt etat\_gong1}~: équation d'état de {\tt gong1} selon \cite{c88},
  seuls H et He sont pris en compte et supposés totalement ionisés.
  \item{\tt etat\_gong2}~: équation d'état de {\tt gong2} selon \cite{c88},
  seuls H et He4 sont pris en compte, dégénérescence ignorée.
  \item{\tt etat\_mhd}~: équation d'état de Mihalas, Hummer \& Dappen
  ({\em appel à {\tt etat\_eff} en cas de difficulté}) utilise les tables
  en binaire {\tt mhd1.bin} ... {\tt mhd7.bin} construites dans le
  sous-directory {\tt SUN\_STAR\_DATA}, voir \paragr{sec:mhdk}.

 \item{\tt etat\_irwin}~: équation d'état {\tt FreeEOS} de A. Irwin.
  La source en est
  téléchargeable sur {\tt http://freeeos.sourceforge.net}. Les extensions
  {\tt 1, 2, 3, 4, 5, 6} Ex: {\tt etat\_irwin1a}, donnent accès à différents
  niveaux de précision/rapidité et de la cohérence avec d'autres source d'EOS
  et/ou d'Opacité
  existantes. Une possibilité de personnalisation est offerte par
  l'intermédiaire
  d'un fichier de nom {\tt opteos} disposé dans l'environnement du calcul. Ces
  options sont décrites dans la routine {\tt etat\_irwin.f} et dans le fichier
  (texte) {\tt etat\_irwin\_explik.f} du sous-directory {\tt SOURCE}. L'option
  1: {\tt etat\_irwin1} est la plus précise et la plus lente, l'option 5:
  {\tt etat\_irwin5}, est la plus rapide et la moins précise. Cette équation
  d'état
  utilise la bibliothèque dont la créaton est décrite \paragr{sec:irwin}.
  {\em Ce sont les directives de Y.Lebreton qui ont permis l'implantation
  de cette EOS dans \cesam2k.}
  \index{etat\_irwin.f}\index{etat\_irwin\_explik}

 \item{\tt etat\_opal}~: équation d'état {\tt opal} 1991,
   ({\em  appel à
   {\tt etat\_eff} en cas de difficulté}) utilise les tables
   en binaire {\tt eos\_opal*.bin} construites dans le sous-directory
   {\tt SUN\_STAR\_DATA}, voir \paragr{sec:eosopal}.
   \item{\tt etat\_opalZ}~: équation d'état de {\tt opal} 2001,
   ({\em appel à {\tt etat\_eff} en cas de difficulté}) utilise les
   tables en binaire {\tt eos\_opalZ*.bin} construites pour {\bf Z fixé}
   dans le sous-directory {\tt SUN\_STAR\_DATA},
   \cfparagr{sec:eosopal}. Par rapport à la version 1991, les
   données de la version 2001 ont été étendues et améliorées, des bugs
   ont été corrigés. Les données ne permettent pas de dépasser $X=0.8$.
   Les fichiers de données 2001 et 1991 sont incompatibles.
   \item{\tt etat\_opalX}~: équation d'état de {\tt opal} 2001, se
   différenciant de {\tt etat\_opalZ} par l'extension du domaine de
   l'abondance d'hydrogène, $X\in [0,1]$, dans le but de calculer
   des modèles sans zone convective
   externe étendue {\em avec diffusion}. Le prix à payer est $Z\equiv 0$,
   ce qui n'est pas très pénalisant, l'équation d'état étant peu sensible à
   $Z$. {\tt etat\_opalX} utilise {\bf la}
   table en binaire {\tt eos\_opalX} construite pour {\bf Z=0 fixé}
   dans le sous-directory {\tt SUN\_STAR\_DATA},
   \cfparagr{sec:eosopal}.

   \item{\tt etat\_opal5}~: Equation d'état similaire aux précédentes.
   La construction de la table, pour une valeur de
   métallicité Z donnée, est effectuée dans le sous directory
   {\tt OPAL2005\_tables.gz}, \cfparagr{sec:eosopal5}.
   Dans le fichier de données, il ne sera pas nécessaire d'indiquer que la table de nom
   {\tt op5\_0174.bin} (ici $Z=0.0174$) est ou non compressée.\index{etat\_opal5}

   \item{\tt etat\_saha}~: Equation d'état tabulée à Z variable
   communiquée et adaptée à \cesam2k par
   V.A. Baturin \&  A.V. Oreshina du Sternberg Astronomical Institute, Lomonosov
   Moscow State. \cfparagr{sec:eosaha} \index{etat\_saha}

 \end{itemize}
  \item {\tt F\_EOS='eos\_opal\_250.bin',7*' '}~: noms des fichiers
  d'équation d'état (au plus 8). Ne {\bf pas} indiquer le nom des fichiers pour l'EOS
  {\tt etat\_saha}.
\end{itemize}
 \index{etat\_eff.f}\index{etat\_gong1.f}\index{etat\_gong2.f}
 \index{etat\_mhd.f}
 \index{etat\_opal.f}\index{etat\_opalZ.f}\index{etat\_opalX.f}
 \index{equation d'état}\index{nom\_etat}

\section{Arguments de la NAMELIST {\tt NL\_OPA}}\label{sec:nlopa}
 \index{f\_opa}\index{NL\_OPA}\index{opa\_int\_zsx.f}\index{opacité conductive}
 \index{opa\_yveline\_jorgen.f}\index{opa\_cond\_iben}\index{opa\_cond\_mestel}
 \index{opa\_opal2\_co}\index{opa\_opalCO.f}\index{opa\_opal2\_cno}\index{H12}
 \index{opacité}\index{opa\_gong.f}\index{opa\_houdek12.f}\index{opa\_yveline.f}
 \index{opa\_yveline\_lisse.f}\index{nom\_opa}\label{sec:nopa}
 \index{nom\_opa\_cond}
 \begin{itemize}
 \item {\tt NOM\_OPA}~: nom de la routine de calcul des opacités moyennes de
 Rosseland~:
 \begin{itemize}
  \item{\tt opa\_gong}~: opacités simplifiées (Kramers amélioré)
  selon \cite{c88}.
  \item{\tt opa\_houdek12}~: opacités de Houdek version 12,
  interpolation par rational B-spline; utilise les
  tables en binaire construites dans le sous-directory
  {\tt SUN\_STAR\_DATA/H12}, voir \paragr{sec:opahd}. Les opacités Houdek
  contiennent {\sl ad initio} des opacités conductives.
  \item{\tt opa\_int\_zsx}~: Opacités OPAL93+Kurucz raccord d'Yveline,
  interpolations linéaires, très utiles pour des tests.
  \item{\tt opa\_opal2\_co, opa\_opal2\_cno}~: Opacités OPAL avec corrections
  pour $Z>0.1$, \cf \paragr{sec:opal2}
  \item{\tt opa\_opalCO}~: Opacités OPAL avec corrections pour $Z>0.1$ ne
  portant que sur C et O.
  \item{\tt opa\_yveline}~: opacités OPAL+Alexander interpolation et
  raccord d'Yveline; utilise les tables
  en binaire {\tt opa\_yveline*.bin} construites dans le sous directory
  {\tt SUN\_STAR\_DATA}, voir \paragr{sec:opayv}.
  \item{\tt opa\_yveline\_jorgen}~: identique à la précédente pour les opacités
  \citet{ags05} avec la correction de Jorgen (astro-ph 0811.100v1 6 Nov 2008)
  qui permet de retrouver le "modèle S".
  \item{\tt opa\_yveline\_lisse}~: opacités OPAL+Alexander raccord
  d'Yveline, interpolation linéaire ou lissage; utilise les tables
  en binaire {\tt opa\_yveline*.bin} construites dans le sous-directory
  {\tt SUN\_STAR\_DATA}, voir \paragr{sec:opayv}.
 \end{itemize}
 \item {\tt F\_OPA='opa\_yveline.bin','/H12/OPINTPATH',6*' '}~:
  noms des fichiers d'opacité (au plus 8). Pour les opacités Houdek, voir les
  commentaires au début de la routine {\tt opa\_houdek12} de la SOURCE,
 en cas de difficulté les opacités {\tt opa\_yveline} étant utilisées par défaut.
\item {\tt NOM\_OPA\_COND}~: nom de la routine de calcul de l'opacité conductive.
 \begin{itemize}
  \item{\tt opa\_cond\_iben}~: opacité conductives de \citet{i75}
  \item{\tt opa\_cond\_mestel}~: opacité conductives de Mestel \citep{cg68}
  \item{' '}~: Opacité conductive ignorée.
\end{itemize}
\end{itemize}
Lorsque la métallicité $Z > 0.1$, Ex: au centre lors de la combustion de
l'hélium, la matière est supposée totalement
ionisée, \cesam2k simplifie le calcul de l'opacité en utilisant
le formalisme free-free de Kramers, \citet{kw91} (17.5-17.6),
avec la correction de Compton, \citet{cg68} (§16.a, table 16.1).

\section{Arguments de la NAMELIST {\tt NL\_NUC}}\label{sec:nlnuc}
\index{mitler}\index{NL\_NUC}\index{nom\_nuc\_cpl}\index{pp3.f}
\index{pp1.f}\index{nom\_nuc}\index{ppcno9Fe.f}\index{ppcno10K.f}
 \index{ppcno9.f}\index{ppcno10.f}\index{ppcno10Fe.f}\index{ppcno11.f}
 \index{ppcno12.f}\index{ppcno10BeBFe.f}\index{ppcno12BeBFe.f}
 \index{ppcno12Be.f}\index{ppcno12Li.f}\index{ppcno3a9.f}\index{ppcno3aco.f}
 \index{réactions thermonucléaires}\index{ppcno3a12Ne.f}\index{ppcno3acos.f}
\begin{itemize}
 \item {\tt NOM\_NUC}~: nom de la routine de calcul des taux des
 réactions thermonucléaires et d'initialiation de la composition
 chi\-mi\-que~:
 \begin{itemize}
  \item{\tt pp1}~: calcul simplifié du cycle PP (ne permet pas la diffusion
  mi\-cros\-co\-pi\-que).
  \item{\tt pp3}~: réactions PP, 3 éléments H, He3, He4 avec H2, Li7, Be7 à
  l'équilibre; a pour vocation les tests de mise au point,
  permet la diffusion mi\-cros\-co\-pi\-que, intervalle de tabulation :
  [1MK,80MK]..
  \item{\tt ppcno9}~: réactions PP+CNO, 9 éléments, H2, Li7, Be7 à
  l'équilibre, intervalle de tabulation :[1MK,40MK].
  \item{\tt ppcno9Fe}~:  réactions PP+CNO, 9 éléments+Fe, H2, Li7, Be7 à
  l'équilibre, intervalle de tabulation :[1MK,80MK].
  \item{\tt ppcno10}~: réactions PP+CNO, 10 éléments, H2, Be7 à
  l'équilibre, intervalle de tabulation :[0.5MK,80MK].
  \item{\tt ppcno10Fe}~: réactions PP+CNO, 10 éléments+Fe, H2, Be7 à
  l'équilibre, intervalle de tabulation :[0.5MK,80MK].
  \item{\tt ppcno10K}~: réactions PP+CNO, 10 éléments+K, H2, Be7 à
  l'équilibre, intervalle de tabulation :[0.5MK,80MK].
  \item{\tt ppcno10BeBFe}~: réactions PP+CNO, 10 éléments + Li6, Be9,
   B11, Fe avec H2 et Be7 à l'équilibre, intervalle de tabulation
:[0.5MK,80MK].
  \item{\tt ppcno11}~: réactions PP+CNO, 11 éléments, Be7 à
  l'équilibre, intervalle de tabulation :[0.5MK,80MK].
  \item{\tt ppcno12}~: réactions PP+CNO, 12 éléments, intervalle de tabulation :
  [1MK,80MK].
  \item{\tt ppcno12Be}~: réactions PP+CNO, 12 éléments + Be9, intervalle de
   tabulation :[0.5MK,80MK].
  \item{\tt ppcno12Li}~: réactions PP+CNO, 12 éléments + Li6, intervalle de
   tabulation :[0.5MK,80MK].
  \item{\tt ppcno12BeBFe}~: réactions PP+CNO, 12 éléments + Li6, Be9,
   B11, Fe, intervalle de tabulation :[0.5MK,80MK].
  \item{\tt ppcno3a9}~: réactions PP+CNO+3$\alpha$, 9 éléments,
  H2, Li7, Be7 à l'équilibre, intervalle de tabulation :[1MK,80MK]..
  \item{\tt ppcno3a12Ne}~: réactions PP+CNO+3$\alpha$+carbone, 12 éléments
   dont Ne{22}, H2, Li7, Be7 à l'équilibre, intervalle de tabulation
:[1MK,800MK].
  \item{\tt ppcno3aco}~: réactions PP+CNO+3$\alpha$+carbone+néon+oxygène,
  17 éléments H2, Li7, Be7 à l'équilibre, intervalle de tabulation :[1MK,3GK].
  \item{\tt ppcno3acos}~: réactions PP+CNO+3$\alpha$+carbone+néon+oxygène,
  19 éléments H2, Li7, Be7 à l'équilibre, intervalle de tabulation :[1MK,3GK].
 \end{itemize}
 L'étendue de l'intervalle de tabulation est fixée en fonction de la zone
 d'efficacité du réseau nucléaire retenu.
 \index{nom\_nuc\_cpl}\index{NACRE}\index{Cau-Fow}\index{Adelb}\index{NACRU}
 \index{NACRE+}
 \item {\tt nom\_nuc\_cpl}~: nom de la compilation des réactions
  thermonucléaires~:
 \begin{itemize}
  \item {\tt Adelb}~: Compilation \citet{al98}.
  \item {\tt Cau-Fow}~: Compilation de \citet{cf88}.
  \item {\tt NACRE}~: Compilation de NACRE \citep{a99}.
  \item {\tt NACRE+}~: Compilation de NACRE \citep{a99} avec correction du
  3$\alpha$ selon \citet{o09}.
  \item {\tt NACRU}~: Compilation de NACRE updatée.
 \end{itemize}
 \item {\tt MITLER=.TRUE.}~: effet d'écran selon \citet{mi77}.
\end{itemize}

\section{Arguments de la NAMELIST {\tt NL\_ATM}}
\index{NL\_ATM}\index{tau\_max}\index{lim\_ro}\label{sec:nlatm}
\begin{itemize}
 \item {\tt NOM\_ATM}~: nom de la routine de restitution de l'atmosphère~:
 \index{restitution de l'atmosphère}\index{nom\_atm}\index{lim\_atm.f}
 \index{lim\_gong1.f}\index{lim\_tau1.f}
 \begin{itemize}
  \item {\tt lim\_atm}~: restitution de l'atmosphère à partir d'une
  loi $T(\tau)$,
  \item {\tt lim\_gong1}~: restitution simplifiée de l'atmosphère,
  cas de {\sc gong1},
  \item {\tt lim\_tau1}~: restitution simplifiée de l'atmosphère,
  atmosphère mono-couche.
 \end{itemize}
 \item {\tt NOM\_TDETAU}~: nom de la loi $T(\tau)$ utilisée pour
 la restitution de l'atmosphère~:
 \index{nom\_tdetau}\index{edding.f}\index{hopf.f}\index{k5750.f}\index{MARCS.f}
 \index{k5777.f}\index{roger00.f}\index{roger02.f}\index{roger05.f}
 \index{roger10a.f}\index{hsra.f}
 \begin{itemize}
  \item {\tt edding}~: loi $T(\tau)$ totalement radiative d'Eddington.
  \item {\tt hopf}~: loi $T(\tau)$ totalement radiative de Hopf.
  \item {\tt hsra}~: loi $T(\tau)$ empirique, non totalement radiative de HSRA.
  \item {\tt K5750}~: loi $T(\tau)$ non totalement radiative, dérivée
  du modèle d'atmosphère solaire de Kurucz, $T_{\rm eff}=5750$\,K.
  \item {\tt K5777}~: loi $T(\tau)$ non totalement radiative, dérivée
  du modèle d'atmosphère solaire de Kurucz, $T_{\rm eff}=5777$\,K.
  \item {\tt MARCS\_Z+0.00}~:  lois $T(\tau)$ non totalement radiatives,
  dérivées des modèles d'atmosphère de MARCS. Elles couvrent l'intervalle
  3000K $<=$ Teff $<=$ 8000K. Des tables\footnote{Codifiées MARCS\_Zsx.yy,
  Exemple: MARCS\_Z-2.00.} pour les
  métallicités $\rm\left[Z/X\right]$:\\
$ -3.00,\ -2.00,\ -1.50,\ -0.75,\ -0.50,\ -0.25,\ +0.00,\ +0.25,\ +0.75,
 \ +1.00$\\
 et des gravités ($\log_{10}\rm g$) recouvrant
l'intervalle $[+3.0,+5.0]$ avec un pas de $+0.5$ sont
disponibles sur simple demande auprès de Bernard.Pichon@oca.eu.\\
\index{MARCS}
Une table supplémentaire: {\tt MARCS\_Z+0.00\_NEW}, calculée avec
$\rm\left[Z/X\right]=0.0$, pour les conditions conditions de faible
 gravité et de forte température effective (atmosphères d'AGB) est aussi
 disponible.\index{MARCS}
\item {\tt piau1, piau3}~: lois $T(\tau)$ non totalement radiatives, dérivées
de modèles d'atmosphère de Kurucz (ATLAS12), de composition solaire As05({\tt
piau1}) ou As09({\tt piau3}), popI.\index{piau1/3}\index{conv\_cgm}
Pour la cohérence, utiliser la routine de convection Canuto Goldman Mazitelli
({\tt conv\_cgm}, $\alpha=0.5$) ApJ 473, 550, 1996. Les tables nécessaires à
 l'exploitation sont disponibles sur simple demande auprès de
Laurent.Piau@cea.fr, elles sont placées dans le sous-directory
{\tt SUN\_STAR\_DATA/PIAU}.

\noi Domaines couverts Teff / logg:
 \begin{verbatim}
 Teff	 logg 0.5  1.0  1.5  2.0  2.5  3.0  3.5  4.0  4.5  5.0  5.5
 4000        +    +    +    +    +    +    +    +    +    +    +
  |          +    +    +    +    +    +    +    +    +    +    +
  |          +    +    +    +    +    +    +    +    +    +    +
 6400        +    +    +    +    +    +    +    +    +    +    +
\end{verbatim}

\item {\tt roger+00}~: loi $T(\tau)$ non totalement radiative, dérivée de
modèles
 d'atmosphère de Kurucz, de métallicité $\rm [Fe/H]=0.0$.Pour l'exploitation
\item {\tt roger+02}~: loi $T(\tau)$ non totalement radiative, dérivée
de modèles d'atmosphère de Kurucz, de métallicité $\rm [Fe/H]=+0.2$.
\item {\tt roger-05}~: loi $T(\tau)$ non totalement radiative, dérivée
de modèles d'atmosphère de Kurucz, de métallicité $\rm [Fe/H]=-0.5$.
\item {\tt roger-10a}~: loi $T(\tau)$ non totalement radiative, dérivée
de modèles d'atmosphère de Kurucz, de métallicité $\rm [Fe/H]=-1.0$
avec éléments $\alpha$ "enhanced".
\end{itemize}
\index{tau\_max}\index{lim\_ro}
\item {\tt TAU\_MAX=10.d0}~: profondeur optique du raccord avec l'enveloppe,
\item {\tt LIM\_RO=.TRUE.}~: condition limite externe en densité,
sinon en pression.
\end{itemize}

\section{Exploitation sous UNIX \& LINUX}
\index{cesam2k.out}\index{UNIX}\index{LINUX}
Sous UNIX/LINUX, il est conseillé\footnote{Voir aussi \paragr{sec:espt}.}
de créer et de laisser le
module exécutable {\tt cesam2k.out} dans le
sous-directory {\tt SOURCE}. Pour l'exploitation, utiliser un sous-directory
propre à chaque application dans lequel on disposera les fichiers de
données ({\tt mon\_modele.don}), les fichiers d'initialisation tels que
{\tt m020.zams}, et, si besoin, les fichiers de réglage et de personnalisation.

\noi Par exemple, pour calculer des modèles de $\xi$\,Hya ($3M_\odot$) de la ZAMS
jusqu'à $400$\,My, saisir{\tt create\_evol\_space KSI\_HYA}. Le sous-directory
{\tt KSI\_HYA} est crée. Il contient tous les éléments nécessaires à
l'exploitation de \cesam2k ainsi que des README d'explication.
Une autre possibilité est de:
\begin{itemize}
\item Créer un sous-directory {\tt KSI\_HYA},
\item Y copier les fichiers {\tt m020.zams}\footnote{Utiliser des fichiers
d'initialisation ASCII ou binaire de $3M_\odot$ si on en dispose.}
et {\tt mon\_modele.don} du sous-directory {\tt EXPLOIT}.
\end{itemize}
Puis:
\begin{itemize}
\item Renommer le fichier de données {\tt mon\_modele.don} $-->$ {\tt ksi\_hya.don}.
\item Adapter la physique et les paramètres de ce fichier au problème de $\xi$\,Hya.
\end{itemize}
Lancer ensuite deux calculs sucessifs
d'un modèle préliminaire de ZAMS avec {\tt (AGEMAX = 0.d0)} adapté aux paramètres
et à la physique adoptés~:
\begin{enumerate}
\item le premier avec une précision réduite {\tt (PRECISION='NP')} et
en utilisant le modèle d'initialisation en ASCII,
\item et le second avec la précision requise par exemple {\tt (PRECISION='AV')} en utilisant le modèle binaire {\tt ksi\_hya\_B.hom}
 obtenu à l'étape précédente.
\end{enumerate}
Conserver le fichier permanent
{\tt ksi\_hya\_B.hom} obtenu à l'issue de cette dernière étape dans
le fichier {\tt ksi\_hya\_ini\_B.hom}. C'est avec ce dernier que pourront être
initialisés les calcul des modèles d'exploitation.

\noi Dans le fichier {\tt ksi\_hya.don} indiquer l'âge à atteindre~:
(AGEMAX = 400.d0) et lancer l'exécution~: {\tt cesam2k.out} en
l'initialisant avec {\tt ksi\_hya\_ini\_B.hom}.

\noi Pour faire par la suite des modèles en changeant un ou
plusieurs paramètres (âge, overshoot, \etc), il suffit d'adapter le fichier
{\tt ksi\_hya.don} à chaque cas particulier.

\subsection{Fichiers ASCII et binaires d'initialisation}
\index{fichiers d'initialisation}\index{définition de la ZAMS}
Une évolution peut être initialisée par un modèle de pré-séquence principale
ou par un modèle de main-séquence {\bf homogène}. Ces modèles sont sous forme de
fichiers ASCII ou binaires; ces derniers étant des sous-produits d'évolutions
ultérieures. Des fichiers ASCII sont disposés dans le sous-directory {\tt EXPLOIT}
et sont implémentés dans le directory de travail si la procédure
{\tt create\_evol\_space} est utilisée, \cfparagr{sec:tfonc}. L'initialisation
PMS permet de faire évoluer la composition chimique dès que le réseau nucléaire
est activé ($T\sim 1MK$). A l'âge de la ZAMS, époque où la production d'énergie
d'origine nucléaire égalise celle d'origine gravitationnelle, la composition
chimique n'est plus homogène alors qu'elle l'est dans modèle d'initialisation
de ZAMS homogène qui, de ce fait, n'a pas de sens physique en soi. Ces conditions
non physiques étant relaxées après quelques pas temporels, typiquement une
vingtaine. Les fichiers binaires permanents d'extensions {\tt .pms} et {\tt .zams}
créés lors d'initialisations ou d'évolutions avec PMS peuvent servir à des
initialisations ou à des reprises
ultérieures; il en est de même des fichiers binaires permanents d'extension
{\tt .hom} créés à l'issue d'initialisations sur la ZAMS. De
façon générale ces fichiers binaires servent à optimiser les initialisations.
\index{extension .hom}  \index{extension .pms} \index{extension .zams}

\subsection{Exemple d'évolution avec PMS}\label{sec:tfonc1}
On se place dans le directory {\tt TEST} créé au \paragr{sec:tfonc}.\\
Saisir {\tt cesam2k.out}\index{PMS}
\begin{tt}\begin{verbatim}
--------------------
CESAM speaks a bit of english if you include in the
working directory a file named 'langue' with the
statements : &NL_LANGUE
langue='english'
/
cf. aide_mem2k, chapter Personnalisation
-------------------

Pour arrèter : taper 0 puis RETURN
Pour poursuivre une évolution : taper 1 puis RETURN
Pour initialiser un modèle de ZAMS : taper 2 puis RETURN
Pour initialiser un modèle de PMS : taper 3 puis RETURN
3  <==
Le modèle initial de PMS est-il donné en binaire ? o/n
n  <==
CESAM utilise un modèle initial de PMS en ASCII
entrer le nom du fichier ASCII du modèle PMS initial
Exemples: 2d-2.pms(Tc=0.2MK), 5d-4.pms(Tc=0.5MK)
 8d-5.pms(Tc=1.0MK), 4d-2.pms (pour M*>10Msol)
 5d-4.pms  <==
CESAM utilise le modle ASCII : 5d-4.pms
entrer l'identificateur du modèle :
Exemple : mon_modele
identificateur des fichiers du modèle : mon_modele  <==

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

          MODELE DE STRUCTURE INTERNE calculé par CESAM2k version : V3.7.24

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

         Début du calcul le : 17 Décembre 2012 à 17h59

         NAMELISTS du fichier: mon_modele.don

 &NL_CESAM NOM_CHEMIN=~/SUN_STAR_DATA/

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

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


-----------------------------------
calcul du modèle initial de PMS
entrer la constante de contraction c_iben, exemples :
pour Tc=100 000 K, entrer c_iben= 2.d-2
pour Tc=500 000 K, entrer c_iben= 5.d-4
pour Tc=1 000 000 K, entrer c_iben= 8.d-5
pour M*>10 Msol, entrer c_iben= 4.d-2
 5.d-4  <==
----- Paramètres pour le modèle quasi-statique ----

saut de la fonction de répartition par couche = 5.000E-02
modification du nombre de couches si la constante de
répartition varie de plus de  5%

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

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

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

Pré Séquence Principale premier modèle : C= 5.000E-04
Rext= 1.652E+01, Lext= 6.893E+01, Tc= 5.010E+05
ok? entrer o/n, pour arrèter entrer s
o  <==
calcul d'un nouveau modèle avec 1.1 C

Intégration du modèle quasi-statique de PMS(début)
 .............................................

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

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

Pré Séquence Principale second modèle : C= 5.500E-04
Rext= 1.714E+01, Lext= 7.322E+01, Tc= 4.838E+05, dt= 9.701E-04
 ok? entrer o/n, pour arrèter entrer s
 o <==
Pas de chute de planètoïdes
 &NL_DEVICE
H=6.00000000,DH=1.50000000,LD=8.50000000,DL=2.00000000,XLEFT=1.79999995,
YBOT=1.39999998,DEVICE=/xw
                                          /

-----Restitution de l'atmosphère (début)-----------
Lrac= 7.322E+01, Rrac= 1.714E+01, Teff= 4.082E+03, grav= 9.325E+01, tau*=
6.454E-01
atmosphère iter.  1 err. max. 1.5E-05, var: Pt, couche 59, corr 1.0E+00
atmosphère iter.  2 err. max. 1.0E-09, var: Pt, couche 66, corr 1.0E+00
atmosphère iter.  3 err. max. 7.8E-16, var: M , couche 31, corr 1.0E+00
------- Restitution de l'atmosphère (fin) ------

*********
âge= 0.000E+00, LogTeff= 3.611E+00, LogL/Lsol= 1.865E+00, LogR/Rsol= 1.235E+00
Log g= 1.968E+00, Pc= 1.384E+11, Tc= 4.838E+05, Roc= 2.138E-03, Xc= 7.000E-01
ePP/eNUC=  0%, eCNO/eNUC=  0%, e3a+C+O/eNUC=  0%, eGRAV/eNUC=  100%, Yc=
2.799E-01
Var. rel. de masse :  0.000E+00, M*= 1.000E+00Msol, modèle de pré-série
principale
Modèle hors la zone instable / outside the instability strip (P= 3.016E+00day)
*********

-- Intégration du modèle quasi-statique (début)----

Modèle n°   1, pas temporel estimé, dt= 9.701E-04

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

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

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

Y.Lebreton OBSPM, P.Morel OCA. Décembre 2008, mixtures+opacités AGS05
Y.Lebreton OBSPM, P.Morel OCA. Novembre 2009, etat_irwin
T.Guillot, P.Morel OCA. Mai 2012, implantation de sorties *.csv

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

 Type <RETURN> for next page:
\end{verbatim}\end{tt}


\subsection{Limitation des sorties}\index{limitation des sorties}
\label{sec:baratine}\index{baratine}\index{blabla}
Les informations qui apparaissant {\em on line} sur le moniteur permettent de
suivre en détail le déroulement des calculs.
Elles peuvent s'avérer non relevantes pour une exploitation "batch" et encombrer
inutilement. A cet effet, la variable logique
{\tt baratine} du module {\tt mod\_donnees}, lorsqu'elle est {\tt .FALSE.},
permet de dérouter sur les fichiers
{\tt mon\_modele\_static, mon\_modele\_atmos, mon\_modele\_evol} la plupart des
les informations concernant le déroulement des calculs en ce qui concerne,
respectivement, la résolution de l'équilibre quasi-statique, la restitution de
l'atmosphère, l'évolution de la composition chimique et de la vitesse
angulaire.
Pour ce faire, il est nécessaire ce coder {\tt baratine=.FALSE.} dans le module
{\tt mod\_donnees} et d'effectuer une recompilation partielle de ce module,
\cfparagr{sec:exploit}. Quand la nature cette limitation de sorties n'est
qu'occasionnelle, il n'est pas nécessaire d'intervenir dans le module {\tt
mod\_donnees} et de refaire une compilation; il suffit de disposer, dans
l'environnement du calcul, le fichier {\tt blabla} décrit au
\paragr{sec:blabla}.

\section{Contenu du package EXPLOIT }
L'exploitation de \cesam2k nécessite l'utilisation d'outils
astrophysique, comme par exemple le calcul de la valeur du rapport $[Z/X]$
d'une mixture, ou encore de transformer une donnée photométrique et une
parallaxe en luminosité. Des programmes fournissant ces données sont regroupés
dans le package {\tt EXPLOIT} avec des programmes d'exploitation des résultats
des programmes de tracés par exemple ainsi que des shell-script de leur gestion
comme la calibration solaire. Ces programmes font appel aux
bibliothèques du sous directory {\tt SOURCE}.
Les exécutables ont été créés globalement lors de l'exécution du Makefile
{\tt make\_expl} \cfparagr{sec:mexec}, le Makefile {\tt  Make\_exe} permettant
de les créer individuellement.

\subsection{Outils astrophysiques}
Ces outils sont les suivants, une
description en est donnée dans leur sources respectives~:
\begin{description}
\item[calc\_nu\_reac:] calcul de taux de réactions et de leur dérivée/T
\item[calc\_alpha\_enhanced:] facteur d'amplification pour les opacités
\item[calc\_bv\_tycho:] B, V, Vh à partir des observations de Tycho
\item[calc\_corr\_bol:] corrections bolométriques et indices de couleur
\item[calc\_dnu0:] calcul de $\delta \nu_0$ et de $\nu_{\rm max}$
\item[calc\_dzdy:] $\Delta Z / \Delta Y$ à partir de Y et [Fe/H] initial et
final
\item[calc\_dtl:] calcul de $\Delta\log_{10} T_{\rm eff}$ et $\Delta L$
\item[calc\_fesh:] calcul de [Fe/H]
\item[calc\_fesh\_zsx:] Fe/H $\Longleftrightarrow$ Z/X
\item[calc\_glmt:] calcul de $\log_{10}g(L, T_{\rm eff}, M/M_\odot)$ ou de
 $L(\log_{10}g, T_{\rm eff}, M/M_\odot)$
\item[calc\_L2Mv :] calcul de Mv à partir de $L/L_\odot$ et inverse
\item[calc\_l\_mbol:] $L/L_\odot$ $\Longleftrightarrow$ $M_{\rm bol}$
\item[calc\_l\_r\_teff:] $L/L_\odot, R/R_\odot \Longrightarrow T_{\rm eff}$,
ou $L/L_\odot, T_{\rm eff}\Longrightarrow R/R_\odot$,
ou $R/R_\odot, T_{\rm eff}\Longrightarrow L/L_\odot$
\item[calc\_l\_g\_teff\_m:] $(T_{\rm eff},M/M_\odot,\log_{10}g)\Longrightarrow
L/L_\odot$ ou $(T_{\rm eff},M/M_\odot,L/L_\odot)\Longrightarrow\log_{10}g$
\item[calc\_lum\_teff:] $(\pi,T_{\rm eff},{\rm Ind.Coul.})\Longrightarrow
L/L_\odot$
\item[calc\_mag\_abs:] calcul de la magnitude absolue
\item[calc\_masses:] pour une binaire, masses individuelles avec erreurs à
partir de la somme et de la fraction de masse
\item[calc\_mbol\_lum:] $M_{\rm bol}$ et $L/L_\odot$ à partir de V et du
module de distance
\item[calc\_mod\_dis:] calcul du module de distance
\item[calc\_moyenne\_ponderee:] calcul de la moyenne pondérée et du sigma
\item[calc\_new\_x:] [Fe/H] initial et X pour avoir fesh0 a la fin d'une
évolution
\item[calc\_rltgm:] Détermination de 2 quantités parmi
$T_{\rm eff},\log_{10}g, M/M_\odot,R/R_\odot,L/L_\odot$ à partir des 3 autres
\item[calc\_r\_theta:] diamètre + $\pi$ en mas $\Longrightarrow R/R_\odot$
\item[calc\_som\_mag:] pour une binaire, calcul de la magnitude de A+B
\item[calc\_teff:]  $T_{\rm eff}$ à partir de calibrations photométriques
\item[calc\_vsini:] $R/R_\odot, V \sin i, i \Longrightarrow$ période de
rotation en jours, radians/sec et Hz
\item[calc\_XisH:] calcul des [Xi/H]
\item[calc\_XYZ:] {X, Y $\Longrightarrow$ Z et Z/X} ou {X, Z $\Longrightarrow$ Y
et Z/X} ou {Y, Z $\Longrightarrow$ X et Z/X}
\item[calc\_z\_sur\_x:] calcul de [Z / X ] pour différentes mixtures
\end{description}

%genereT_bibli_list-dbg outils list_outils
%exeT_test-dbg program

\subsection{Programmes d'exploitation}
D'autres programmes du sous directory {\tt EXPLOIT} permettent une
exploitation des fichiers ASCII ou binaires créés par \cesam2k~:
\index{calib2k\_pms.f}\index{des2k\_abon.f}\index{des2k\_hr.f}
\index{des2k\_osc.f}
\index{des2k\_dhve.f}\index{des2k\_abontc.f}\index{des2k\_abonts.f}
\index{des2k\_ZC.f}\label{sec:pexploit}
\index{calib2k\_zams.f}
\begin{itemize}
 \item {\tt calib2k\_pms}~: Formation du fichier de données {\tt
mon\_modele.don}
 pour calibration d'un modèle solaire avec PMS.
 \item {\tt calib2k\_zams}~: Formation du fichier de données {\tt
mon\_modele.don}
 pour calibration d'un modèle solaire initialisé sur la ZAMS.
 \item {\tt des2k\_ZC}~: Dessin de l'évolution des zones convectives en
 fonction du temps.
 \item {\tt des2k\_abon}~: Dessin des abondances en fonction de la masse
 ou du rayon.
 \item {\tt des2k\_abontc}~: Dessin de l'évolution de X, Y, Z au centre,
 en fonction du temps.
 \item {\tt des2k\_abonts}~: Dessin de l'évolution de X, Y, Z en surface,
 en fonction du temps.
 \item {\tt des2k\_bin, des2k\_bin\_tds}~: Dessin des variables quasi-statiques
 à partir d'un fichier binaire.\index{des2k\_bin.f}
 %\item {\tt des2k\_dhve}~: Dessin des coefficients des équations de la diffusion
% du moment cinétique.
 \item {\tt des2k\_diff\_spl}~: Dessin des différences entre deux modèles avec
 interpolation par B-splines.\index{des2k\_diff\_spl.f}
 \item {\tt des2k\_grad}~: Dessins des différents gradients d'un modèle.
 \index{des2k\_grad.f}
 \item {\tt des2k\_hr / des2k\_lr}~: Dessin du diagramme HR à partir d'un fichier d'extension {\tt .HR} / {\tt .LR} .
 \item {\tt des2k\_opa}~: Dessin des différences relatives entre deux tables
 d'\opa.\index{des2k\_opa.f}
 \item {\tt des2k\_osc}~: Dessin de 1, 2 ou 3 modèles en fonction du rayon,
 de la masse avec des possibilités de zoom.
 \item {\tt des2k\_rot}~: Dessins séparés des variables de la diffusion du
 moment cinétique.
 \item {\tt des2k\_rot\_ext}~: Dessin du profil, en fonction du temps, de la
 vitesse angulaire de la couche externe en km/s ou en rad/sec.
 \item {\tt des2k\_vaiss}~: Dessin du profil de la \fbv.
 \index{fichier\_vent.f}
 \item {\tt fichier\_vent}~: Programme de construction d'un fichier de
 composition chi\-mi\-que du vent.
\end{itemize}
Dans chacun de ces programmes, des informations concernant leur utilisation sont
données sous la forme de commentaires.

\subsection{Shell-scripts d'exploitation}\label{sec:exploit}
Des shell-scripts {\tt bash} permettent
l'exploitation\footnote{Les options de compilation sont
celles du compilateur LINUX LF95, il convient de les adapter à celles du
compilateur utilisé.} de \cesam2k essentiellement la calibration solaire  et de
l'évolution à partir de la ZAMS ou de la PMS.
\index{evol2k\_zams}
\index{calib2k\_pms}\index{evol2k\_pms}
\begin{itemize}
 \item {\tt calib\_pms}~: procédure de calibration d'un modèle solaire
 initialisé sur la PMS homogène.
 \item {\tt calib\_zams}~: procédure de calibration d'un modèle solaire
 initialisé sur la ZAMS homogène.
 \item {\tt create\_evol\_space}~: création d'un espace de travail contenant les
 fichier nécessaires à l'exploitation de \cesam2k, les fichiers de données et de
 réglage types et des modèles initiaux de PMS et de ZAMS.
 \item {\tt evol2k\_pms}: script pour une  évolution initialisée  sur la PMS
 \item {\tt evol2k\_zams}: script pour une évolution initialisée  sur la ZAMS
 \item {\tt evolsun\_pms} : script pour une évolution solaire initialisée sur la PMS
 \item {\tt evolsun\_pms\_list} : évolution d'une liste de modèles solaires
 initialisés sur la PMS
 \item {\tt evolsun\_zams} : script pour une évolution solaire initialisée sur la ZAMS
 \item {\tt evolsun\_zams\_list} : évolution d'une liste de modèles solaire initialisés sur la ZAMS
 \item {\tt Make\_exe} : makefile pour compilation et création d'un exécutable
 \item {\tt Make\_expl} : makefile pour compilation et création des exécutables d'EXPLOIT
 \item {\tt Make} : makefile pour compilation et création des exécutables de la SOURCE
\end{itemize}

\subsection{Données diverses}\label{sec:dondiv}
\begin{itemize}
 \item {\tt 2d-2\_pms}~: fichier de modèle initial de PMS, avec constante de
 contraction {\tt c\_iben}=$2.10^{-2}$, soit $T_c \sim 1.10^5$ K.
 \item {\tt 4d-2\_pms}~: fichier de modèle initial de PMS, avec constante de
 contraction {\tt c\_iben}=$3.10^{-2}$, soit $T_c \sim 3.10^5$ K.
 \item {\tt 5d-4\_pms}~: fichier de modèle initial de PMS, avec constante de
 contraction {\tt c\_iben}= $5.10^{-4}$, soit $T_c \sim 5.10^5$ K (combustion du
 deutérium).
 \item {\tt 8d-5\_pms}~: $8.10^{-5}$ fichier de modèle initial de PMS, avec constante
 de contraction {\tt c\_iben}=, soit $T_c \sim 1.10^6$ K (amorçage des réactions
 nucléaires du cycle PP).
 \item {\tt blabla}~: exemple de NAMELIST "NL\_BLABLA" pour la limitation des
 écritures lors d'une exécution.
 \item {\tt device}~: exemple de NAMELIST "NL\_device" permettant le cadrage des graphes "on line".
 \item {\tt langue}~: exemple de NAMELIST "NL\_LANGUE" définissant la langue à
 utiliser pour les commentaires produits en cours d'exécution, \cfparagr{sec:lang}.
 \item {\tt m010.zams}~: fichier de modèle initial de ZAMS, de $1.M_\odot$.
 \item {\tt m020.zams}~: fichier de modèle initial de ZAMS, de $2.M_\odot$.
 \item {\tt m050.zams}~: fichier de modèle initial de ZAMS, de $5.M_\odot$.
 \item {\tt m150.zams}~: fichier de modèle initial de ZAMS, de $15.M_\odot$.
 \item {\tt mixture}~: exemple de NAMELIST "NL\_MIXTURE" pour définition d'une
 nouvelle mixture initiale, \cfparagr{sec:mix}.
 \item {\tt mon\_modele\_modif\_mix}~: exemple de NAMELIST "nl\_modif\_mix" de
 modification d'une mixture initiale, \cfparagr{sec:rapab}.
 \item {\tt mon\_modele.don}~: exemple de fichier de données.
 \item {\tt mon\_modele.rap\_iso}~: exemple de NAMELIST "NL\_rap\_iso" de modification
 des rapports isotopiques initiaux, \cfparagr{sec:rapiso}.
 \item {\tt mon\_modele.rg}~:  valeurs par défaut du fichier de réglages, \cfparagr{sec:reg}.
 \item {\tt mon\_modele.vent}~: exemple de de NAMELIST "NL\_VENT" de composition
 chimique du vent, \cfparagr{sec:vent}.
 \item {\tt mon\_modele.zoom}~: exemple de de NAMELIST "NL\_DES" pour adaptation
 du graphe on line.
  \item {\tt opteos}~: exemple de la NAME\_LIST "nl\_oes" permettant d'adapter les arguments de l'équation d'état Irwin, \cfparagr{
sec:pfree}.
 \item {\tt planet}~: exemple de de NAMELIST "NL\_PLANET" de notification des
  caractéristiques de planétoïdes, \cfparagr{sec:cplanet}.
 \item {\tt vson\_rho\_sun\_data}~: profils de la vitesse du son et de la densité
 solaires.

\end{itemize}
