\label{sec:perso}
Bien que l'exécutable soit construit un fois pour toutes, il est possible
d'adapter de façon externe un certain nombre de paramètres. Ces
adaptations s'effectuent, à la demande, par l'intermédiaire de
fichiers {\em ad hoc} à disposer dans le sous-directory où s'effectue le calcul.
Ne seront concernés que les modèles calculés dans cet environnement.
De façon générale, dans le directory concerné, une personnalisation peut
affecter, soit un modèle particulier, soit tous les modèles à calculer
dans ce directory. \cesam2k recherche d'abord si une
personnalisation existe pour le modèle à calculer puis, à défaut, si elle
existe pour tous les modèles du directory. Des exemples de
tels fichiers de personnalisation se trouvent dans
le sous-directory {\tt EXPLOIT} de la distribution. On peut ainsi~:
\begin{itemize}
 \item Avoir des commentaires d'interactivité en une langue différente du
 français\footnote{Seulement en anglais dans la présente version.}.
 \item Utiliser une mixture initiale différente de celles implémentées.
 \item Modifier les rapports d'abondances d'une mixture.
 \item Redéfinir les rapports isotopiques.
 \item Créer un fichier de sortie en ASCII contenant des quantités
 en nombre et ordre différents de ce qui est implémenté.
 \item Utiliser des réglages personnalisés.
 \item Dessiner des boites correspondant à des objectifs à atteindre,
 et adapter les échelles du diagramme HR tracé ``on line''.
 \item Adapter les dimensions des cadres et le nom du device à utiliser.
 \item Détourner sur des fichiers en ASCII une partie des informations
  apparaissant 'on line' sur le moniteur.
\end{itemize}
{\embf Par sécurité}, ainsi qu'il est explicité pour chaque cas particulier,
 certains des fichiers de personnalisation concernant la composition
chimique ne seront pris en compte que si le paramètre
{\tt modif\_chim} de la NAMELIST {\tt nl\_chim} du fichier de données est
{\tt .TRUE.} \cf \paragr{sec:nlchim}.

\section{Langue}\label{sec:lang}\index{langue}\index{anglais}
Une alternative pour avoir des commentaires d'interactivité
en une langue différente du français~:
\begin{enumerate}
 \item Si on désire {\bf toujours} avoir des commentaires
 en une langue différente du français,
 {\bf avant} la compilation, dans la routine {\tt cesam} du
 sous-directory {\tt SOURCE}, commenter ou supprimer les instructions
 (vers les lignes 270-278)~:
 \begin{tt}\begin{verbatim}
      INQUIRE(file='langue',exist=ok)
      IF(ok)THEN
       OPEN(unit=30,form='formatted',status='old',file='langue')
       READ(30,nl_langue) ; CLOSE(UNIT=30)
      ELSE
       WRITE(*,30)
30     FORMAT(/,'Sometimes CESAM can speak english,',/,
     1 'see aide_mem2k, chapter personnalisation ')
       langue='francais'
      ENDIF
\end{verbatim}\end{tt}
et décommenter l'instruction de la ligne suivante (279)~:
 \begin{tt}\begin{verbatim}
c       langue='english'
\end{verbatim}\end{tt}

Les commentaires en anglais seront désormais utilisés pour toute exploitation.
\item Si on désire n'avoir qu'occasionnellement les commentaires
 en anglais, introduire dans l'environnement
 d'exploitation le fichier de nom {\tt langue} que l'on trouvera dans
 le sous directory EXPLOIT, ou encore l'y créer, \cfparagr{sec:dondiv}~:
 \begin{tt}\begin{verbatim}
&NL_LANGUE
langue='english'
/
\end{verbatim}\end{tt}
Dans toutes les exploitations effectuées dans cet environnement
les commentaires apparaitront en anglais.
\end{enumerate}
Par défaut, \ie\ en l'absence de l'une des dispositions précédentes,
les commentaires seront en français.

\medskip
\PB{Dans la version actuelle de \cesam2k, seulement
une partie des principaux commentaires a été traduite en anglais.}

\medskip
\PB{{\bf Call for man power }~: il est facile de transposer dans
une autre langue les algorithmes
créés pour l'anglais ou encore compléter et/ou améliorer la syntaxe des
ces derniers, les volontaires sont bienvenus.}

\section{Mixture}\label{sec:mix}\index{mixture}\index{extension .mix}
\index{option mixture}
Pour utiliser une mixture initiale différente\footnote{Voir aussi
\cfparagr{sec:mix2}.} de celles implémentées,
coder {\tt nom\_abon='mixture'} dans le fichier de données et,
 dans l'environnement, placer un fichier de nom
{\tt mixture} contenant la mixture désirée en Dex. \cesam2k
{\embf refusera} toute mixture dont la normalisation est différente de H=12.
L'exemple suivant se trouve dans le sous-directory {\tt EXPLOIT}, \cfparagr{sec:dondiv}~:
\begin{footnotesize}\begin{tt}\begin{verbatim}
&NL_MIXTURE
ab(1)=12.00d0 , ab(2)=10.d0 , ab(3)=2.69d0  , ab(4)=2.15d0  , ab(5)=2.60d0
ab(6)=4.55d0  , ab(7)=3.97d0  , ab(8)=4.87d0  , ab(9)=1.56d0  , ab(10)=4.08d0
ab(11)=2.33d0 , ab(12)=3.58d0 , ab(13)=2.47d0 , ab(14)=3.55d0 , ab(15)=1.45d0
ab(16)=3.21d0 , ab(17)=1.5d0  , ab(18)=2.52d0 , ab(19)=1.12d0 , ab(20)=2.36d0
ab(21)=-1.17d0 , ab(22)=1.02d0 , ab(23)=0.d0   , ab(24)=1.67d0 , ab(25)=1.39d0
ab(26)=3.5d0  , ab(27)=0.92d0 , ab(28)=2.25d0
/
\end{verbatim}\end{tt}\end{footnotesize}
On trouve, dans la routine {\tt abon\_ini}, l'ordre et l'identification
des éléments par leurs indices, en fait, celle des opacités OPAL.
Dès lors, {\embf tous} les modèles calculés dans ce directory, s'ils utilisent
l'option {\tt nom\_abon='mixture'}, auront des abondances initiales conformes à
celles du fichier {\tt mixture}. Une autre possibilité permet de
personnaliser davantage; elle consiste à donner au fichier {\tt mixture}
le nom générique du modèle avec l'extension {\tt .mix},
Ex~: {\tt mon\_modele.mix}. Dès lors, {\embf seul} le modèle {\tt mon\_modele}
aura la composition chimique initiale désirée.
Avec l'option {\tt mixture}, \cesam2k cherche {\embf en priorité}
 à utiliser le fichier
{\tt mon\_modele.mix}; en son absence, c'est le fichier {\tt mixture} qui
est utilisé; le calcul n'est effectué que si l'un de ces deux fichiers
existe dans l'environnement.\\
Cette option nécessite de coder {\tt MODIF\_CHIM=.TRUE.}, \cfparagr{sec:nlchim}.

\medskip
\PB{Avec l'option {\tt nom\_abon='mixture'}, les abondances relatives des
métaux devraient, en principe, correspondre à celles des opacités,
au moins, faute de mieux, pouvoir être
prises en charge par la routine d'opacité utilisée.}

\section{Rapports isotopiques}\label{sec:rapiso}
\index{rapports isotopiques}
Pour utiliser des rapports isotopiques différents de ceux pris
par défaut par \cesam2k, placer dans l'environnement du calcul, un fichier de nom
{\tt mon\_modele.rap\_iso} dans lequel seront précisés les rapports
isotopiques à utiliser. Si {\embf tous} les modèles de l'environnement
doivent être calculés avec des rapports mo\-di\-fiés, donner au fichier
précédent le nom de {\tt rap\_iso}, \cfparagr{sec:dondiv}. Exemple~:
\begin{footnotesize}\begin{tt}\begin{verbatim}
&NL_RAP_ISO
be7sbe9=1.d-25,
be7sz=1.d-29,
c13sc12=1.10d-2,
h2sh1=3.01d-5,
he3she4=1.1d-4,
he3she4z=4.185d-4,
li6sli7=7.5d-2,
mg25smg24=0.0125,
mg26smg25=0.013,
ne22sne20=6.79d-2,
n15sn14=0.366d-2,
o17so16=0.038d-2
o18so16=0.008d-2
/
\end{verbatim}\end{tt}\end{footnotesize}
Les noms utilisés pour les variables suggèrent les notations des rapports
isotopiques {\tt be7sbe9}~: rapport isotopique \Be7/\Be9.
{\tt he3she4z} est le rapport isotopique \He3/\He4 sur la ZAMS lorsque le
deutérium initial a été transformé en \He3, {\tt be7sz} est l'abondance,
quasi nulle, du \Be7 initial dans Z.\\
{\embf Par sécurité}, ce fichier ne sera pris en compte, que si le paramètre
{\tt modif\_chim} de la NAMELIST {\tt nl\_chim} du fichier de données est
{\tt .TRUE.} \cfparagr{sec:nlchim}.

\section{Personnalisation des rapports d'abondances}\label{sec:rapab}
\index{abondances}\index{rapport des abondances}\index{abondance d'hélium}
\index{abondances initiales}\index{NL\_MODIF\_MIX}
\index{add\_Li}\index{add\_Be}\index{add\_B}\index{add\_F}\index{add\_Ne}
\index{add\_C}\index{add\_N}\index{add\_O}\index{add\_Na}\index{add\_Mg}
\index{add\_Al}\index{add\_Si}\index{add\_P}\index{add\_S}\index{add\_Cl}
\index{add\_Ar}\index{add\_K}\index{add\_Ca}\index{add\_Sc}\index{add\_Ti}
\index{add\_V}\index{add\_Cr}\index{add\_Mn}\index{add\_Fe}\index{add\_Co}
\index{add\_Ni}\index{add\_Z}\index{modif\_xish}
La mixture une fois définie, par exemple {\tt solaire\_gn}
\citep{gn93}, il est possible d'y modifier les abondances.
Pour ce faire, créer dans l'environnement du calcul, un fichier de nom
{\tt mon\_modele.modif\_mix}, dans lequel seront indiquées les modifications
en DeX à apporter. Si tous les modèles de l'environnement doivent être
calculés avec ces modifications, donner au fichier précédent le nom de
{\tt modif\_mix}. Exemple se trouvant dans le directory EXPLOIT, \cfparagr{sec:dondiv}~:
\begin{footnotesize}\begin{tt}\begin{verbatim}
&nl_modif_mix
add_Li=2.d0 ; add_Be=0.d0 ; add_B=0.d0
add_C=-1.0d0 ; add_N=0.d0 ; add_O=-0.3d0 ; add_F=0.d0 ; add_Ne=0.d0
add_Na=0.d0 ; add_Mg=0.d0 ; add_Al=0.d0 ; add_Si=0.d0 ; add_P=0.d0
add_S=0.d0 ; add_Cl=0.d0 ; add_Ar=0.d0 ; add_K=0.d0 ; add_Ca=0.d0
add_Sc=0.d0 ; add_Ti=0.d0 ; add_V=0.d0 ; add_Cr=0.d0 ; add_Mn=0.d0
add_Fe=0.d0 ; add_Co=0.d0 ; add_Ni=0.d0 ; add_Z=-1.d0
/
\end{verbatim}\end{tt}\end{footnotesize}
Les notations suggèrent les quantités à modifier; ainsi {\tt add\_Li=2.d0}
est le nombre de DeX à ajouter à l'abondance du lithium.
{\tt add\_Z=-1.d0} signifie que les
abondances {\embf de tous les métaux} seront diminuées d'un DeX.
{\embf Par sécurité}, ce fichier ne sera pris en compte que si le paramètre
{\tt modif\_chim} de la NAMELIST {\tt nl\_chim} du fichier de données est
{\tt .TRUE.} \cf \paragr{sec:fdon}.

Si on a codé {\tt GARDE\_XISH=.FALSE.} dans le fichier de
données \cfparagr{sec:msh},
le paramètre {\tt add\_Z}, qui agit sur les métaux dans leur
ensemble, n'a pas d'incidence dans la détermination des
rapports métal/Z. A l'opposé, une valeur non nulle de {\tt add\_C}, par
exemple, aura une incidence sur les rapports métal/Z.

\section{Personnalisation de la composition chimique du vent.}
\label{sec:vent}\index{vent.f}\index{perte de masse}\index{gain de masse}
Dans le cas standard, avec la donnée {\tt mdot} non nulle, \cfparagr{sec:fdon},
la composition chimique de la masse perdue ou gagnée par le vent est celle
de la couche la plus externe du modèle.
Il est possible de personnaliser cette composition chi\-mi\-que,
afin qu'elle diffère de celle de la couche externe. Pour ce faire, placer
dans l'environnement du calcul, un fichier de nom
{\tt mon\_modele.vent} dans lequel seront indiquées les fractions de masse
de chaque élément chimique dans l'apport ou la perte de masse.
Si tous les modèles de l'environnement doivent être
calculés avec ces modifications, donner au fichier précédent le nom de
{\tt vent}. Un exemple se trouve dans le sous-directory EXPLOIT, \cfparagr{sec:dondiv}~:
\begin{footnotesize}\begin{tt}\begin{verbatim}
&nl_vent
vt_H=0.7347,vt_He=0.2483,vt_Li=1.033e-08,vt_Be=1.72816e-10,
vt_B=4.859e-09,vt_C=2.899e-03,vt_N=8.493e-04,vt_O=7.885e-03,
vt_F=4.182e-07,vt_Ne=1.768e-03,vt_Na=3.501e-05,vt_Mg=6.736e-04,
vt_Al=6.078e-05,vt_Si=7.434e-04,vt_P=8.198e-06,vt_S=3.704e-04,
vt_Cl=4.924e-06,vt_Ar=7.315e-05,vt_K=3.845e-06,vt_Ca=6.541e-05,
vt_Sc=4.126e-08,vt_Ti=3.041e-06,vt_V=3.888e-07,vt_Cr=1.856e-05,
vt_Mn=1.357e-05,vt_Fe=1.287e-03,vt_Co=3.492e-06,vt_Ni=7.6092e-05
/
\end{verbatim}\end{tt}\end{footnotesize}
Les notations sont suggestives, ainsi {\tt vt\_Li=1.033e-08}
est la fraction masse du lithium dans le vent. Une renormalisation étant
effectuée lors de l'utilisation,
il n'est pas nécessaire de normaliser les fractions de masse figurant dans
ce fichier, certaines peuvent être nulles et même négatives.
Suivant le signe de {\tt MDOT}, la quantité de matière, apportée ou
soustraite, est supposée provenir de la \ZC\ externe qui existe toujours.
Le fichier {\tt vent} peut être construit\footnote{La composition chimique de
l'exemple est celle de la mixture météoritique de \citet{gs98}. Par ailleurs,
il n'est pas nécessaire, ici, de respecter la 'case'.} par le programme
{\tt fichier\_vent} du sous-directory EXPLOIT.

{\embf Par sécurité}, ce fichier ne sera pris en compte, que si le paramètre
{\tt modif\_chim} de la NAMELIST {\tt nl\_fdon} du fichier de données est
{\tt .TRUE.} \cf \paragr{sec:nlchim}.

\section{Chute de planètoïdes}\label{sec:cplanet}
\index{planètoïdes}
Au cours d'une évolution \cesam2k permet
de simuler une chute de planètoïdes sur un intervalle de temps limité.
Pour ce faire, placer
dans l'environnement du calcul, un fichier de nom
{\tt mon\_modele.planet} dans lequel seront indiquées les fractions de masse
de chaque élément de la composition chimique des planètoïdes, le nombre total de
planètoïdes qui sera reçu par l'étoile, l'intervalle de temps concerné
par la chute et le profil de l'intensité de la chute.
Si tous les modèles de l'environnement doivent être
calculés avec ces modifications, donner au fichier précédent le nom de
{\tt planet}. Un exemple se trouve dans le sous-directory EXPLOIT, \cfparagr{sec:dondiv}~:
\begin{footnotesize}\begin{tt}\begin{verbatim}
&nl_planet
vt_H=0.7347,vt_He=0.2483,vt_Li=1.033e-08,vt_Be=1.72816e-10,
vt_B=4.859e-09,vt_C=2.899e-03,vt_N=8.493e-04,vt_O=7.885e-03,
vt_F=4.182e-07,vt_Ne=1.768e-03,vt_Na=3.501e-05,vt_Mg=6.736e-04,
vt_Al=6.078e-05,vt_Si=7.434e-04,vt_P=8.198e-06,vt_S=3.704e-04,
vt_Cl=4.924e-06,vt_Ar=7.315e-05,vt_K=3.845e-06,vt_Ca=6.541e-05,
vt_Sc=4.126e-08,vt_Ti=3.041e-06,vt_V=3.888e-07,vt_Cr=1.856e-05,
vt_Mn=1.357e-05,vt_Fe=1.287e-03,vt_Co=3.492e-06,vt_Ni=7.6092e-05,
ypl=0.d0,zpl=0.9d0,
n_planet=10,profil='rectangle',age_deb=5.d0,age_fin=10.d0,
r_giration=1.d0,t_giration=1.d0
/
\end{verbatim}\end{tt}\end{footnotesize}
Les notations sont suggestives.
Avec {\tt vt\_H} et {\tt vt\_He} {\embf non nuls}, des valeurs non nulles de
{\tt ypl} et/ou {\tt zpl} permettent de modifier les fractions de masse de
l'hydrogène {\tt X}, de l'hélium {\tt Y} et des métaux {\tt Z}
des planètoïdes. Avec le fichier {\tt planet} précédent,
l'abondance en masse de l'hélium sera $Y=0$, celle des métaux $Z=0.9$ et
celle d'hydrogène $X=1.0-0.0-0.9=0.1$. {\tt age\_deb} et {\tt age\_fin} sont
respectivement les âges, en million d'années, du dédut et de la fin de
la chute. {\tt r\_giration} est le rayon de giration en UA et {\tt t\_giration}
le temps de giration, en années. \cesam2k estime
la vitesse angulaire des planètoïdes en admettant qu'ils décrivent un arc de
$\pi/2$ pendant un temps de giration. Une valeur négative du temps de giration,
correspond à un apport négatif de moment cinétique, i.e. rétrograde.
Avec une valeur nulle du temps de giration il n'y a pas d'apport de moment
cinétique.\\
Dans la routine {\tt planetoides.f}, \cesam2k propose 4 profils~:
\begin{enumerate}
 \item {\tt rectangle}~: L'intensité de la chute est constant.
 \item {\tt triangle}~: Le profil de l'intensité de la chute est triangulaire,
 le maximum étant centré au milieu de l'intervalle de temps.
 \item {\tt parabole}~: Le profil de l'intensité de la chute est parabolique,
 le maximum étant centré au milieu de l'intervalle de temps.
 \item {\tt gauss}~: Le profil de l'intensité de la chute est une gaussienne
 limitée à l'intervalle de temps.
 Le maximum est centré au milieu de l'intervalle de temps. L'écart type est
 arbitrairement pris égal au tiers de l'intervalle de temps.
\end{enumerate}
Cette option nécessite de coder {\tt MODIF\_CHIM=.TRUE.}, \cfparagr{sec:nlchim}.

\section{Personnalisation des paramètres de free-eos}
\label{sec:pfree}\index{etat\_irwin.f}
Il est possible de personnaliser les paramètres de l'équation d'état {\tt free-eos}
en plaçant dans l'environnement un fichier de nom {\tt opteos} avec les valeurs
 désirées pour les paramètres, \cfparagr{sec:dondiv}~:
\begin{verbatim}
&nl_eos
ifoption=2,
ifmodified=1,
ifion=-1
/
\end{verbatim}
L'implantation de cette équation d'état est décrite au \paragr{sec:irwin} et
son exploitation dans \cesam2k au \paragr{sec:nletat}. L'utilisation des
paramètres est décrite dans les READMEs de la source
{\tt http://freeeos.sourceforge.net}, {\tt free\_eos-2.2.1.tar.gz}.
Quelques indications sont données dans la routine {\tt etat\_irwin.f} du
sous-directory {\tt SOURCE} de \cesam2k.

\section{Personnalisation du fichier Ascii}\label{sec:ascii}\index{ascii.f}
Pour créer un fichier de sortie en ASCII personnalisé, coder
{\tt nom\_output='ascii'} dans le fichier de données et,
 dans l'environnement, placer un fichier de nom
 {\tt sortie\_ascii}\footnote{Dont un exemple se trouve dans le
 sous-directory {\tt EXPLOIT}.} dans
 lequel sont définis les paramètres créant la structure souhaitée.
 Le codage, {\embf délicat}, est
 décrit dans la routine {\tt output.f}, les indices des variables sont
 identifiés dans la routine {\tt ascii.f} du sous-directory {\tt SOURCE}.
 Dans l'exemple~:
\begin{tt}\begin{verbatim}
12  3
1 2 3 4 5 6 9 10 11 12 13 15
1 2 3
.FALSE.
.TRUE.
.TRUE.
Fichier pour sortie ascii:
\end{verbatim}\end{tt}
les quantités ont les significations suivantes~:
\begin{itemize}
 \item 12 3 $\longrightarrow$ on écrira 12 variables globales et 3
 variables locales,
\item  1 2 3 4 5 6 9 10 11 12 13 15 $\longrightarrow$ indices
 des 12 variables globales (mstar*msol, rtot*rsol, ltot*lsol, z0, x0,
 alpha, X dans ZC, Y dans ZC, d2p, d2ro, age, vsal, w\_rot initial).
\item  1 2 3 $\longrightarrow$ indices des 3 variables locales~: Rayon,
 $\log$(M/Mtot), Température).
\item .FALSE. $\longrightarrow$ la masse sera écrite en fraction de la
 masse totale et non en DeX.
\item .TRUE. $\longrightarrow$ tabulation du centre à la surface.
\item .TRUE. $\longrightarrow$ on ajoute la composition chimique locale à
 la suite des variables locales.
 \item Fichier pour sortie ascii~: texte d'entête du fichier de sortie.
\end{itemize}
Le nom générique du fichier de données est repris pour former celui
du fichier de
sortie en ASCII~: avec {\tt mon\_modele.don} le nom du fichier de sortie
est {\tt mon\_modele-ascii}. Avec le fichier cité comme exemple, on a obtenu~:
\begin{scriptsize}\begin{tt}\begin{verbatim}
Fichier pour sortie ascii: mon_modele-ascii
CESAM2k version 0.0.0.0 lagr colloc  1  2 np no diffus, 06 Juillet 2003 17h46
Physique utilisée: etat_eff, opa_yveline_lisse, conv_jmj, ppcno9, NACRE
solaire_gn, lim_atm, hopf, perte_ext, diffm_mp, difft_nu, ctes_94
 10  H1  He3  He4  C12  C13  N14  N15  O16  O17  Si28
       541        12        13        10        -1
 1.989190000000E+33 6.010741880620E+10 2.478601588572E+33 2.000000000000E-02 7.000000000000E-01
 1.800000000000E+00 7.000000000000E-01 2.800000000000E-01-4.241618755712E+01-2.519592625187E+01
 0.000000000000E+00 0.000000000000E+00
 0.000000000000E+00 0.000000000000E+00 1.359013041915E+07 7.000000000000E-01 8.826927693868E-05
 2.799117307231E-01 3.425410508692E-03 4.128413464627E-05 1.059170894932E-03 4.168107978225E-06
 9.641655715570E-03 3.903146720978E-06 5.824407491460E-03
 1.302042603905E+09 1.175206179724E-08 1.353639070045E+07 7.000000000000E-01 8.826927693868E-05
 2.799117307231E-01 3.425410508692E-03 4.128413464627E-05 1.059170894932E-03 4.168107978225E-06
 9.641655715570E-03 3.903146720978E-06 5.824407491460E-03
           ......................................................

 6.014328252162E+10 1.000000000164E+00 4.543950048106E+03 7.000000000000E-01 8.826927693868E-05
 2.799117307231E-01 3.425410508692E-03 4.128413464627E-05 1.059170894932E-03 4.168107978225E-06
 9.641655715570E-03 3.903146720978E-06 5.824407491460E-03
 6.014636429604E+10 0.000000000000E+00 4.543482782113E+03 7.000000000000E-01 8.826927693868E-05
 2.799117307231E-01 3.425410508692E-03 4.128413464627E-05 1.059170894932E-03 4.168107978225E-06
 9.641655715570E-03 3.903146720978E-06 5.824407491460E-03
\end{verbatim}\end{tt}\end{scriptsize}
Comme pour les abondances initiales \cf \paragr{sec:mix}, il est possible
de ne créer ce fichier de sortie personnalisé, que pour un modèle particulier,
en donnant au fichier de sortie {\tt sortie\_ascii} le nom du modèle avec
l'extension {\tt .ascii}, exemple~: {\tt mon\_modele.ascii}. S'il existe dans
l'environnement, ce seront les paramètres de ce dernier fichier qui seront pris en
compte pour {\embf le} modèle de nom {\tt mon\_modele}. En son absence,
\cesam2k cherchera à utiliser ceux du fichier {\tt sortie\_ascii}, à défaut
de ce dernier, il n'y aura pas de sortie ASCII.

\section{Réglages}\label{sec:reg}\index{réglages}

\PB{Les arguments du fichier de réglage sont indiqués das le fificier
{\tt arguments\_des\_reglages.f} des sous directories {\tt SOURCE} et
{\tt EXPLOIT}.}
\medskip

Les réglages, prédéfinis dans la routine {\tt cesam.f} en fonction du paramètre
{\tt PRECISION} du fichier de données {\tt mon\_modele.don}, peuvent être
personnalisés en codant~: {\tt PRECISION='rg'}, et en disposant dans
l'environnement un fichier de nom~:
{\tt reglages}\footnote{Dont un exemple se trouve dans le sous-directory
{\tt EXPLOIT}.} comportant les réglages à utiliser.
{\bf Ces valeurs prévalent sur celles prédéfinies.}
Exemple d'un tel fichier {\tt reglages}~:
\index{NL\_RLG}\index{rg}
\begin{tt}\begin{verbatim}
&NL_RLG
m_qs=1,
m_ch=2,
m_rot=2,
m_tds=2,
m_ptm=2,
ordre=4,
precic=1.d-4,
precix=1.d-3,
precit=0.15d0,
psi0=0.1d0,
loc_zc=1.d-3,
dtmax=100.d0,
dtmin=1.d-10,
dt0=10.d0,
d_lum=0.1d0,
d_press=0.1d0,
d_ray=0.1d0,
d_temp=0.1d0,
age_car=0.1d0,
ini0=3,
n_atm=75,
kipp=.TRUE.,
en_m23=.FALSE.,
ctel=0.d0,
ctep=-1.d0,
ctem=15.d0,
cter=0.d0,
ctet=0.d0,
dn_fixe=0.05d0,
dpsi=0.05d0,
mu_saha=.TRUE.,
ajuste=.FALSE.,
lisse=.FALSE.,
npt_lisse=1,
q0=0.0d0,
l0=0,
new_bv=.TRUE.,
fcv=.FALSE.,
des_instb=.FALSE.,
no_discon=.FALSE.,
n_min_ZC=5,
l_fac=.TRUE.,
d_zc=5
/
\end{verbatim}\end{tt}
Les significations des paramètres sont les suivantes~:
\begin{description}
\item[m\_qs=1]~: Ordre des B-splines pour les variables quasi-statiques.
\index{m\_qs}
\item[m\_ch=2]~: Ordre des B-splines pour l'interpolation de la composition
chimique.\index{m\_ch}
\item[m\_rot=2]~: Ordre des B-splines pour l'interpolation du moment cinétique.
\index{m\_rot}
%\item[m\_rl=8]~: Ordre des B-splines pour le lissage des coefficients de la
%diffusion du moment cinétique.
\item[m\_tds=2]~: Ordre des B-splines pour interpolation de l'énergie
graviphique.\index{m\_tds}
\item[m\_ptm=2]~: Ordre des B-splines pour interpolation de la perte de masse.
\index{m\_ptm}
\item[ordre=4]~: Ordre du schéma d'intégration des réactions nucléaires.
avec {\tt rk\_imps}.\index{ordre}
\item[precic=1.d-4]~: Précision sur l'intégration temporelle de la composition
chimiqie (réactions nucléaires)\index{precic}
\item[precix=1.d-3]~: Précision sur les itérations Newton-Raphson des
intégrations spatiales.\index{precix}
\item[precit=0.15d0]~: Variation maximale/masse pour l'intégration temporelle
 de la composition chimique.\index{precit}
\item[psi0=0.1d0]~: Constante de répartition à assurer.\index{psi0}
\item[loc\_zc=5.d-3]~: Précision de la localisation des limites \ZRZCs.
\index{loc\_zc}
\item[dtmax=100.d0]~: Pas temporel maximum.\index{dtmax}
\item[dtmin=1.d-10]~: Pas temporel minimum.\index{dtmin}
\item[dt0=10.d0]~: Pas temporel initial.\index{dt0}
\item[d\_lum=0.1d0]~: Variation temporelle maximale en luminosité.\index{d\_lum}
\item[d\_press=0.1d0]~: Variation temporelle maximale en pression.\index{d\_press}
\item[d\_ray=0.1d0]~: Variation temporelle maximale en rayon.\index{d\_ray}
\item[d\_temp=0.1d0]~: Variation temporelle maximale en température.\index{d\_temp}
\item[age\_car=0.1d0]~: Estimation du temps nécessaire à l'amortissement des
conditions initiales physiquement incohérentes.\index{age\_car}
\item[ini0=3]~: Nombre d'itérations Newton-Raphson avec réestimation
de la composition chimique, de la vitesse angulaire et des limites \ZRZCs.
\index{ini0}\label{sec:kipp}
\item[n\_atm=75]~: Nombre de couches dans l'atmosphère restituée.\index{n\_atm}
\item[kipp=.TRUE.]~: Utilisation de l'approximation de Kippenhahn.\index{kipp}
\item[en\_m23=.TRUE.]~: Variables utilisées: $\ln P$.\index{en\_masse}
\item[ctel=0.d0]~: Facteur de répartition en luminosité.\index{ctel}
\item[ctep=-1.d0]~: Facteur de répartition en pression.\index{ctep}
\item[ctem=15.d0]~: Facteur de répartition. en masse.\index{ctem}
\item[cter=0.d0]: Facteur de répartition. en rayon.\index{cter}
\item[ctet=0.d0]~: Facteur de répartition. en température.\index{ctet}
\item[dn\_fixe=0.05d0]~: Taux limite de variation du nombre de couches
nécessitant un réajustement de la grille
 fixe pour l'interpolation de la composition chimique.\index{dn\_fixe}
 \item[dpsi=0.05d0]~: Taux limite de variation de la fonction de répartition
nécessitant une modification du nombre de couches.\index{dpsi}
\item[mu\_saha=.TRUE.]~: on supposera le milieu totalement ionisé pour le
calcul du poids molécualaire moyen.\index{mu\_saha}
\item[ajuste=.FALSE.]~: ajustement de T, X central, He core, Rext (sauf Teff).
\index{ajuste}
\item[lisse=.FALSE.]~: la composition chimique est lissée par diffusion
turbulente.\index{lisse}
\item[npt\_lisse=1]~: nombre de points à droite et à gauche pour le lissage des
discontinuités de composition chimique.\index{npt\_lisse}
\item[q0=0.1d0]~: Dans les fichiers ASCII de sortie, un point est placé à
$q0>0$ fois l'espacement entre les deux premiers points.\index{q0}
\item[l0=4]~: Dans les fichiers ASCII de sortie,
ajout de l0+1 points entre lim\_ZC et lim\_ZC+1, de
l0 points entre lim\_ZC+1 et lim\_ZC+2..............................
 1 point entre lim\_ZC+l0 et lim\_ZC+l0+1
et par symétrie entre lim\_ZC et lim\_ZC-l0.\index{l0}
\item[new\_bv=.TRUE.] : Calcul de la \fbv\ utilisant
$\varphi=\dsdx{\ln\rho}{\ln\mu}$.\index{new\_bv}
\item[fcv=.FALSE.]~: on ne forcera pas la convergence.\index{fcv}
\item[des\_instb=.FALSE.]~: la zone d'instabilité n'apparaitra pas dans le graphe.\index{des\_instab}
\item[no\_discon=.FALSE.]~: on tiendra compte des discontinuités.\index{no\_discon}
\item[n\_min\_ZC=5]~: nombre minimum de couches dans une ZC.\index{n\_mi\_ZC} ,
\item[l\_fac=.TRUE.]~: on rapprochera à mieux que $d\_zc\%$ chaque discontinuité
d'un bord de couche.\index{l\_fac}
\item[d\_zc=5]~: taux de rapprochement des discontinuités d'un bord de couche.
\index{d\_zc}
\end{description}
L'ajustement des constantes de répartition {\embf est délicat}. Les valeurs
par défaut sont {\tt ctel=0, ctep=-1, ctem=15, cter=0, ctet=0}. Prendre
garde que {\tt ctep} et {\tt ctet} doivent être {\embf négatives}, car elles
affectent respectivement la pression et la température, quantités
décroissantes du centre à la surface.

Comme pour les abondances initiales \cf \paragr{sec:mix}, il est possible
de caractériser
ces réglages en donnant au fichier des réglages, le nom du modèle avec
l'extension {\tt .rg}, exemple~: {\tt mon\_modele.rg}. S'il existe dans
l'environnement, ce seront les paramètres de ce fichier qui seront pris en
compte pour {\embf tous les} modèles de nom {\tt mon\_modele}. En son absence, ce
seront les quantités du fichier {\tt reglages}, si aucun de ces fichiers n'est présent
dans l'environnement le calcul est interrompu.

\section{Limites des shells source}\index{shells source}\index{HHe}
\label{sec:HHe}
Dans le fichier de données, \cfparagr{sec:nlcesam}, en ajoutant au nom du type
de dessin ({\tt nom\_des})
l'extension {\tt HHe}, exemples: ({\tt des\_mHHe, des\_mCHHe}) \cesam2k
 constitue un
fichier ASCII, de nom {\tt mon\_modele.mHHe}, contenant, pour chaque modèle,
le muméro du modèle, l'âge, les limites en masse des shells source H et He
({\tt FORMAT(i4,es22.15,2es10.3)}. Exemple:
\begin{verbatim}
   0 0.000000000000000E+00 0.000E+00 0.000E+00
   1 1.457003778976599E-05 0.000E+00 0.000E+00
   2 2.622606802157879E-05 0.000E+00 0.000E+00
   3 3.555089220702902E-05 0.000E+00 0.000E+00
   4 4.301075155538921E-05 0.000E+00 0.000E+00
   5 5.034857965187844E-05 0.000E+00 0.000E+00
   6 5.830431518424666E-05 0.000E+00 0.000E+00
   .
   .
   .
 438 1.273514666553580E+02 8.460E-01 4.876E-01
 439 1.273587701858585E+02 8.461E-01 4.916E-01
 440 1.273658006341733E+02 8.461E-01 4.983E-01
 441 1.273725838461919E+02 8.463E-01 5.049E-01
 442 1.273791284635452E+02 8.457E-01 5.088E-01
 443 1.273854522544298E+02 8.459E-01 5.154E-01
 444 1.273915631736387E+02 8.460E-01 5.193E-01
 445 1.273974741112376E+02 8.462E-01 5.257E-01
\end{verbatim}

\section{Zoom}\label{sec:zoom}\index{zoom}
Il est possible de personnaliser le dessin {\tt 'en\_masse'}
du diagramme HR en y
dessinant un rectangle cible à atteindre par le trajet évolutif, ainsi que
de modifier l'échelle du diagramme, en centrant celui-ci autour de la valeur
du point cible à des distances précisées en DeX de $\log T_{\rm eff}$ et
$\log L$. On peut aussi tracer une cible de [Fe/H] à atteindre, dans
le cadre du dessin des abondances en fonction de la masse. Pour ce faire,
coder {\tt nom\_des='zoom'} dans le fichier de données et,
dans l'environnement, placer un fichier de nom {\tt mon\_modele.zoom}, dans
lequel sont introduites les données nécessaires. Exemple d'un tel fichier
 {\tt mon\_modele.zoom}\footnote{Exemple se trouvant dans le
 sous-directory {\tt EXPLOIT}.}~:
\begin{tt}\begin{verbatim}
 &NL_DES
 teff_des=6530.,
 dteff_des=50., -30.,
 zoom_t=0.2, -0.1,
 l_des=6.844,
 dl_des=0.6, -0.4
 zoom_l=0.1, -0.2,
 fesh_des=0.05,
 dfesh_des=0.03, -0.04,
 logteff_max=4.6,
 logteff_min=3.5,
 logl_max=5.,
 logl_min=2.5
 /
\end{verbatim}\end{tt}
\begin{itemize}
 \item Il n'y a pas de personnalisation si {\tt teff\_des} et/ou
 {\tt l\_des} est négatif ou nul.
 \item Le zoom n'est effectif qu'avec {\tt zoom\_t} et {\tt zoom\_l}
  non nuls.
 \item Pour ne pas dessiner de cible en métallicité, coder
 {\tt fesh\_des=100}.
 \item Les valeurs des abscisses et ordonnées maximales \ie\
 {\tt logteff\_max (min), logl\_max (min)} {\embf prévalent}
 sur celles des
 cibles; dans l'exemple précédent, les échelles logarithmiques
 en température effective et en luminosité
 variront respectivement dans les intervalles $[4.6,3.5]$ et $[2.5,5.]$.
 \item Pour {\bf seulement} adapter les échelles du diagramme HR coder,
 par exemple~:
 \begin{tt}\begin{verbatim}
 &NL_DES
 teff_des=0.,
 dteff_des=0., 0.,
 zoom_t=0., 0.,
 l_des=0.,
 dl_des=0., 0.
 zoom_l=0., 0.,
 fesh_des=0.,
 dfesh_des=0., 0.,
 logteff_max=3.6,
 logteff_min=2.5,
 logl_max=2.,
 logl_min=-2.5
 /
 \end{verbatim}\end{tt}
 \item Pour {\bf seulement} cadrer le diagramme HR autour de la cible, coder~:
 \begin{tt}\begin{verbatim}
 &NL_DES
 teff_des=6530.,
 dteff_des=50., -30.,
 zoom_t=0.2, -0.1,
 l_des=6.844,
 dl_des=0.6, -0.4
 zoom_l=0.1, -0.2,
 fesh_des=100.,
 dfesh_des=0., 0.,
 logteff_max=0.,
 logteff_min=0.,
 logl_max=0.,
 logl_min=0.
 /
\end{verbatim}\end{tt}
\end{itemize}
La signification des variables est la suivante~:
\begin{description}
\item[teff\_des=6530.]~: Température effective du centre de la cible.
\item[dteff\_des=50., -30.]~: Ecarts en température effective.
\item[zoom\_t=0.2, -0.1]~: Extension en température effective autour de la
cible en  DeX.
\item[l\_des=6.844]~: Luminosité du centre de la cible.
\item[dl\_des=0.6, -0.4]~: Ecarts en luminosité.
\item[zoom\_l=0.1, -0.2]~: Extension en luminosité autour de la cible en  DeX.
\item[fesh\_des=0.05]~: Cible à atteindre en métallicité.
\item[dfesh\_des=0.03, -0.04]~: Ecarts de la cible en métallicité.
\item[logteff\_max=4.6]~: $\log_{\rm 10}$ de l'abscisse maximale
du diagramme HR.
\item[logteff\_min=3.5]~: $\log_{\rm 10}$ de l'abscisse minimale
du diagramme HR.
\item[logl\_max=5.]~: $\log_{\rm 10}$ de l'ordonnée maximale
du diagramme HR.
\item[logl\_min=2.5]~: $\log_{\rm 10}$ de l'ordonnée minimale
du diagramme HR.
\end{description}
Il est possible de
personnaliser davantage, en donnant au fichier {\tt zoom}
le nom attaché au modèle avec l'extension {\tt .zoom},
exemple~: {\tt mon\_modele.zoom}. Dès lors, les dispositions précédentes ne
seront effectives {\embf que pour} le modèle {\tt mon\_modele}.
Avec l'option {\tt zoom}, {\sc cesam} cherche à utiliser le fichier
{\tt mon\_modele.zoom}, en son absence, le fichier {\tt zoom} est
utilisé.

\section{Device}\label{sec:device}\index{device}
Pour adapter les dimensions des cadres et le nom du device à utiliser,
placer dans l'environnement du calcul, un fichier de nom {\tt device}, dans
lequel sont introduites les données nécessaires.
Exemple d'un tel fichier {\tt device}\footnote{Exemple se trouvant dans le
sous-directory {\tt EXPLOIT}.}~:
\begin{tt}\begin{verbatim}
&NL_DEVICE
h=6.5,
dh=1.5,
ld=8.5,
dl=2.,
xleft=1.8,
ybot=1.4,
device='/xw'
/
\end{verbatim}\end{tt}
La signification des variables est la suivante (unités en cm)~:
\begin{description}
\item[h=6.5]~: Hauteur des cadres.
\item[dh=1.5]~: Espace vertical entre les cadres.
\item[ld=8.5]~: Largeur des cadres.
\item[dl=2.]~: Espace horizontal entre les cadres.
\item[xleft=1.8]~: Distance du dessin à la limite gauche de l'écran.
\item[ybot=1.4]~: Distance du dessin à la limite inférieure de l'écran.
\item[device='/xw']~: Nom du device.
\end{description}
On peut aussi modifier directement les valeurs de ces paramètres dans
le module
{\tt mod\_donnees} du sous-directory {\tt SOURCE}.

\section{Limitation des écritures}\index{limitation des écritures}
\index{limitation des sorties}
\label{sec:blabla}\index{blabla}\index{baratine}
Des informations permettant de suivre en détail le déroulement des calculs
apparaissant "on line" sur le moniteur.
Elles peuvent s'avérer non relevantes pour une exploitation "batch" et encombrer
inutilement. \cesam2k a la possibilité d'en détourner la plus grande partie
vers des fichiers qui pourront, au besoin, lus "off line",
\cfparagr{sec:baratine}. Quand la nature cette limitation de
sorties n'est qu'occasionnelle, il suffit de disposer, dans
l'environnement du calcul, le fichier suivant de nom {\tt blabla}~:
 \begin{tt}\begin{verbatim}
&NL_BLABLA
baratine=.FALSE.
/
\end{verbatim}\end{tt}
On trouvera ce fichier dans le sous-directory {\tt EXPLOIT}.

