10/11/2008 (TG):
* Modification of the Makefile to include automatic test routine: make TEST or make test
* Added test_cepam.bash
* Added new subroutine Coeur/modelini.f (from resoutp.f)
* Cleanup of directories [svn: branches/convergence]

25/05/2008: 
* Corrected bug in derivative of lim_xa3t10 (was unstable & never used).
* Modified slightly test_lim_ext.f

17/05/2008: 
* Corrected a bug in etat_hhe: namelists were using yz_surf and yhe_atm instead of
  yzsurf and yhe. Strangely, this created an error when reading "etat_hhe.dat". 
* Created test_etat in tests (use Make test_etat)

16/05/2008:
* Makefile changed: added space after then at end of line (pbs with new term)
* rwbin changed: possibility to write/read ASCII files (.basc)
* src/cepasc2bin.f changed: transforms bin to .basc and vice versa
* com/ch2asc and com/ch2bin: 2 bash scripts to transform all files .basc <-> .bin
* a2b_EOS: strangely there was a logT_ instead of logT_1 => corrected
* Note g95 doesn't allow for variable formats (<n>f5.2)

12/05/2008:
* Makefile changed to include g95 compilation
* eos.f modified (contained a loop with a real variable)
* constantes.f modified (read with strings in format not allowed in g95)

03/04/2008:
* Pracmax added to LIM_NOYAU: limiting core pressure limit. Was 1d15, now 1d18.

18/11/2007:
* New subroutine LIM_FGRID.f for Uranus-Neptune atmospheric conditions
* STATICP_M.F modified small bug in the activate_core option

16/11/2007:
* ETAT_XYZ.F modified to check for either cepam_etat_xyz.don or cepam_etat.don
* OPA_F.F modified for a small bug in the handling of the rosseland opacity file

18/06/2007:
* OPA.F modified to read file "cepam_opa.don" with opacity multiplier

03/04/2007:
* Makefile modified in order to avoid saving .o files.

29/03/2007:
* RESOUTP.f modified so that there is no reinitialisation when dR/dt is positive
	may have to be improved on a case by case basis.

27/03/2007:
* CEPAM.f modified to include new routine call:
* EVOL_TEQ.f new subroutine to allow modifications of the equilibrium
 	temperature with time from optional file cepam_evol_teq.don

26/03/2007:
* STATICP_M.f modified to allow the possibility to have a core cooling
	as a function of temperature difference core-envelope.
	This is controlled by the file "cepam_lcentre.don"

26/03/2007:
* There was a convergence problem when you would change e.g. the 
	external boundary condition (e.g. from P=10 bar to P=1 bar). 
	Added a test that if the case, the initial model is
	reinterpolated (very roughly). 
=> RESOUTP.f modified accordingly. 


25/03/2007:
* Modified some small writing problems when mplanet<>0:
	mnoyau*mjup is the mass of the core in grams
	mplanet*mjup is the mass of the planet in grams
	Improved that in the following routines:
	*STATICP_M.f: small write statement
	*THERMOS.f: at the center: epsilon(1)=ljup/mjup/MPLANET
	*LIST_AUTRE.f: minor
	*COLNOY.f: we were saving mnoyau/mplanet; mnoyau now. 
	*LIM_NOY.f
	*EQNOY.f	
=> VERSION 3.0.4



23/05/2006:
* CEPAM.f: On ecrit le fichier d'evolution _ev.des a partir des
	premieres donnees uniquement (permet de savoir quand un modele
	n'a pas fonctionne).
* STATICP_M.f: Correction d'un petit bug lorsque activate_lcore=1 et mnoyau=0

12/12/2005:
* RESOUTP.f: Modification de la procedure de recherche de modele
	initial avec possibilite de demander un rayon initial fixe.

24/10/2005:
* STATICP_M.f: Ajout de la lecture de namelists dans le fichier "cepam_lcentre.don"
	Ces namelists ajoutent la possibilite d'integrer l'evolution
	planetaire en modifiant les conditions limites de luminosite a
	l'interface noyau/enveloppe

31/10/2004:
* FIGURES, FIGURES2...etc: Implementation d'un fichier d'option "cepam_figures.don"
	S'il existe, il permet d'enclencher le un systeme d'options de sortie
	Format: ligne 1> IOPTF / ligne 2> 1
	Options: 0: defaut/ 1: ecriture du fichier contrib<n>.dat

06/06/2004:
* CEPAM: Modification de la procedure de gestion des noms
	ATTENTION: les scripts automatiques doivent etre modifies en
	consequence !
* RWDES, RWBIN: Petites modifs (Attention: character*31 nom_fich)

25/08/2003:
* OPTMOD: Ajout d'un call sub_cep(3,...) -> creation des fichiers .don .bin
* RWDES: Ajout d'un common planetes
* CEPAM: Fichier .asc (unit=2) ouvert plus tot (sinon creation d'un fort.2)

22/08/2003:
* OPTIMISE: Ajout d'une option "calcul d'optimisation simple" (repondre "x")
* LISTP: N'ecrit plus le fichier .des -> RWDES
* RWDES: Changement des sorties du fichier .des (- de redondance).

11/08/2003:
* EQUATIONS: if (init) enleve
* Travail sur test_poly / repertoire results/Poly
* FIGURES

23-28/07/2003:
* CEPAM3.0.3
* Modifs sur magma/mambo (Santa-Cruz): 
* ETAT_HHE (commentaires pour cepam_etat.don mis a jour)
* CEPAM modele binaire: evolution.bin au lieu de evol.bin -si age_list(1)<0
        lecture des namelists precede rwbin, rwdes
* Ajout de la procedure rwdes: lecture d'un fichier ascii (.des)
* RESOUTP: pour tenir compte de rwdes
* STATICP_M: suppression des pauses
* RWDES: nouvelle subroutine
        

02-03/07/2003:
* CEPAM3.0.2
* Ajout de ierr dans gausdn, sbsp1dn, resoutp, cepam...etc.
* Optimise est bcp plus robuste: rappel de cepam lorsque ierr<>0.
* Optimise inclue les cas de Jupiter et Saturne
* lancement d'optimisations multiples sur suppe: voir repertoire
Cepam/don pour tous les fichiers a transferer. Pour executer, faire
csh don/suppe_multi.csh

10-18/06/2003:
* BUG: Dans compchim, on imposait en surface Y(Helium)=YZ(2) au lieu
de Y(Helium)=(1-Z)YZ(2). D'ou un decalage par rapport aux modeles
precedents.
* Modification de figures: tests de regularite du modele
* Figures3: on repart en general d'une solution iteree. On teste pour
les NaN avec "if (s2(i).eq.s2(i))"
(Faire attention au fait que les variables mm0 et ss0 contiennent
les valeurs de la masse, omega...etc. Il faut les reinitialiser a
chaque appel dans optimise!). 
* Attention: Figures4 n'a pas ete modifiee... A FAIRE

1/06/2003:
* BUG: conservation de omega dans "equations": modification de la
convergence des equations lors de la minimisation (rplanet<>rjup)
* Interpolation du potentiel centrifuge en 1-m: bcp + stable (1-m^2/3
serait encore mieux...)

30/05/2003: 
* CEPAM3.0.1
*Modif de l'interpolation de la composition chimique, du TdS en mu
*Interpolation du potentiel gravitationnel en 1-r au lieu de 10-r
(cosmetique!)
*Test d'une solution polytropique: precision limitee par la C.L. (Qqs
pbs de convergence a resoudre peut-etre?), mais < 8e-5 (avec 100,500 pts).
*Ajout de la routine etat_poly: equation d'etat polytropique. Quelques
modifs dans thermos pour permettre le calcul lorsque delta=0.

Mai 2003:
* optimise.f a ete modifie: creation d'un repertoire Optimise. 
Convergence par calcul de Jacobien (Newton-Raphson), puis simplex. 
* colnoy.f modifie pour plus de robustesse, on se base sur le calcul
d'un noyau 2 couches a densite constante. Bon comportement. 
* lecture/ecriture des modeles: on utilise une nouvelle subroutine
rwbin.f Ceci est fait dans CEPAM plutot que dans resoutp. 
* Lecture des namelists dans lit_nl.

Avril, Mai 2003: 
********* Passage a CEPAM3 *************
* CEPAM3.0.0
* Modification des fichiers d'entree .don. Introduction des equations
d'etat des elements lourds.

12/01/2003: modif de lim_noyau (appel a equation pour dro/dr), de figures.

09/01/2003: petite modif dans resoutp (definition de modele)
	    ptt modif dans staticpm (luminosite au centre -> 
		il faudrait modifier cette facon de faire)
	    actualisation de jacob_cep.f, optimise.f

07/01/2003: 
* VERSION 2.0.2
* Compilation de bin/jupiter a partir d'un fichier Makefile
* Ajout d'une EOS SESAME/Ross: etat_div.f

18/05/2000:
* fichier contrib4.dat: fonctions de contribution des moments gravitationels.
(dans figures4.f)

12/12/1999:
* limites de tableaux passes de 500 a 600 dans pp2, pp2d, pp2dd (etait a 60! dans pp2dd)
* RESOUTP: test sur le rayon inclus aussi pour dt>0

09/12/1999:	
* Ajout du sous-prog lim_xa3 qui utilise les dernieres conditions aux
limites d'apres Hubbard et al. (utilise le fichier
data/Limites/xgm5afil3pp.dat)
* Modif de lim_x pour afficher un message d'avertissement qd
gravite<100
* Ajout d'une inversion avec age=0 (dans cepam.f)
* Conditions d'arret de l'integration: dt<0 et R >= 6 Rjup (dans cepam.f)
* Fichiers d'input et output par defaut: zeus-ini et zeus (ds cepam.f)
* Changement de pas temporel qd dt<0 (ds updatep.f)