c*********journal des modifications/implantations dans CESAM2k/SOURCE***

	CHARACTER (len=8), PARAMETER, PUBLIC :: version='V3.12.26'

!Signification du numéro de version : Va.b.c
!a augmente si les fichiers binaires de reprise *.pms, *.rep, *.dat...
! changent même si certains restent compatibles, b et c sont alors mis à 0
!b augmente si l'un des fichiers, données (*.don), réglage, sortie
! (*.osc)... change, c est alors mis à 0
!c augmente s'il y a des modifications d'algorithmes (bug, nouvelles
!  implantations, suppressions, ...) pouvant affecter les résultats

!Journal des aménagements

! 01.11.17 ----> CESAM_V3.12.26 (svn 550)
! 1 - Ajout du commentaire optionnel dans la routine sortie.f, aménagements des
!     divers appels à sortie

! 04.10.17 ----> CESAM_V3.12.25 (svn 545)
! 1 - corrections d'un bug concernant les dérivées / ro dans etat_irwin
! 2 - aménagements de divers détails mineurs

!21.09.17 ----> CESAM_V3.12.23  (svn 553)
! 1 - Make, Make_exe, Make_expl : personnalisation des link suivant les noms des
!     machines. Introduction en export du nom de le machine OSnom

! 14.09.17 ----> CESAM_V3.12.22  (svn 552)
! 1 - cesam.f rectification d'un bug pour le calcul de la gravité,
!     incidences dans list et mod_donnees
! 2 - lim_zc simplification de l'usage de l_fac, incidences dans cesam.f

! 08.09.17 ----> CESAM_V3.12.21  (svn 551)
! 1 - Ajout d'un SAVE dans control_dpt
! 2 - divers détails dans SOURCE et EXPLOIT

! 25.08.19  ----> CESAM_V3.12.20 (svn 550)
! 1 - Dans SOURCE, dans mod_opa suppression de z14xcotrin21, dans opa.f suppression
!     des appels à opa_opal2_co et opa_opal2_cno par l'intermédiaire de opa_opal2
! 2 - Dans mod_donnees: initialisation de nom_chemin à '~/SUN_STAR_DATA'

! 29.07.17 ----> CESAM_V3.12.19 (svn 539)
! 1 - Dans SOURCE, les lois T(tau) de L.Piau gr1interpatmas59cgm.f et
!     gr3interpatmas59cgm.f sont respectivement renommées piau1.f et piau3.f
!     adaptation dans AT1PopIas59tdetau.f et AT3PopIas59tdetau.f
! 2 - lim_zc: ajouts de print idm6 à divers endroits

! 03.07.17 ----> CESAM_V3.12.18 (svn 507)
! 1 - Dans SOURCE/cesam.f, list.f csv_lb : ajout du calcul de p0 dans la fonction P0
!     et écritures. Ajout de p0 dans le HR
! 2 - Réorganisation des README dans SOURCE et EXPLOIT

! 28.06.17 ----> CESAM_V3.12.17 (svn 506)
! 1 - Dans SOURCE/ecrit_nl adjonction d'une entrée OPTIONAL permettant d'actualiser
!     la longueur de mélange, les abondances en X et Y, la masse et la vitesse de
!     rotation
! 2 - dans EXPLOIT introduction de routines de calibration

! 23.05.17 ----> CESAM_V3.12.16 (svn 505)
! 1 - Dans la SOURCE,: petites corrections dans corr_bess.f, corr_bol_lejeune,
!     etat_opal5Z. Dans EXPLOIT: petits aménagements dans des2k_dvson, Make_expl,
!     mon_modele.don. Introduction du shell-script create_evol_space.

! 13.05.17  ----> CESAM_V3.12.15 (svn 504)
! 1 - Introduction de l_pertmt calcul de le perte de masse due aux réactions
!     nucléaire. Ce calcul est effectif quand le subscript $ est mis à la fin
!     du nom de la perte de masse Ex: no_pertem$.
! 2 - Suppression du contrôle des variations de P et T quand il y a perte de masse


! 25.04.17  ----> CESAM_V3.12.14 (svn 500)
! 1 - utilisation cosdd et sindd dans fleche, COSD et SIND n'étant pas toujours
!     reconnus par gfortran

! 19.04.17  ----> CESAM_V3.12.13 (svn 495)
! 1 - Introduction dans la SOURCE des modules d'exploitation et d'outils astro
!     et des routines correspondantes, adaptation de Make, Makefile de création !     de la SOURCE
! 2 - Création du Makefile Make_expl de construction des exécutables
!     d'exploitation (calc fesh.out : calcul de Fe/H)
!

! 18.03.17  ----> CESAM_V3.12.12 (svn 495)
! 1 - Corrections du Makefile de la SOURCE (Make)

! 06.03.17  ----> CESAM_V3.12.11 (svn 492)
! 1 - Simplification du Makefile de la SOURCE  (Make)

! 24.02.17  ----> CESAM_V3.12.9 (svn 490)
!1 - Utilisation de règles implicites dans le makefile Make de la SOURCE de cesam

! 19.02.17 ----> CESAM_V3.12.4 (svn 487)
! 1 - Suppression de tous les Makefiles et de core.mak. Remplacement par le Makefile
!     Make permettant les diverses compilations avec divers compilateurs.

! 04.02.17 ----> CESAM_V3.12.3 (svn 486)
! 1 - Remplacement de HOUDEK par H12 dans les règles de links des Makefiles

! 30.01.17 ----> CESAM_V3.12.1 (svn 484)
! 1 - Remplacement de opa_houdek9 par opa_houdek12, adaptation de core.mak

! 19.01.17 ----> CESAM_V3.11.2 (svn 483)
! 1 - Introduction du core.mak dans les Makefiles
! 2 - Mise en place de opa_houdek12.f

! 07.01.17 ----> CESAM_V3.11.1(svn 482)
! 1 - Amémagements / créations des Makefiles: Make2k Make2k-dbg Makeifort
!     Makeifort-dbg MakeT MakeT-dbg

! 21.09.16 ----> CESAM_V3.11.0 (svn 479) simplifications + nettoyage
! 1 - modification du fichier de réglage
! 2 - suppression de variables l_demi, avec_tds
! 3 - suppression de routines actu_esp, Lnuc_tds,
! 4 - suppression de la luminosité dans le contrôle de la répartition
! 5 - contrôle du déplacement des limites ZR/ZC vers un bord de couche
! 6 - valeurs des constantes et de la fonction de répartitions fixées

! 08.12.15 -----> CESAM_V3.10.9 (svn 478)
! 1 - Divers détails dans Z_read2_saha et evol
! 2 - Implémentation du package etat_opal5 procuré par Anna

! 01.12.15 -----> CESAM_V3.10.8 (svn 475)
! 1 - Implémentation du package etat_saha.f selon Gryaznov etal. procuré par Anna & Baturin.

! 07.11.15 -----> CESAM_V3.10.7 (svn 473)
! 1 - cesam.f, resout.f, des_m.f, des_r.f : augmentation du nombre max de modèles,
!     modification du format de l'écriture de model_num i4 --> i5
! 2   Aménagements mineurs dans rq_reac.f et adjonction d'une sortie OPTIONAL ecranw
! 3   Adjonction de la routine control_dPT.f du module mod_static: contrôle des
!     variations temporelles de P et T

! 12.09.15 -----> CESAM_V3.10.6 (svn 472)
! 1 - mod_mumerique: suppression de a_eq_b_*. a_ne_b_*.f, eiggen.f
! 2 - mod_numerique: ajout de v_propres.f, suppression de eggen.f
! 3 - mod_variables, inter.f: arrêt de la dichotomie à epsi=1.d-6
! 4 - mod_nuc, rq_reac.f : introduction de l'écran fort
! 5 - mod_static: mise à jour de static_m23.f
! 6 - divers endroits: corrections mineures

! 03.09.15 -----> CESAM_V3.10.5 (svn 471)
! 1 - Ajout de svn_cov.f dans mod_numerique (oubli pour V3.10.4)

! 03.09.15 -----> CESAM_V3.10.4 (svn 470)
! 1 - cesam.f: kipp=.TRUE. utilisation de la simplification de Kippenhahn par défaut
! 2 - corrections mineures dans osc_nadia, svd_cmp, static_m13

! 22.08.15-----> CESAM_V3.10.3 (svn 469)
! 1 - ajout dans mod_numerique, svd_cmp.f : recherche des valeurs singulières,
!     eiggen .f : calcul des valeurs propres. des module procedure a_eq_b, a_ne_b
! 2 - abon_ini.f : amélioration de l'algorithme de personnalisation des abondances initiales
! 3 - inter.f : modification du critère d'arrêt de la dichotomie
! 4 - ajout dans mod_numerique a_eq_b.f , a_ne_b.f , simq.f , maehly.f , svdcmp.f

! 27.07.15----> CESAM_V3.10.2 (svn 468)
! 1 - cesam.f : suppression de la garde pour arrêt sur l'âge

! 11.07.15 ----> CESAM_V3.10.1 (svn 467)
! Nombreux aménagements motivés par la tentative de franchir le flash de l'hélium
! pour le soleil. Le fichier des réglages a été aménégé.

! 08.02.15 ----> CESAM_V3.9.2 (svn 466)
! 1 - Création des MAKEFILEs Make2k, Make2k-dbg, MakeT, MakeT-dbg pour compilation lf95,
!     Makeifort, Makeifort-dbg pour compilation ifort
! 2 - Introduction de control_LPRT: limitation des variations temporelles des MAX
! 3 - Transformation de l'INCLUDE tabul_TdS.f en sous-programme interne de resout.f (une
!     correction dans un INCLUDE n'est pas vue par Make)
! 4 - lim_zc.f les limites ZR/ZC trop externes (Mzc/Mtot < 1.d-7) ne sont plus affinées,
! 5 - Les ZC trop externes (Mzc/Mtot < 1.d-7) ne sont plus mélangées
! 6 - Nombreuses corrections othographiques UTF-8 --> 8859-1

! 01.12.14 ----> CESAM_V3.9.1 (svn 462)
! 1 - Addition de |e_grav| dans Lnuc.f ---> Lnuc_tds.f permettant de tenir compte de
!     l'énergie graviphique dans la fonction de répartition. La luminosité
!     "gravo-nucléaire" ainsi obtenue, toujours croissante, est utilisée dans la
!     fonction de répartition. Cette disposition est
!     introduite par un des suffixes G ou N de l'argument "PRECISION",
!     Ex:PRECISION='PRG'; G(gravité+nucléaire) ou N(nucléaire seul)
! 2 - Adaptation de la routine espmt.f: adaptation de ctel en fonction des dispositions
!     précédentes
! 3 - Les itérations d'initialisation de ZAMS sans énergie graviphique (dt=0) sont
!     poursuivies jusqu'à obtention d'un nombre de couches compatible avec la
!     répartition.
! 4 - Introduction de la variable LOGICAL l_actu définissant les coefficients de
!     répartition à actualiser dans la routine actu_esp.f. Par défaut,
!     les coefficients de répartition pour la pression et la masse sont respectivement
!     fixés à -1 et 15. Avec les options "N" ou "G" le coefficient de la luminosité
!     est variable et fixé à l'inverse du maximum de la luminosité "gravo-nucléaire".
!     Les coefficients de répartition sont nuls pour la température et le rayon.
! 5 - Le nombre maximal de couches est fixé dans le fichier de données dans la
!     limite de 5000. Lorsque ce nombre maximal de couches est atteint, la valeur du
!     coefficient de répartition utilisée peut dépasser sa valeur initiale fixée à 0.03
! 6 - Remplacement de la routine coeff_esp.f par actu_esp.f. Suppression des routines
!     control_PT.f et control_LPRT.f
! 7 - Pour les modèles avancés, création des précisions av, ce, He, utilisant les
!     variables ln P, ln T, R^2, L^2/3 et M^2/3
! 8 - Introduction du tableau logique l_actu permettant de choisir les coefficients de
!     répartition actualisables au cours de l'évolution; addition d'une ligne dans le
!     fichier des réglages.
! 9 - Création d'un fichier de sortie réduit en lieu et place du fichier HR. Obtenu
!     ajoutant le suffixe "R" à l'argument NOM_OUTPUT du fichier de données
!     Ex: NOM_OUTPUT=osc_adiaR; extension *.HR --> *.LR Ex: soleil.LR. Création de la
!     routine d'exloitation lit_lr.f et du programme de dessins des2k_lr.f

! 08/07/14 ----> CESAM_V3.8.5 (svn 460+1)
! Principale implémentation: contrôle de la répartition en tenant de la "luminosité
! nucléaire": intégrale sur l'* du débit epsilon_nuc d'énergie nucléaire
! 1 - evol.f déplacement des masses mc_tmp au plus près des limites convd/convf des ZC
! 2 - création de la routine PUBLIC de mod_evol calc_grad : calcul des gradients
! 3 - création du programme  des2k_Lrad : dessin des luminosités en utilisant les données
!     d'un fichier binaire (directory EXPLOIT)
! 4 - lim_zc.f, resout.f,  mod_static.f: introduction de l'entier PRIVATE rap_max du module
!     mod_static : maximum des rapports de localisation des limites ZR/ZC. Dans resout.f la
!     localisation des limites sera poursuivie tant que celles-ci ne seront pas à moins
!     de 5% d'un point de grille.
! 5 - coll_qs.f : amélioration de la gestion des itérations et des corrections NR
! 6 - diverses redéfinition de ini0, ini1, iter_max dans cesam.f et resout.f
! 7 - Création de la FUNCTION PUBLIC Lnuc.f des module mod_static et mod_exploit: calcul de
!     l'intégrale de l'énergie nucléaire du centre à un point de l'espace
! 8 - Création de la routine PUBLIC coeff_esp.f du module mod_static : i) actualisation des
!     coefficients d'espacement, ii) du facteur de répartition en fonction de la
!     variation relative du maximum de luminosité. Utilisation de Lnuc pour le contrôle
!     de la répartition
! 9 - création de la routine actu_nqs.f: ajustement du nombre de couches
! 10 - création des routines ecrit_static.f et plot_static.f dessins pour debug. Adjonction
!     du flag de contrôle Kdes_stat défini dans lit_nl.f
! 11 - rectification d'un bug (,,) dans coll_qs.f
! 12 - rectifications du codage de lettres accentuées (é, è, ê, à etc...)

! 11/05/14----> CESAM_V3.8.4 (svn 457)
! 1 - cesam.f : par défaut, dessin de la bande d'instabilité avec PRECISION= CE, ce
! 2 - lim_zc.f: reprise du calcul du facteur de répartition fac en utilisant les
!     coefficients d'espacement ce qui s'est avèrè plus efficace que leur détermination
!     à partir des positions des points de grilles

! 20/04/14----> CESAM2k_V3.8.4 (svn 456)
! 1 - des_m.f, des_r.f : écriture du nombre de couches dans le plot
! 2 - resout.f : si err. max. est en 1, grad_rad=0 dans les écritures
! 3 - codage UTF8 --> 8859-1 dans diverses routines
! 4 - suppression de la routine inutilisée left_right.f du module mod_numerique

! 01/02/14----> CESAM2k_V3.8.3 (svn 455)
! 1 - mod_numerique: implantation de la routine générique shell
! 2 - coll_qs.f : aménagements et réduction des ALLOCATEs
! 3 - control_LPRT.f : contrôle des variations de T pour toutes les couches

! 06/12/13----> CESAM2k_V3.8.2 (svn 454)
! 1 - resout.f encadrement du nombre de couches par [n_min,n_max]
! 2 - resout.f iter_max0=100 pour les initialisations PMS et ZAMS

! 03/12/13----> CESAM2k_V3.8.1 (svn 453)
! 1 - cesam.f: correction de bugs dus à une erreur d'utilisation de espm.f
! 2 - mod_static.f, espmt.f devient PRIVATE dans mod_static
! 3 - gr3interpatmas59cgm.f, tho.f, gr1interpatmas59cgm.f, z14xcotrin21.f,
!     trho_4000.f, opa_opalCO.f corrections D+3 pour les formats F

! 01/12/13----> CESAM2k_V3.8.0 (svn 450)
! 0 - les items d_grav et ro_test ont été supprimés dans le fichier de réglages
! 1 - list.f: remplacement de TdS par Eg=-(dU+PdV)/dt; indication sur la première
!     colonne de chaque couche de sa nature radiative ou convective et si elle
!     est en compression ou en extension
! 2 - resout.f: implantation des INCLUDEs ajustements.f et calc_TdS.f
! 3 - lim_zc.f: simplification du calcul du facteur de répartition fac. Elimination
!     de la dichotomie pour les limites ZC/ZC trop externes
! 4 - Améliorations de l'algorithme de dichotomie de inter.f
! 5 - suppression de d_grav et ro_test, faisant double emploi avec le contrôle
!     sur les variations temporelles de P et T.

! 15/06/13----> CESAM2k_V3.7.34 (svn 449)
! 1 - lit_nl.f suppression de la perte de masse interne due à mc^2 avec pertm_ext
!     et mdot=0
! 2 - evol.f suppression de la possibilité d'abondances < 1.d-50
! 2 - Adaptations mineures dans coll_qs.f, mod_static

! 22/04/13----> CESAM2k_V3.7.33 (svn 448)
! 1 - lim_zc.f : rectification d'un bug lié au déplacement des limites ZR/ZC
! 2 - Divers ménagements mineurs

! 12/04/13----> CESAM2k_V3.7.32 (svn 447)
! 1 - static_m13.f : déplacement des définitions de dtn et dpn
! 2 - cesam.f : augmentation de la valeur par défaut de psi0=0.05 --> 0.07

! 01/04/13----> CESAM2k_V3.7.31 (svn 446) (svn 445 a été contournée)
! 1 - static_m13.f, resout.f : suppression de la forme 1 du TdS de Kippenhahan
! 2 - lim_zc.f amélioration de la gestion des abscisses réelles new(jlim)
! 3 - resout.f, coll_qs.f, reout.f : suppression de iter_stag
! 4 - coll_qs.f : aménagements dans le calcul de l'amortissement NR et pour
!     corrections en LnP et LnT

! 02/03/13----> CESAM2k_V3.7.29 (svn 444)
! 1 - cesam.f: suppression de la precision He qui n'apporte rien de significatif.
! 2 - aménagements secondaires divers
! 3 - paramètres pour les stades avancés: L, TdS, T n'étant pas toujours des
!     des fonctions monotones, mises à 0 de ctel, ctet, cter étant laissés à 0

! 13/01/13 ----> CESAM2k_V3.7.28 (svn 442)
! 1 - list.f: affichage du gamma de la convection en lieu et place de cv,
!     incidences dans cesam.f et ecrit_ascii.f
! 2 - lit-nl.f, lim-zc.f, mod_donnees.f : introduction du suffixe 'S' pour les
!     arguments de la NAMELIST NOM_CONV='conv_cm_rezaS' il y suppression des
!     placements de couche sur chaque limite ZR/ZC

! 21/12/12 ----> CESAM2k_V3.7.26 (svn 441)
! 1 - cesam.f: suppression d'écritures de debug après la formation du fichier.rep
! 2 - corrections mineures d'écritures

! 20/12/12 ----> CESAM2k_V3.7.25 (svn 437) (EXPLOIT + NOTICE--> 439)
! 1 - control_LPRT.f : contrôle des variations de T limité aux 50 premières
!     couches
! 2 - lit_binaire.f : ajout de lim et dts dans la liste de lecture
! 3 - routines de diffusion microscopique: introduction d'un coefficient de
!     diffusion turbulente dans les ZC égal à 1/3 lm Vc. Effectif en affectant
!     le subscript V au nom de la routine de diffusion turbulente Ex: difft_nuV
!     introduction dans lit_nl.f, adaptation dans mod_exploit.f, mod_donnees;
!     introduction de d_cv dans mod_evol.f, evol.f, eq_diff_chim.f
!     création de d_convec.f : calcul du coefficient de diffusion
! 4 - changement d'unité du temps caractéristique de convection jours --> secondes

! 26/11/12 ----> CESAM2k_V3.7.24 (svn 436)
! 1 - lit_nl.f, cesam.f : corrections diverses relatives aux sorties de
!     modèles ASCII et/ou binaires intermédiaires

! 24/11/12 ----> CESAM2k_V3.7.23 (svn 435)
! 1 - resout.f, csv_tr.f, mod_static.f, mod_variables.f, t_convec.f:
!     aménagements pour les sorties "TR" (S_entro, S_0atm, tot_conv,
!     MI_atm, M_iner), déplacement des allocations de cesam-->resout

! 14/11/12 ----> CESAM2k_V3.7.22 (svn 434)
! 1 - resout.f: calcul des moments d'inertie entre les limites ZC/ZR
! 2 - cesam.f: ajout d'une virgule après les entropies dans l'en-tête
!     des TR csv; utilisation de kipp=.TRUE. avec la précision PR

! 04/11/12 ----> CESAM2k_V3.7.21 (svn 433)
! 1 - création de la routine t_convec.f: calcul du temps caractéristique
!     de convection d'une ZC (en jours)
! 2 - resout.f, cesam.f: aménagement du contenu du fichier
!     mon_modele_TR.csv, ajout des moments d'inertie aux limites ZR/ZC,
!     des temps caractéristiques de convection des ZC, de l'entropie
!     spécifique aux limites
! 3 - création des routines calc_entropie.f: calcul de l'entropie en
!     utilisant l'EOS etat_irwin dont un argument de sortie a été ajouté
!     (disposition qui n'impose pas nécessairement de coder
!     NOM_ETAT='etat_irwin' dans le fichier de données mon_modele.don)
! 4 - rconv_mean.f : calcul du rayon situé à 1/2 longueur de mélange de
!     la base d'une ZC
! 5 - Energie graviphique: remplacement de la formule 4.25  par 4.27 de
!     Kippenhahn & Weigert, incidences dans resout.f, static_m13.f,
!     inchangé dans static_m23.f pour cohérence avec les anciennes versions

! 27/10/12 ----> CESAM2k_V3.7.20 (svn 431)
! 1 - cesam.f: adaptation de paramètres pour la précision He

! 22/09/12 ----> CESAM2k_V3.7.19 (svn 430)
! 1 - pertm.f : introduction de la perte de masse 'pertm_dominik':
!     "stationary dust-driven winds" Dominik & al. 1990 A&A, 240, 365
!     et de 'pertm_reim-domi' si la perte dominik > reimers
! 2 - tabul_nuc.f, taux_nuc.f : implémentation du taux du 3alpha de
!     Nguyen, Numes & Thomson, astro-ph arXiv: 1209.44999v1 22/09/12
!     appelé par NACRE& ou NACRU&
! 3 - resout.f : en cas de difficulté, dt--> 0.75 dt au lieu de dt/2
! 4 - lim_zc.f : simplification du calcul des corrections fac(i) permettant
!     de déplacer les limites ZR/ZC sur un point de discrétsation en masse
!     Adaptations dans static_m13.f et static_m23.f

! 19/09/12 ----> CESAM2k_V3.7.18 (svn 429)
! 1 - cesam.f: calcul de dr/dt en km/s, list.f: écriture de dr/dt en lieu
!     et place de Pgaz/Ptot, écriture de dRext/dt. Adaptation des
!     paramètres pour la précision He.
! 2 - evol.f correction d'un bug pour l'extraction de T pour le calcul
!     de la dégénérescence, limitation du pas temporel (10 ans) à
!     l'approche du flash de l'hélium en cas de dégénérescence
! 3 - resout.f: amélioration des écritures des solutions intermédiaires
!     pour debug

! 30/08/12 ----> CESAM2k_V3.7.17 (svn 426)
! 1 - lim_zc: reprise de l'ancienne gestion de fac(i)
! 2 - Ajustement des réglages pour les stades avancés avec tentative,
!     pour 1Msol PopI, d'atteindre l'amorçage du 3alpha (succès) et de le
!     dépasser (échec); pour utilisation coder:  precision='He'

! 20/08/12 ----> CESAM2k_V3.7.16 (svn 425)
! 1 - Introduction du TdS dans la fonction de répartition
! 2 - Ajustements divers pour les stades avancés
! 3 - lim_zc: réaménagement de la répartition quand l'extension, fac(i), de
!     la couche i est extérieure à [ 2/3 , 2 ]

! 18/08/12 ----> CESAM2k_V3.7.15 (commit S3)
! 1 - lit_binaire.f: possibilité d'utiliser la physique du *.don qui,
!     éventuellement, peut différer de celle utilisée pour le calcul du modèle
!     dans le but de tester l'effet d'une autre physique (Ex. avec des2k_bin)
! 2 - Réaménagements du module mod_numérique
! 3 - lim_zc.f: redistribution des couches si fac en dehors de [0.3,1.7]
! 4 - lit_nl.f, cesam.f, mod_donnees.f: introduction de la possibilité de
!     former des fichiers intermédiaires d'oscillations et/ou de reprise
! 5 - Perte de masse: rassemblement des formules de perte de masse dans
!     la routine générique pertm. Suppression des routines de perte de masse.
!     S'il y a perte de masse, prise en compte systématique de la
!     perte locale de masse due à e=mc^2.
!     Avec nom_pertm = no_pertm il n'y a aucune perte de masse.
! 6 - Introduction du TdS dans la fonction de répartition avec, par défaut,
!     ctel=0.05 et cter=0.05
! 7 - Ajustements divers pour les stades avancés


! 14/07/12 ----> CESAM2k_V3.7.14 (commit S2)
! 1 - Création de la routine PUBLIC dege_elec.f du module mod_numerique
! 2 - resout.f et evol.f: écriture de la valeur de la dégénérescence
!     électronique et utilisation de sa variation temporelle maximale
!     pour contrôle du pas temporel
! 3 - tabulation vth: ajout de la valeur de la dégénérescence
! 4 - resout.f & lim_zc.f remplacement de no_rep par new_rep pour
!     le réaménagement de la répartition
! 5 - aménagements mineurs dans cesam.f, update.f resout.f, lmix.f,
!     lim_zc.f, dgrad.f
! 6 - Fixation de la constante de répartition psi0=0.06, les variations
!     dues à l'évolution étant compensées par l'augmentation/diminution
!     du nombre de couches
! 7 - dnunl.f: diminution de l'ordre de la spline d'interpolation
! 8 - static_m13.f : on tente de converger avec r et m^1/3 < 0
! 9 - cesam.f: réaménagement de ini0 (nb. iter. NR avec réestim. de lim.
!     ZR/ZC et comp. chim) ce qui permet d'atteindre 0% pour la localisation
!     des limites ZR/ZC et assurer une meilleure convergence
! 10 -resout.f: ajustement de la dépendance de iter_max en fonction de ini0

! 25/06/12 ----> CESAM2k_V3.7.13 (commit S2)
! 1 - resout_chim.f: SAVE mg(ordre des intégrations GAUSS)
! 2 - cesam.f: précision AV, suppression du lissage, psi0 à 0.035
! 3 - introduction de la routine de perte de masse pertm_reimers.f
! 4 - control_LPRT.f: suppression du contrôle des variations temporelles
!     de pression s'il y a perte de masse (Pb. aux couches externes)
!     aménagements dans pertm_ext, mod_static, static_m13, pert_msol
! 5 - pertm_tot.f: correction dans le cas m^1/3

! 19/06/12 ---> CESAM2k_V3.7.12 (commit S4)
! 1 - Rectification d'un bug dans lim_atm.f teff et grav pouvant ne pas
!     être updatés correctement
! 2 - cesam.f: Suppression d'un appel inutile à la restitution de l'atmosphère
! 3 - cesam.f, resout.f: modif. de l'utilisation de la convergence forcée (fcv)
!     Réduction à 10 du nombre max. d'itérations pour les évolutions
! 4 - cond_mestel.f correction de la formulation
! 5 - evol.f réaménagements dans l'estimation de la précision
! 6 - rk_imps.f : rectification des commentaires de sortie en cas de variation
!     trop grandes de T ou de ro

! 02/06/12 ---> CESAM2k_V3.7.11 (commit S1)
! 1 - Création de la routine csv_lb.f pour sortie csv selon des desiderata
!     de L Bigot Implications dans lit_nl, mod_cesam, mod_exploit, cesam
! 2 - Création de la routine write_rg.f: constructeur de fichiers de réglage *.rg
! 3 - Menus aménagements dans l_mix, l_mixt, csv.tr, resout

! 19/05/12 --> CESAM2k_V3.7.10 (commit S1)
! 1 - Création des routines ecrit1, ecrit1 d'écritures pour condenser cesam.f
! 2 - A chaque itération, calcul des variables destinées aux sorties ASCII
!     à chaque pas temporel, en particulier nu0, deltanu0.
! 3 - Création de la routine csv_tr.f en remplacement de vais-tr.f
! 4 - Suppression de gamma_atm de la liste d'arguments de list.f
! 5 - Rectification de bugs dans vais_tr, suppression des arguments.
!     Introduction de lstar dans mod_variable, utilisation dans atm, vais_tr

! 09/05/12 --> CESAM2k_V3.7.9 (commit S3)
! 1 - Implantation des fichiers csv*.f du module mod_numerique pour les écritures
!     ASCII facilement lisibles par les logiciels du type IDL (à l'initiative de
!     T.Guillot)
! 2 - Implantation de la routine csv_tr.f pour écritures personnalisées
!     des fichiers ASCII *.TR et aménagements dans cesam.f
! 3 - Modification du libellé de l'extension *.TR  --> *_TR.csv
! 4 - Suppression de l'écriture de l'en-tête par csv_write
! 5 - Normalisations de rv_zc et mv_zc
! 6 - Réactualisation de r_zc et m_zc dans evol.f à l'issue de la mise en place des
!     discontinuités de comp. chim.
! 7 - Dans k5750 et K5777 recadrage de l'input tau

! 22/04/12 --> CESAM2k_V3.7.8 (commit S4)
! 1 - Mise en place du l'INCLUDE vais-tr.f estimant numériquement dv_tr(lim)
!     dérivées / ln r de la fréquence de Brunt-Väissälä aux limites ZR/ZC
!     quantités écrites dans le fichier mon_modele.HR
!     incidences dans cesam.f, lim_zc.f, mod_exploit.f, mod_donnees.f,
!     mod_variables.f, write_nl.f, lit_hr.f
! 2 - Création d'une ligne spéciale pour Tristan dans le *.HR;
!     après les abondances ligne des oscillations nu0,dnu02,dnu13,a
!     éventuellemnt (subscript T dans nom_output) ligne pour Tristan:
!     mmtI,c_tr=vais(2),pt(1),t(1),ro(1),dv_tr(1:lim),ro_tr(1:lim)
! 3 - Création d'un fichier de sortie spécial pour Tristan mon_modele.TR
!     formé et écrit dans l'INCLUDE vais_tr
! 4 - Corrections de bugs dans vais_tr

! 05/04/12 --> CESAM2k_V3.7.7 (commit S1)
! 1 - Mise en place de dn_sort : avec, dans le fichier de données *.don,
!     dtlist < 0,  Ex: dtlist=-50.d0,  le fichier output sera formé tous
!     les 50 modèles, incidences dans lit_nl.f, mod_cesam.f, cesam.f,
!     ecrit_ascii.f, osc_adia.f, osc_invers.f, osc_nadia.f, ascii.f
! 2 - Ajout de Pc, Tc, Roc dans le fichier HR
! 3 - write_nl.f, l130 : suppression de l'écriture sur l'unité 2

! 24/03/12 --> CESAM2k_V3.7.6 (commit S0)
! 1 - Dans resout.f : augmentation, au delà de ini0, du nombre d'itérations
!     avec lim_zc, evol, pertm si le rayon total varie de plus de 5d-3
!     entre 2 itérations
! 2 - Dans resout.f : amélioration de la logique pour arrêt quand
!     Tstop ou Rstop est traversé

! 14/03/12  --> CESAM2k_V3.7.5 (commit S5)
! 1 - Sortie du dernier modèle avec |-n_max| couches lorsqu'il y a arrêt
!     avec Tc, Rext, Xc, fixé pour tous les types de précision, éventuellement
!     avec ajustement de la localisation
! 2 - Introduction dans cesam.f du C de Tristan (c_tr) dans le fichier HR et lit_hr.f
! 3 - Le contrôle des variations LPRT est effectué que si l'état de mélange
!     est identique au point de contrôle aus temps t et t+dt. Idem pour les
!     variations du TdS dans static_m13 et static_m23. Introduction de lmix_t,
!     x_mix_t, n_mix_t
! 4 - Mise en place dans lit_nl.f d'un critère pour rendre impossible la
!     présence simultanée de plusieurs critère d'arrêt avec n_max < 0.
! 5 - resout.f: mise en place de critères et d'écritures pour permettre
!     le reprise d'une évolution avec n_max < 0 et d'arrêts lorsque, dans
!     le modèle repris, le critère d'arrêt est déjà dépassé.
! 6 - correction dans lit_nl d'un bug lié à new_bv

! 16/01/12  --> CESAM2k_V3.7.4 (commit S2)
! 1 - SUN__DATA/HOUDEK/lib : réécriture des COMMONs a, b, c, d, e, f
!     dans itoh_coeff_cf.f dont la forme d'origine était refusée par le
!     compilateur pathf95
! 2 - correction d'un bug de lecture f7.3 --> f6.3 dans alex05bext.f
! 3 - limitation des écritures de sortie de table dans dopintc, dopintf
!     et dopints
! 4 - ajout de bsp2dn.f (interpolation 2D) dans mod_numerique
! 5 - Introduction du suffixe A dans NOM_OUTPUT, imposant l'utilisation
!     de l'ancienne formulation de la fréquence de Brunt-Vaissala.
!     Implications dans CESAM.f et lit_nl. L'ancienne formulation reste
!     utilisée avec les précisions sa et SA
! 6 - Introductions: teff_sol, mbol_sol, init_ctes_phys dans mod_donnees,
!     f_logf, bsp2dn dans mod_numerique

! 25/12/11  --> CESAM2k_V3.7.3 (commit S1)
! 1 - resout.f: aménagement du contrôle du nombre de couches pour introduire la
!     possiblité d'augmenter/diminuer le nombre de couches en utilisant
!     l'extension M et la valeur de n_max lors des reprises.
! 2 - Mise en sommeil de la pression turbulente et de la diffusion du moment
!     cinétique
! 3 - lit_binaire.f: suppression de l'affectation de l'unité 30 au fichier de
!     données

! 17/12/11  --> CESAM2k_V3.7.1 (commit S2)
! 1 - Simplification et réorganisation des types de précision: les types de
!     précision en minuscules (Ex: pr) utilisent les variables quasi-statiques
!     ln P, ln T, (R/R_sol)^2/3, (L/L_sol)^2/3, (M/M_sol)^2/3, les types de
!     précision en majuscules (Ex: PR) utilisent les variables quasi-statiques
!     ln P, ln T, R/R_sol,  L/L_sol, (M/M_sol)^1/3.
!     mod_donnees.f: extension de precision'char*6), cesam.f: introduction de
!     suffix(4) pour différencier les extensions A, L, M, O
! 2 - suppression de l_piau et grad_piau

! 18/10/11 --> CESAM2k_V3.7.0 (commit S6)
! 1 - Ecriture dans le fichier HR de  nu0, dnu02, dnu13, a (demande de L.Bigot)
!     Passage de ces données dans les arguments du SSP list.f
!     Aménagements dans les programmes d'exploitation de EXPLOIT
! 2 - remplacement de lisse_sum.f par lisse_sumd.f limitant le lissage au
!     voisinage des limites ZR/ZC
! 3 - lim_zc.f : Initialisations de r_ov, r_zc, m_zc à -1d10, au lieu de -HUGE(1.d0)
! 4 - resout.f : formation du moment d'inertie mmtI, variable PUBLIC de mod_static.f
!     cesam.f : écriture de mmtI à la suite des variables de sismologie dans le
!     fichier HR, lit_hr.f : lecture de mmtI
! 5 - mod_cesam : déclaration PUBLIC de des_instab.f
! 6 - pp3.f: ligne 206 correction d'un <cr>

! 27/09/11 --> CESAM2k_V3.6.2 (commit S2)
! 1 - Aménagements dans opa_compton et opa
! 2 - Améliorations des paramètres pour les stades avancés
! 3 - lim_zc : avec overshoot, on s'assure que deux limites successives ZR/ZC
!     sont séparées d'au moins n_min_ZC couches

! 25/07/11 --> CESAM2k_V3.6.1 (commit S1)
! 1 - Aménagements dans saha.f : optimisations, sortie (optional) de l'énergie
!     d'ionisation totale
! 2 - taux_nuc.f : Introduction du taux de Cacioli et al. arXiv:1107.4514v1 22/07/11
!     pour la réaction n°11 N15(p,g)O16 en tant que "NACRU" (NACRE updatée)

! 25/06/11 --> CESAM2k_V3.6.0 (commit S1), fichier de réglages modifié
! 1 - Remplacement+simplifications de la routine control_LPT.f par control_LPRT.f
!     Introduction/adaptation dans cesam.f des variables de controle des variations
!     de Lext, P, T, Rext introduction dans le fichier des réglages.
! 2 - Mise en place de no_discon permettant de supprimer les discontinuités de
!     composition chimique. Introduction dans le fichier de réglage de no_discon,
!     d_ray et n_min_ZC
! 3 - Aménagements et regroupemens des précisions ar & av, cr & ce

! 23/05/11 --> CESAM2k_V3.5.4 (commit S2)
! 1 - Mise en place de nouveaux isotopes
! 2 - Implantation de nouvelles réactions nucléaires pour les stades avancés
! 3 - Introduction des neutrons
! 4 - Création de la routine ppcno3an.f cycles PP, CNO, 3 alpha, C, Ne avec neutrons

! 26/04/11 --> CESAM2k_V3.5.3 (commit S2)
! 1 - Directory EXPLOIT: dessin en indice de couche dans des2k_bin, divers ajouts
! 2 - Aménagements des précisions av et ar
! 3 - Directory EXPLOIT: création de la routine test_nuc_bin: test du jacobien du
!     réseau nucléaire d'un modèle repris en binaire
! 4 - Dans static.f: contrôle du pas temporel sur L, P, T, introduction de d_lum,
!     d_press, d_temp, supression des contrôles sur Tc, T_eff qui deviennent caduques

! 22/01/11 --> CESAM2k_V3.5.2 (commit S8)
! 1 - Correction m_qs --> ord_qs dans pertm_tot.f
! 2 - Aménagements du calcul de précision dans resout_rota.f
! 3 - Amélioration de l'utilisation de inter.f: remplacement de m23 et r2 par
!     m_stat et r_stat, m23_ou_r2 devient m_ou_r (en a1), en_masse devient en_m23,
!     static_m et static_r deviennent respectivement static_m23 et static_m13
! 4 - Nettoyages divers, suppressions de tableaux obsolètes, mise à jour des
!     routines d'exploitation.
! 5 - coll_atm.f: précision de la restitution de l'atmosphère >= 1d-9
! 6 - suppression momentannée(?) de la diffusion du moment cinétique. Aménagements
!     pour la conservation locale du moment cinétique.
! 7 - Mise à jour de routines d'exploitation

! 14/01/11 --> CESAM2k_V3.5.1 (commit S1)
! 1 - Suppression du contrôle de la variation d'abondance si le point au temps
!     t+dt se trouve dans une ZC. Suppression de lmix_t
! 2 - Formation de R2 et M23 dans coll_qs, suppressions dans resout.
! 3 - Suppression de la routine inutile integrales.f

! 07/12/10 --> CESAM2k_V3.5.0 (commit S2)
! 1 - Création et introduction de la routine lisse_sum.f: lissage par moyenne
!     glissante intégrée sur 2*npt_lisse points. Introduction dans evol.f
! 2 - Suppression de no_discon et du lissage par diffusion. Incidences dans evol.f
! 3 - Modification du fichier des réglages
! 4 - Adaptations diverses dans cesam.f, evol.f, base_chim.f, mod_données
!     ce lissage n'est pas encore implémenté dans la rotation.

! 07/10/10 --> 28/11/10 CESAM2k_V3.4.2 (commit S6)
! 1 - Implantation de la loi T(tau) 3D de L.Piau + ajustements + aménagements
! 2 - Création et aménagements de la précision 'ci' Cépheides avec r, l, m^1/3
! 3 - Création de list_progr_exploitation.f: liste des programmes d'exploitation
! 4 - Création de la logical function lmix_t(nu)=.TRUE. si nu est dans une ZC au
!     temps t utilisé dans evol.f pour l'estimation de la précision aux points
!     sautant de ZR<-à->ZC entre t et t+dt

! 07/10/10 CESAM2k_V3.4.1 (commit S2)
! 1 - Implantation de la loi T(tau) 3D de L.Piau

! 17/07/10 CESAM2k_V3.4.0 (commit S1)
! 1 - Début d'implantation des variables principales ln P, ln T, r, l, m^1/3
!     permettant d'avoir l < 0. Il n'y a pas d'oscillation pour L quand r --> 0.
!     Cet ensemble de variables remplace "en R" : en_masse=.FALSE. Les équations
!     de structure sont calculées dans static_r.f. Cette version n'est pas
!     totalement opérationnelle.
! 2 - Corrections de bugs, aménagements divers pour static_r (m^1/3, version non
!     opérationnelle).

! 19/05/10 CESAM2k_V3.3.7 (commit S3)
! 1 - Implantation du module mod_bp_for_alecian_new pour le calcul des
!     accélérations radiatives, nouvelle formulation. Des données à disposer
!     dans le sous directory SUN_STAR_DATA sont disponibles sur demande auprès
!     de B. Pichon.
! 2 - Suppression des routines alecian1 et alecian2 et suppression de kap et
!     dkapx dans l'appel à frad
! 3 - Correction de bugs liés aux dimensions de g_rad et dg_rad dans
!     diffm_br
! 4 - resout_chim.f : Détermination du nombre de points pour l'intégration de
!     Gauss en fonction de m_ch
! 5 - saha: calcul des énergies d'ionisation (OPTIONAL)
! 6 - Implantation des accélérations radiatives dans diffm_br, mise à jour
!     de la notice
!
! 27/03/10 CESAM2k_V3.3.6 (commit S3)
! 1 - Déplacement du contrôle du pas temporel en ce qui concerne Teff, Tc et L
!     de evol à static
! 2 - Ecriture, dans mon_modele.lis, de C/N en masse et en nombre à surface
! 3 - Mise en place des routines T(tau) piau1 piau3 de L.Piau, les tables
!     pour leur exploitation sont disponibles sur demande auprès de
!     laurent.piau@cea.fr.

!02/02/10 CESAM2k_V3.3.5 (commit S12)
! 1 - Création de la routine ppcno3a14: ppcno3a9 sans élément à l'équilibre
! 2 - lim_zc.f, ligne 455, lim == 0 --> lim <= 0
! 3 - corrections de bugs dans ppcno12, ppcno3a14, opa
! 4 - opa_yveline: utilisation de inside pour éviter les sorties de table
! 5 - Création du nom de dessin des_mS, permettant d'afficher les abondances
!     de surface et non les valeurs maximales dans la cadre de des_m
! 6 - Pour les Cépheides, hors de la zone d'instabilité, limitation du pas temporel
!     à 0.03My  si 1.2d8 < Tc < 1.5d8. Introduction de la matrice P dans least_sq.f
!     Dessin de la zone d'instabilité avec des_m
! 7 - Introduction de la rotation de JpZh: omega_jpz, corrections de divers dégats
! 8 - Introductions de log10_teff, log10_teff_t, d_teff pour l'introduction d'un
!     contrôle du pas temporel en fonction des variations de Teff, cas de Céphéides.
!     Introduction d'un contrôle du pas temporel sur les variations de L.
!     Incidences sur le fifichier réglages
! 9 - Création de la routine cepheid_instb.f: zônes d'instabilité, modifications
!     dans cesam, evol, des_m, pulsation, des2k_hr.

!27/12/09 CESAM2k_V3.3.4 (commit S5)
! 1 - Création du fichier de sortie ASCII mon_modele.mHHe limites de H et He
!     shells burning. Pour qu'il soit créé, dans le *.don coder
!     nom_des='des_mHHe' ou 'des_mCHHe' ou 'des_rHHe'
! 2 - Création de la routine ppcno3a10: ppcno3a9+Li7
! 3 - Création de la routine pulsation.f du module mod_evol:
!     détermine si un modèle est dans la zone d'instabilité des Céphéïdes.
!     Utilisation dans ceam.f pour création du fichier mon_modele.cephe, list.f
!     pour écritures, evol.f pour contrôle du pas temporel. Création du type de
!     précision 'ce' pour Céphéïdes. En codant PRECISION='ceO' on forme le
!     fichier ASCII de chaque modèle dans la zone d'instabilité.

!05/11/09 CESAM2k_V3.3.2 (commit S5)
! 1 - mod_numerique: implantation de la routine least_sq.f résolution d'un
!     système linéaire surdéterminé par moindres carrés
! 2 - introduction de la compilation de réac. nuc. NACRU (NACRE updatée)
!     utilisant des taux updatés et, par défaut, les taux origines de NACRE
!     incidences dans taux_nuc et tabul_nuc, introdution du taux modifié de
!     N14(p,g)O15(e+ nu)N15 cf. Yveline Lebreton
! 3 - Introduction dans f_rad.f de frad_test permettant des tests de Burgers
!     avec accélérations radiatives
! 4 - Dans opa.f, utilisation de opa_compton dès que T > 4.d7, totalement ionisé
! 5 - Systématisation de l'écriture identifiant l'appel quand il y a sortie de
!     table dans bsp1dn
! 6 - difm_br.f : rectification d'une erreur concernant le signe du gradient de
!     pression de radiation qui est opposé à l'accélération radiative
! 7 - Introduction des accélérations radiatives dans diffm_mp
! 8 - evol.f : ajout contrôle du pas temporel portant sur la variation du flux

!05/11/09 CESAM2k_V3.3.1
! 1 - etat_irwin.f: l'extension "1a" est remplacée par "6"
! 2 - list.f, cesam.f: le modèle est-il dans la zône d'instabilité? affichage
!     de la période
! 3 - opa_yveline_lisse.f: en cas de sortie de table on se place aux bords

!17/07/09 CESAM2k_V3.3.0
! 1 - Le fichier de données est modifié:
!   a) Suppression de lim_jpz du fichier de données *.don
!   b) Introduction de la possibilité d'un arrêt sur Rtot quand Rtot est traversé
!    dans un sens quelconque. Aménagements dans
!    lit_nl.f, mod_donnees, mod_variables, cesam.f, resout.f, update
!   c) La précision 'aj' (ajuste) est supprimée, l'ajustement des arrêts sur Tc, Xc,
!    Rtot.. est effective pour tous les types de précision en incluant un "A",
!    dans la désignation de la précision dans le fichier de données *.don
!    Exemples: PRECISION='saA', 'rgA', 'prA' etc... sans ajustement continuer à
!    coder 'sa', 'rg', 'pr' etc...
!    modifications dans lit_nl.f, cesam.f, des_m.f, fcmax.f, lim_zc.f,
!    lit_binaire.f, mod_donnees.f
! 2 - mod_atm, coll_atm : nom_atm-->var_atm (confusion avec l'atm de restitution)
! 3 - coll_atm : amélioration de l'algorithme des corrections
! 4 - eq_atm : aménagement des tests de dérivation
! 5 - cesam.f: diminution du psi0, ordre=4 pour certaines précisions
! 6 -des_m/r : aménagements d'échelles pour le diagramme HR
! 7 -Simplification de l'usage de la chaîne de caractères nom_qs, aménagements
!     dans coll_qs.f, mod_donnees.f, cesam.f
! 8 - Aménagements pour le calcul de Ex. dans toutes les routines de réactions
!      thermonucléaires pp*.f
! 9 - cesam.f introduction du fichier *.cone et commentaires de combustion du néon
! 10 -mu_mol.f : ajout en sortie de mue, aménagements dans thermo, dgrad, difft_smc
! 11 -thermo.f: ajout de mu, mue, dmu_x en sortie, aménagements dans cesam,
!     ecrit_ascii, lim_zc, static_m, static_r
! 12 -taux_nuc.f : tabulation NACRE+ avec corr. de Ogata pour 3alpha réaction 15
!      He4(2a,g)C12
! 13 -Introduction de precic pour arrêt des itérations NR pour la composition
!      chimique / gramme, incidences dans rk_imps.f, mod_donnees.f, resout_chim.f,
!      cesam.f, fichier réglages
! 14 -evol.f : contrôle des variations de composition chimique basé sur la variation
!      max de l'abondance / gramme; incidences sur evol.f, rk_imps.f, cesam.f
! 15 -mvt_dis  -- > no_dis incidences dans cesam.f, mod_donnees.f, evol.f,
!      fichier "réglages"
! 16 -Suppression de ab_min et de fmin_abon, incidences dans mod_donnes.f, cesam.f,
!      evol.f, opa_opal_CO.f, rk_imps.f, resout_chim.f, pp*.f, fichier réglages
! 17 -Création du nom de dessin des_mC, permettant d'afficher les abondances
!      centrales et non les valeurs maximales dans la cadre de des_m; incidences sur
!      des.f et des_m.f
! 18 -Introduction de l_demi dans le fichier de réglages
! 19 -Introduction des fichiers arguments_du_point_don.f et arguments_des_reglages.f
!      dans le directory SOURCE pour faciliter leur construction
! 20 -mod_opa.f et routines d'opacité: introduction du LOGICAL dehors signalant les
!      sorties de table, diminution du pas temporel dans certaines situations
! 21 -Identification de l'origine de l'appel en cas de sortie de table d'opacité
! 22 -Création des routines resout_lisse_dif.f et eq_difg.f pour le lissage par
!      diffusion/galerkin de la composition chimique. Utilisation suivant le type
!      de précision ou en ajoutant la terminaison L au type de précision
!      Exemples: prL, saL, etc..
! 23 -Réorganition du fichier de réglages
! 24 -Exemples de fichier de données "*.don" et de réglages "*.rg" dans le
!      sous-directory SOURCE
! 26 -tabul_nuc.f: mise à jour de la description des réseaux
! 27 -abon_ini.f: introduction des tables d'abondances de Asplund Grevesse, Sauval
!      & Scott 09 arXiv:0909.0948v1 du 04.09.09, mise à jour de rapports isotopiques
! 28 -mod_numerique.f implantation des fonctions génériques entre ('entre_d.f,
!      entre_s.f, entre_i.f) et inside(inside_d.f, inside_s.f)
! 29 -evol.f et lim_zc.f amélioration de la gestion des limites ZR/ZC
! 30 -cesam.f: introduction de age_car, âge caractéristique à partir duquel un
!     modèle est supposé relaxé des conditions initiales. Aménagement des types
!     de précision
! 31 -Introduction de m_zc23 dans mod_variables, calcul direct de m_zc23 dans lim_zc
!     simplification de l'exploitation de m_zc23 dans evol
! 32 - Introduction des routines d'opacité conductive opa_cond_iben <-- kappa_cond
!      et opa_cond_mestel+cond_mestel, modifications dans mod_opa.f, lit_nl.f,
!      write_nl.f, mod_exploit.f, mod_donnees.f, opa.f
!      Ex: NOM_OPA_COND='opa_cond_mestel
! 33 - Implantation de l'EOS etat irwin, explications dans etat_irwin_explik.f

!07/07/09 CESAM2k.V3.2.13
! list.f : correction d'un bug concernant Yc, ajout de Xc, écriture de Z/X_sol
! Dans pp1.f on impose mu_saha=.FALSE.
! Dans diffw*.f: corrections de commentaires
! evol, resout:  modif. limitation pas temporel en fonction de delta ln Tc

!15/11/08 CESAM2k.V3.2.12
! 1 - bsp_gal.f, newspl_gal.f : optimisation, introduction de mg
! 2 - eq_ini_rota4.f, initialise_rota4 : initialisation de U par eq. diff. lin.
! 3 - eq_diff_rota4.f : condition limite externe sur Tau au lieu de U=0
! 4 - coeff_rota4.f, eq_diff_rota4.f : multiplication par nu des coefficients de
!     l'équation de U
! 5 - remplacement des routines conv_cm_reza.f et conv_cgm_reza.f, (corrections
!     communiquées par Reza Samadi)
! 6 - fcmax.f: le cas precision = 'mx' n'était pas envisagé
! 7 - eq_diff_chim.f: mcz_ext --> nuzc_ext
! 8 - evol.f: déplacement des définitions de nuzc_ext et mzc_ext,
!     suppression des "bricolages techniques" pour la base de la rotation
! 9 - mod_atm: extension des tables de Marcs introduction/modifications de
!     trho_4000.f, marcs.f, mod_atm.f (B.Pichon), correction d'une coquille
! 10 - pertm*.f: reformulation des pertes/gains de masse, pertm_waldrom.f
!      correction d'un (petit) bug dans la détermination de la luminosité.
! 11 - newspl_gal.f: appels à bsp1dn limités à l'intervalle [x(1),x(nx)],
!      augmentation du nombre d'abscisses pour intégration de gauss
! 12 - ppcno10.f : Li7 est hors équilibre
! 13 - diffus.f: suppression des "bricolages"
! 14 - eq_diff_rota4.f : optimisation de la condition limite de Phi en surface
!      condition limite dU=0 au centre
! 15 - cesam.f : aménagement des commentaires de la logique de baratine,
!      grille_fixe est imposé avec diffusion du moment angulaire
!      Il est possible de poursuivre un calcul avec des résaux nucléaires
!      différents si les éléments chimiques pris en compte sont identiques Ex
!      Exemple: ppcno9 et ppcno3a9. Amélioration du calcul des pourcentages
!      d'énergie PP, CNO..GRAV, incidences pour les écritures et dans list.f
! 16 - lim_zc.f : adaptation du changement de grille fixe
! 17 - abon_ini.f : introduction des mixtures photosphérique et météoritique
!      de Asplund Grevesse Sauval 05, table crée par Yveline Lebreton
! 18 - calib2k_zams.f, calib2k_pms.f introduction de la possibilité de modifier
!      le Z/X de calibration
! 19 - implantation de diffw_toul.f : coefficients de diffusion simplifiés pour
!      la diffusion du moment cinétique
! 20 - marcs.f, roger.f : écriture de la valeur utilisée de [Fe/H]
! 21 - Aménagements dans la notice
! 22 - module mod_exploit : des2k_abon.f introduction du dessin de grad_mu,
!      création de des2k_dvson.f
! 23 - Suppression des routines ctes*.f, les initialisations sont effectuées
!      directement dans ini_ctes. Les noms des ensembles de constantes
!      i.e. ctes_94 restent inchangés. Création d'ensembles de constantes
!      ctes_31 (de Toulouse), ctes_ba (Basu & Antia) et ctes_gaia.
! 24 - lit_nl.f : possibilité de faire le calcul avec rot_0 "ie.sans rotation"
!      bien que la vitesse angulaire initiale w_rot soit (par erreur) non nulle
! 25 - Introduction de la correction de Jorgen (astro-ph 0811.100v1) coder
!      NOM_OPA=opa_yveline_jorgen
! 26 - write_nl.f : utilisation de grad_ovi/s pour le signe de ovshti/s
! 27 - Atmosphères roger & marcs : utilisation des valeurs limites en cas de
!      sorties de table en teff et/ou logg, adaptations dans tdetau.f,
!      trho.f, trho_4000.f
! 28 - Introduction de tacho_ext/full.f: diffusion sous/sur et sous les ZC
!      selon les Toulousains, pour Li7 solaire calibration, utilisation avec
!      difft_nut_ext/full

!20/09/08 CESAM2k.V3.2.11
! 1 - déplacement de coeff_vth de mod_etat --> mod_evol
! 2 - diffw_mps.f prise en compte du statut OPTIONAL de dhv
! 3 - coeff_rota4.f correction sur chi_t
! 4 - coeff_rota.f, coeff_rota4.f, eq_diff_rota4.f suppression des dérivées
!     dhv et dfrl inutilisables
! 5 - lit_nl.f : suppression de l'utilisation de la condition limite JPZ
!     pour la rotation (il conviendrait de passer à la version V3.3...)
! 6 - mise à jour des SCRIPTS_BASH

!09/08/08 CESAM2k.V3.2.10
! 1 - ecrit_ascii.f zbar --> z_bar ligne 366
! 2 - bsp_dis.f n --> nx ligne 81
! 3 - cesam.f psi0 : 0.07 -->0.08 par défaut,
!    Krot=3,4 : l_demi=.FALSE., ini0=3
! 4 - coeff_rota.f : dfrl(Krot,.. --> dfrl(3,..
! 5 - diffw.f dhv en OPTIONAL, coeff_rota3 suppression de dhv
! 6 - eq_diff_rota4: equation de psi dans les ZC
! 7 - evol : simplification de la gestion des ZC/ZR

!14/07/08 CESAM2k.V3.2.9
! 1 - cesam.f : aménagement des réglages pour la précision 'av'
! proportions d'énergies en % de l'énergie nucléaire totale
! ajustement des réglages avec Krot=3,4
! 2 - evol.f : amélioration de l'approche de la fin de la MS,
!  introduction de l_demi dans les précisions (cesam.f, mod_donnees.f)
!  amélioration de l'algorithme imposant le mélange radiatif de l'enveloppe
! 3 - lim_zc.f : aménagement de la suppression des limites ZR/ZC trop proches
! 4 - diffw.f : suppression de deff, dv, dh >= 50

!06/05/08 CESAM2k.V3.2.8
!1 - list : attribution de SAVE à fes, bug signalé par B.Pichon
!2 - cesam : introduction de la possibilité de poursuivre avec diffusion
! un modèle initialement sans diffusion et réciproquement (souhait de
! F.Thévenin). Limitation du pas temporel maximal suivant la valeur de la masse
! initiale. Introduction de l_demi. Création du fichier ASCII 4d-2.pms pour
! initialisation PMS de modèles M>10Msol, et ZAMS m=15Msol, le fichier ASCII
! 4d-2.pms et m150.zams sont mis dans le sous directory EXPLOIT
!3 - routines de réactions thermonucléaires : suppression de rot_solid inutile
!4 - eq_diff_chim : réintroduction de Deff avec rotation, voir commentaires
!5 - lim_atm : adjonction d'une loi P(tau) pour Teff > 12000
!6 - mod_donnees : augmentation de pnzc=10. Introduction de l_demi
!7 - evol : introduction de demi, simplification de la logique de suppression
! d'une ZR externe, déplacement de la logique des mvt_dis
!8 - opa.f : la limite pour les opacités compton est ramenée à de 7 à 6.D7
!9 - upddate : suppression de la limitation du pas temporel qui existe dans evol

!01/05/08 CESAM2k.V3.2.7
!1 - etat_opalX/Z : correction de dut=duro*drot+eos(4)*1.e6 au lieu de eos(8)
!    erreur signalée par Yveline, Joséphina, Laurent
!2 - ecrit_ascii correction d'un bug signalé par Yveline. Plantage lors de la
! formation d'un fichier ASCII pour un modèle totalemnt convectif avec ajout
! d'un point au voisinage du centre
!3 - coeff_vth :  adaptation au cas en_masse=.FALSE.

!25/04/08 CESAM2k.V3.2.6
!1 - cesam.f : prise en compte de blabla
!              aménagement des paramètres pour modèles en eulérien
!2 - static_r : remise en fonction
!3 - lim_zc : correction d'un bug affectant la suppression des noyaux radiatifs
!             ou convectifs d'étendue insignifiante

!20/03/08 CESAM2k.V3.2.5
!1 - ascii.f : correction d'un bug signalé par B.Pichon :
! ligne 144 lire evar(2,:)=EXP(evar(2,:)) et non evar(2,:)=10.d0**evar(2,:)
!2 - le fichier journal change de nom, il devient journal.f

!10/02/08 CESAM2k.V3.2.4
!1 - Rotation : rassemblement des routines resout_rota3/4 dans resout_rota.f.
! Création de la routine jacob_rota.f; conditions limites identiques pour les
! formalismes tz97 et mz04; mise en place de l'option lim_jpz pour les deux
! formalismes. Introduction de ord_rot=m_rot+r_qs pour diffusion du moment
! cinétique, ord_rot=m_rot sinon; adaptation de lit_binaire.f
!2 - Suppression de la possibilité d'utiliser les anciens fichiers de données.
!3 - Mise à jour du MAKEFILE du sous-directory EXPLOIT
!4 - Correction d'un commentaire dans marcs.f (cf. B.Pichon).

!01/01/08 CESAM2k.V3.2.3
!1 - cesam.f : Création de la précision 'aj' (ajuste) précision réaliste ('pr')
! avec ajustement du dernier pas temporel pour obtenir  Tc, He core, Xc à
! la fin d'une évolution.
!2 - mod_numerique : création de la routine bsp_gal, permettant la
! détermination des coefficients splines d'interpolation de n fonctions connues
! en tout point. Le but étant d'obtenir une estimation des dérivées.
! Création de la routine de changement de base newspl_gal.
! Ces routines utilisent le formalisme des éléments finis Galerkin.
!3 - mod_evol : création de la routine tab_vth, tabulation de grandeurs
! thermodynamiques sur le base de la composition chimique. Le but étant le
! calcul des dérivées spatiales de ro, mu.. d'où vaissala. Création de la
! routine coeff_vth permettant le calcul des quantités à tabuler.
! Les tableaux vth et vth_t sont des éléménts de mod_variables.
! Allocation, exploitation de ces tableaux à divers endroits : cesam,
! ecrit_ascii, evol, coeff_rota4.
! 4 - Aménagements pour l'atmosphère. Introduction des lois t(tau) hsra et marcs
! ces dernière sont disponibles sur demande auprès de B.Pichon). Correction
! de bugs dans trho et roger, améliorations de fesh+00.data et fesh+02.data
!(B.Pichon)
! lim_atm : rétablissement de la limite tau_max=20 avec loi T(tau) non
! purement radiative

!01/10/07 CESAM2k.V3.2.2
!1 - Implantation de la routine mu_mol (Module mod_etat) de calcul de diverses
! quantités liées à µ. Introduction dans la routine dgrad pour le critère de
! Ledoux, et implantation dans thermo, lim_zc, for037. Introduction de mu_mol
! dans le MAKEFILE du directory EXPLOIT.
!2 - Amélioration de la définition de phi, incidences dans coeff_rota3/4,
! eq_diff_chim, cesam. et ecrit_ascii sur la fréquence de BV avec new_bv
! Utilisation de new_bv=TRUE pour toutes les précisions sauf 'sa'.
! Dans les fichiers de sortie ASCII, s'il n'est pas tenu compte de la pression
! turbulente, replacement de Ptot/Pgaz=1 [var(23,i)], par grad_mu=dln mu/dln P.
!3 - Le critère de Ledoux peut être utilisé avec et sans convection.
!4 - Adjonction de g, msol, rsol et lsol dans les glob des fichiers ASCII de
! sortie.
!5 - Dans lim_zc : réduction du nombre d'allocations, déplacement de la
! déallocation de new toujours alloué (ligne 1157,remarque de B.Pichon).
! Déplacement du calcul et des allocations de v_son et r_son
! (remarque de B.Pichon).
!6 - Implémentation de la semi-convection : routine difft_smc du type difft.
! Introduction dans le MAKEFILE du directory EXPLOIT.
!7 - Dans dnunl suppression des allocations (remarque de B.Pichon).
!8 - Dans cesam, update, (reglages) introduction de la variable de "précision"
! d_tc limitant la variation de LOG10 température centrale sur un pas temporel,
! Correction d'une erreur pour le calcul BV dans la cas new_bv=.TRUE..
!9 - Dans lit_nl : suppression de la possibilité de l'utilisation des anciens
! fichiers de données (suggestion de B. Pichon).
!10 - Sous-directory EXPLOIT création du programme de dessin des2k_grad,
! des2k_diff_spl, amélioration de f037_2k, et de divers programmes de dessin.
!11 - Mise à jour de la NOTICE.
!12 - Sous-directory SCRIPTS, création de scripts en BASH.

!14/07/07 CESAM2k.V3.2.1
!Dans cesam :
!1- Avec diffusion du moment cinétique (Krot=3, 4) afin de laisser deff jouer
! pleinement son rôle, re_nu est fixé à 0
!2- Introduction du facteur fmin_abon, défini par la précison ou les réglages,
! fmin_abon permet de régler ab_min dans le but d'une amélioration
! du pas temporel lorsque H brûle en couche. Précédemment
! fmin_abon=0.01 était fixé.
!3- Dans ecrit_ascii, ajout de points dans les couches centrales,
! localisation suivant les paramètres q0 et l0
!4- Dans mod_evol, evol : introduction de coll_rot et colpnt_rot pour la
! détermination des points de collocation pour la rotation. Un décentrage des
! points au milieu des l'intervalles permet de stabiliser la solution avec
! m_rot=2, 4,..pair
!5- Dans EXPLOIT : introduction de coll_rot et colpnt dans le MAKEFILE
!6- Dans evol : amélioration de l'algorithme de l'adaptation des la ZC externe.
!7- Dans coeff_rota3/4 phi_e ne peut être nul, + schéma implicite

!01/07/07 CESAM2k.V3.2.0
!Dans le fichier reglages : ajout de l0 points à disposer autour des
! discontinuités dans les fichiers de sortie ASCII, ajout de la variable
! new_bv calcul de la fréquence BV utilisant Phi (d ln ro / d ln mu),
! avec, faute de mieux, Phi=1
!Dans ecrit_ascii : adaptation pour ajout des points
!Dans bsp_dis et noeud_dis, adaptation au formalisme de cesam2k
! du traitement des discontinuités dans les routines
!Dans evol : s'il y a retrait du coeur convectif
! 1- sans diffusion, interpolation linéaire de la comp.chim. dans la
! la zone de retrait,
! 2- avec et sans diffusion suppression de la discontinuité au niveau du
! raccord ZC/ZR.
! Mise en réserve du traitement de l'augmentation de l'abscisse inférieure
! des ZC externes.
!Dans cesam : m_ch=2 (interpolation linéaire) de la composition chimique
! pour toutes les précisions
!Dans mod_numerique implantation de max_local : détermination des maxima locaux
! pour normalisation dans les dessins
!Dans EXPLOIT, modification des programmes des2k_vaiss et des2k_abon, mise en
! place de max_local dans le MAKEFILE

!14/06/07 CESAM2k.V3.1.3
!Ajout de la routine ctes_94m, identique à ctes_94 avec des valeurs entières
!pour les masses atomiques.
!Dans la routine ecrit_ascii, ajout éventuel d'un point à une distance q0 du
!centre. Ajout de q0 dans le fichier reglage.
!Dans le sous-directory EXPLOIT :
! Installation du zoom dans le programme des2k_abon.
! Modification du zoom dans le programme des2k_vaiss.

!25/05/07 CESAM2k.V3.0.2
!Dans eq_diff_rota3/4, resout_rota3/4 :
! aménagement des conditions limites
!Dans sous-directory EXPLOIT :
! update des jacobiens pour les calibrations en tenant compte de la diffusion du
! moment cinétique, programmes calib2k_pms, calib2k_zams
!Dans le sous-directory NOTICE :
! correction de quelques points, des conditions limites pour la rotation

!03/05/07 CESAM2k.V3.0.1
!Dans cesam.f :
! correction d'un bug concernant les proportions pp, CNO, 3a etc..
! m_rot=2 pour la précision 'np'
!Dans coeff_rota4 : élimination du cas nui <= 0
!Dans coeff_rota3 : introduction de phi/delta*gad_mu dans dgrad (remarque de Andy)
!Dans eq_diff_rota3 et eq_diff_rota4 :
! correction d'un bug concernant la dérivée de bs(2)
! calcul direct de Lambda, reprenant (idée de Phi)
!Mise en place de la variable logique baratine permettant de
! détourner une grande partie des commentaires 'on line' vers les fichiers
! mon_modele_static, *_atmos, *_evol pour, respectivement, l'équilibre
! quasi-statique, l'atmosphère, l'évolution de la composition chimique et
! la rotation (idée de B.Pichon).
!Dans lit_nl introduction de la possibilité de lire un fichier .don
! ancien et un fichier .don sans diffusion ni rotation (idée de B.Pichon).
!Création du fichier pgplot_factice.f permettant d'éviter l'utilisation du
! PGPLOT pour les dessins "on line"

!03/05/07 CESAM2k.V3.0.0
!Pour la rotation, remplacement de la méthode des éléments finis par celle
!de collocation
!Pour des initialisations, les fichiers binaires *.hom, *.pms sont utilisables

!03/02/07 CESAM2k.V2.5.2
!Pour la rotation : mise en place de formalismes identiques dans ZR et ZC
!avec, dans ZC, des coefficients de diffusion >> 1

!04/02/07 CESAM2k.V2.5.1
!Modification du fichier reglages
!Dans evol.f :
! Suppression de demi
!Dans cesam.f :
! Avec Krot=3,4 on ignore les discontinuités de composition
! chimique et on les lisse mvt_dis=.FALSE., lisse = .TRUE.
! pour la diffusion du moment cinétique utilisation de la formulation approchée
! de mu(Krot=3,4), mu_saha=.FALSE. (sauf avec precision='rg')
! Pour précision réaliste & corot (pr & co),
! fonction d'espacement limitée à Ln P, ajustement de psi0 pour avoir un nb de
! couches de l'ordre de 1000 pour une 1.5Msol sur la ZAMS
! addition de ajuste et lisse dans le fichier reglages
!Dans base_chim.f :
! Utilisation d'une base continue non dérivable
!Dans base_rota.f :
! Utilisation d'une base continue non dérivable
!Dans print_ctes.f :
! permutation de Li6 et Li7 dans la liste d'écriture
!Dans coeff_rota3.f, coeff_rota4.f, eq_diff_poisson.f :
! allocation du tableau ion, (cf. B.Pichon)

!19/01/07 CESAM2k.V2.4.6
!Dans coeff_rota4
! Addition de pi dans cte5_0
! Changement de signe de cte_0(4) Bug signalé par Yveline
! Corrections des dln_mu Bug signalé par ANDY
! Addition de mu pour les tracés :
!  Augmentation à 30 de nb_coeff dans mod_evol
!  Modification de des2k_rot dans EXPLOIT
!Dans resout_rota4, ligne 345 2 au lieu de iv

!11/12/06 CESAM2k.V2.4.5
!Dans coeff_rota3 et 4
!Introduction d'une approximation de dln ro / dln mu
!Utilisation de l'approximation numérique de d ln P/d nu
!Dans cesam.f, énergies pp, cno, 3a+C+O, grav en Lsol
!Dans evol permutation de l'ordre des intégrations de diffusion
!composition chimique puis moment cinétique
!Dans coeff_rota4, utilisation de chi_T
!Dans resout_rota4 suppression du controle des corrections NR

!06/12/06 CESAM2k.V2.4.4
!Aménagement des tests de dérivation dans coeff_rota4, eqdiff_rota4
!Correction d'Yveline : dans cesam.f, utilisation de -TdS au lieu de ABS(TdS)
!pour l'estimation de l'énergie graviphique.

!25/11/06 CESAM2k.V2.4.3
!Amélioration de la logique du dtmax dans cesam.f
!Suppression des discontinuités dans la base de la rotation
!les variables deviennent continues non dérivable ==> suppression du
!rétablissement des continuités dans resout_rota3/4.
!Augmentation du nombre de variables dans le fichier ASCII *_coeff-rota.dat
!Aménagement du programme des2k_rot du sous directory EXPLOIT pour le tracé des
!des variables de la rotation
!Limitation à 50My du pas temporel s'il y a diffusion du moment cinétique
!et limitation à 10% la variation du pas temporel
!Augmentation/diminution du nombre de couches limité à 5% (au lieu de 10%)

!27/10/06 CESAM2k.V2.4.2
!rectification lnt923" --> lnt932  réaction 5 de NACRE dans taux_nuc.f
!rectification des commentaires concernant mu_e, var(15) est mu_e
!Mise à jour de la notice et de l'aide mémoire

!12/09/06 CESAM2k.V2.4.1
!Introduction de la variable logique nc_max permettant d'imposer le nombre
!maximum de couches n_max pour le calcul du dernier modèle : nc_max=n_max < 0
!Dans cesam.f, introduction de la variable coox et du fichier *.coox
!combustion de l'oxygène
!Dans opa_houdek9 appel à opa_opal2_co en cas de sortie de table (pour
!l'atmosphère opa_opal2_co, utilise opa_yveline)

!22/08/06 CESAM2k.V2.4.0
!Introduction des réactions nucléaires de la combustion du carbone
!C12(C12,..).. C12(O16,..) et de l'oxygène O16(O16,..).
!Intervalle de tabulation des taux des réactions fixé à 1MK
!Création de la précision 'av' avec des aménagement permettant d'atteindre
!les stades avancés
!Création des routines ppcno3aco  ppcno3acos de la combustion de H à O
!Introduction du vecteur iter_qs permettant d'adapter les variables contrôlées
!pour la résolution de l'équilibre quasi-statique
!Création de la routine opa_compton et utilisation dès que Z > 0.098
!Ecriture de Teff dans le dessin
!Calcul des poids statistiques dans saha
!Dans des_m abondances centrales avec 'av'

!05/07/06 CESAM2k.V2.3.3
!Introduction des réactions nucléaires de NACRE pour le 3alpha et C12(a,g)O16
!Début de la combustion du carbone mise à 6d8K, tabulation pour ppcno3ac10
!repousée à T9
!Création de la précision 'av' : long runs pour les longues évolutions
!Réaménagement du module mod_nuc
!Mise en place de SAVE dans etat_opal* (remarque de JP Marques)

!28/06/06 CESAM2k.V2.3.2
!Modification de la localisation des limites ZR/ZC dans lim_zc
!Pour le calcul des coefficients de diffusion du moment cinétique, utilisation
!de Omega et U au temps t

!13/06/06 CESAM2k.V2.3.1
!Adaptation de l'expression approchée de \phi
!Possibilités de dessin on line et off line des coefficients de la diffusion
!du moment cinétique

!24/05/06 CESAM2k.V2.3.0
!Introduction du formalisme de diffusion du moment cinétique
!selon Mathis & Zahn 2004 occasionnant diverses adaptations du fichier de
!données et de réglages.
!Création d'une routine de dessin des variables de la diffusion du moment
!cinétique
!Utilisation de l'ancienne formule de la fréquence de Brunt Vaissala
!Correction d'un bug dans les routines PPCNO12, PPCNO12Be, PPCNO12BeBFe,
!PPCNO12Li
!Accélération de la convergence de la diffusion du moment cinétique

!24/04/06 CESAM2k.V2.2.0
!Restructuration permettant différents calculs de la vitesse angulaire
!Modification du fichier de données
!Introduction du formalisme de Matis & Zahn 2004 (début)
!Introduction de la conservation locale du moment cinétique

!11/04/06 CESAM2k.V2.1.0
!Vérification du Jacobien de ppcno3a12Ne (bug)
!Redéfinition de scale dans rkimps
!Suppression, dans cesam, de Kipp=.TRUE. pour les modèles après la séquence
!principale
!Suppression, dans cesam, de n_max de la NAMELIST nl_rlg

!20/03/06 CESAM2k.V2.0.8
!Aménagements dans etat_opal, ZFSinterppeos, opal_ascii_bin, calib2k
!Suppression des SUM dans opa_opal2 (bug)
!Compléments de formules de réactions nucléaires NACRE
!Aménagements mineurs dans des_m, des_r, resout, z14xcotrin21

!20/03/06 CESAM2k.V2.0.7
!Création du programme d'exploitation des2k_opa.f
!Création de la routine générique coeff_rota.f et des routines coeff_rota_saha.f,
!coeff_rota_ioni.f, coeff_rota_z16.f
!Dans ces routines mise à 0 de chi_T mal calculé avec les données dont on dispose
!Création de la routine difft_sun.f
!Addition de la viscosité cinématique au coefficient Dv
!Création de la routine coeff_rota_ioni.f
!Introduction de mini, valeur minimale de Dv, dans le module mod_evol
!Dans lit_nl, avec diffusion du moment cinétique, on impose D_turb >= 50
!Introduction du nom de la routine de calcul des coefficients de rotation dans
!le type de précision et dans le fichier réglages
!Introduction dans mod_donnees de la variable logique ecrit_rot conditionnant
!l'écriture du fichier mon_modele_coeff_rota.dat pour dessin des coefficients
!de rotation, addition de ecrit_rot dans le fichier reglages
!Création de la routine difft_gab.f
!Mise à jour des paramètres de précision et de l'aide mémoire
!Aménagements du test de sécurité MODIF_CHIM dans abon_ini.f vent.f,
!planetoides.f

!02/02/06 CESAM2k.V2.0.6
!Suppression de w_form du fichier de données, mis dans le fichier reglage
!Création de la routine générique coeff_rota.f appelant l'une des deux routines
!coeff_rota_saha.f et coeff_rota_z16.f
!Mise en place des chutes de planétoïdes : modification de la composition
!chimique de la ZC externe et apport/retrait de moment cinétique, adaptation du
!fichier planet

!01/02/06 CESAM2k.V2.0.5
!Limitation de la source du vent à la ZC externe, suppression du
!paramètre p_vent du fichier vent, simplification du traitement du vent,
!aménagement du programme fichier_vent.f et des fichiers exemple.vent et vent
!Mise en place des chutes de planétoïdes, création de la routine planétoïdes et
!des fichiers exemple.planet et planet
!Ajonction de la masse terrestre dans les fichiers ctes85 et ctes94
!Ajonction de W_FORM dans le fichier "reglages"
!En cours :
!Création, suppression, aménagement et mise au point des divers routines et
!programmes concernés par la diffusion du moment cinétique, principalement :
!cesam, ecrit_rota, diffus,
!resout_rota, resout_chim, coeff_rota, eq_diff_rota, eq_diff_chim, des2k_rot

!15/12/05 CESAM2k.V2.0.4
!Introduction de l'argument optionnel duale de la routine newspl
!Construction et aménagements de diverses routines pour la diffusion du moment
!cinétique

!20/10/05 CESAM2k.V2.0.3
!Aménagements mineurs dans inter
!Correction du calcul de d ln l / d ln m, et d ln ro / d m^2/3
!Possibilité de calcul direct des coefficients de diffusion du moment cinétique
!Introduction de tab_coeff_rota
!Introduction dans resout de la variable et fonction logique cmax et fcmax pour
!l'utilisation du nombre maximum de couches avant de sortir
!Suppression du dessin de ro sans diffusion (discontinuités)
!Dans evol, sans diffusion,léger lissage par contour de la composition chimique
!pour lisser le retrait des ZC
!Ajout de lim et de model_num à la fin des fichiesr binaires

!13/10/05 CESAM2k.V2.0.2
!Calcul direct de Deff dans eq_diff_chim
!Permutation de l'ordre diffusion du moment cinétique <==> diffusion des
!éléments chimiques
!Création du programme de dessin des2k_coeff_rota
!Ecriture du numéro du modèle dans des_m et des_r
!Dans les fichiers de sortie ASCII, ajout de 20 points de grille de part et
!d'autre des limites ZR/ZC pour affiner le profil de la fréq.BV
!En abscence de diffusion amélioration de la formulation de la fréq.BV

!05/10/05 CESAM2k.V2.0.1
!Amélioration du choix de no_croiss dans noein
!Amélioration d'écritures et introduction de no_croiss dans linf
!SAVE pour les quantités ***0 des conditions limites de static_m

!01/10/05 CESAM2k.V2.0.0
!SAVE dans opa_yveline, etat_opalX, etat_opalZ
!Mise en place des numéros des modèles, sorties de tous les modèles en ASCII
!et .rep avec leur numéro, conservation du numéro dans les fichiers binaires
!Choix de grad_ad ou grad_rad dans les zones overshootées

!20/09/05 CESAM2k.V1.1.15
!Adjonction de v dans l'expression de teta dans les tests de dérivation de
!static_m et static_r
!Création de la routine ppcno3a13Ne22

!01/09:05 CESAM2k.V1.1.14
!Correction de dgravr dans thermo et thermo_atm
!Facteur 2/3 sur l'accélération centrifuge dans coll_atm et eq_atm
!Aménagements des équations relatives à la diffusion du moment cinétique
!Création des programmes de dessin des2k_dhve, des2k_rot, des2k_bin

!30/08/05 CESAM2k.V1.1.13
!Suppression de commentaires dans les modules mis dans la notice
!Suppression de nom_elem en dp des définition de mod_nuc
!Suppression de la variable pmw du module mod_donnees
!Suppression des tableaux xlim_rot, xcin et xcint du module mod_evol

!03/08/05 CESAM2k.V1.1.12
!Inversion de la chronologie du journal
!Corrections de bugs signalés par A.Moya dans coeff_rota :
!C12=1 et signes - pour C16,17,18
!Rétablissement de d2U/dnu2=0, et équation de diffusion de Omega dans ZC
!Déplacement de l'allocation de frot, dfrot dans lim_zc
!Annulation de C15*, C8 et C9 dans coeff_rota

!27/06/05 CESAM2k.V1.1.11
! coeff_rota, utilisation de rho ie. sans passer par l'équation d'état,
! pour cohérence avec dln ro
! des_m, augmentation du nb. de chiffres significatifs pour les abondances max
! diffm_mp changement de signe de l'accélération centrifuge
! diffm_mp & diffm_br coefficient de l'accélération centrifuge
! aménagements dans resout_rota, eqdiff_rota, diffus, coef_rota
! création du programme de dessin des2k_dhve du sous-directory EXPLOIT

!16/06/05 CESAM2k.V1.1.10
! Suppression du fichier *.atm pour initialiser ZAMS ou PMS
! Tracé de ro dans des_m et des_r
! Abondances des éléments au centre dans list
! Augmentation de m_rot --> 4 dans cesam.f pour tous les réglages
! Après la ZAMS on impose l'approximation de Kippenhahan

!13/06/05 CESAM2k.V1.1.9
! Correction C12(a,g)O16 dans ppcno3ac10
! Elargissement des dessins du HR jusqu'à nb nmax modèles
! etat_opal SAVE pour la variable iri
! Lxchim(nchim) dans etat_opal et etat_ceff
! lit_nl, nb_max_models pour lit_nl_2korg
! Réduction du pas temporel à la fin de la ZAMS
! Augmentation du nombre de couches après la TAMS, He et C burning
! Définition de z_table=z0 dans opa_opal2
! diffw routine PUBLIC de mod_evol
! Création de lit_binaire dans mod_exploit
! Création du programme des2k_dhve dans EXPLOIT

!03/06/05 CESAM2k.V1.1.8
! Rectification d'une virgule dans mod_exploit
! Suppression de la référence à compg(ihe4,1) dans list
! Save de cte1 dans colatm
! Mise du numéro de version dans journal

!31/05/05 CESAM2k.V1.1.7
! Mise en service du programme des2k_rot du sous-directory EXPLOIT
! SAVE et allocation des tables de données dans opa_yveline

!26/05/05 CESAM2k.V1.1.6
! Restriction de l'utilisation de ln ro = bp(7,:) au cas avec diffusion
! et ord_qs  2 à cause de la discontinuité de ro aux limites ZR/ZC
! Introduction du numéro de version par un include dans mod_donnees

!20/05/05 CESAM2k.V1.1.5
! Corrections dans lim_gong1, lim_tau1, add_ascii, opa_yveline lisse,
! pour réimplantation des modèles de GONG

!14/05/05 CESAM2k.V1.1.4
!Sous directory SOURCE:
! Correction du bug dans la formule de normalisation dans evol
! lire	chim(1:nchim,i)=chim(1:nchim,i)/norm
! et non  chim(1:nchim,:)=chim(1:nchim,:)/norm
!Sous directory EXPLOIT:
! Création des programmes de dessin des2k_abonts et des2k_abontc
! Suppression du programme des2k_abont
! Dessin de X, Y, Z dans des2k_abon

!05/05/05 CESAM2k.V1.1.3
! Correction d'un bug lié au calcul de Teff si n_atm=1 ie. dans les
! cas GONG1 et GONG2
! Calcul ("exact") de la féquence de Brunt-Vaissala en utilisant ln ro
! Suppression de contour

!02/05/05 CESAM2k.V1.1.2
! Amélioration des algorithmes gérant les arrêts sur t_stop et x_stop
! Implantation de l'arrêt sur He_core

!27/04/05 CESAM2k.V1.1.1
!  Correction de bugs engendrés par ihe4=-100 et lvent=.TRUE.avec PP1

!23/04/05 CESAM2k.V1.1.0
! Introduction dans resout de l'arrêt sur x_stop
! Introduction dans resout de l'arrêt sur t_stop