CHARACTER (len=7), PARAMETER, PUBLIC :: version='V2.4.4' !Signification du numéro de version : Va.b.c !a augmente si les fichiers binaires de reprise *.pms, *.rep, *.dat... changent !b augmente si l'un des fichiers de données ou de personnalisation change !c augmente si modification d'algorithmes (bug, nouvelle implantation, ...) !Journal des aménagements de CESAM2k !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 T > 0.007T9 !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