Equilibre radiatif-convectif
----------------------------

Possibilités :
 + Temperature de surface imposée, ou modèle thermique de sol (en mettant 
soil_model=n ou y dans gcm.def) 
 + surface océanique ou continentale (natsurf dans lmdz1d.def)
 + beta (=evap/[evap. pot.]) interactif (bucket) ou imposé (en mettant 
qsol0>0 dans physiq.def, ce qui impose qsol = qsol0)
 + refroidissement radiatif imposé (iflag_radia=0 dans physiq.def, et 
derniere colonne de lscale.inp.001 = cooling rate en K/s) ou interactif 
(iflag_radia=1 dans physiq.def et derniere colonne de lscale.inp.001 =0).

Options du SCM et nature de la simulation définies par :
   lmdz1d.def
   run.def (durée de la simul = nday) 
   gcm.def
   physiq.def

Etat initial :
   prof.inp.001

Conditions de grande échelle (vent géostrophique ... ; ici zéro) :
   lscale.inp.001

Pas de forçages (convergences de chaleur et d'humidité)

Pas de forçage du soulèvement au début :

-------------------------------------------------------------------------

execution du SCM :
  soit lmdz1d.e > lmdz1d.lis
  soit xqt.x 
Le script d'execution xqt.x permet de renommer les fichiers de sortie (par 
exemple, "xqt.x 1" copiera le fichier histmth.nc dans histmth_1.nc) et 
modifie ces fichiers de façon que l'axe vertical soit orienté avec les 
pressions croissant vers le bas et que la valeur "indéfinie" soit 10^35 (et 
non zero, comme dams le fichier originel).

------------------------------------------------------------------------

Autre script qui peut etre utile : compile.x permet de compiler le SCM et 
de renommer l'éxécutable.