Problemes de stabilité numérique :
    On passe à Ri=0 dans le calcul de Kz (zri=0., dans vdif_k.F90)
    On ajoute un facteur z/(l0+z) devant le calcul de l'entrainement de thermcell_2002.F90 
          le deuxieme changement est davantage pour eviter d'avoir quelque chose
          de trop grand dans la couche de surface


Radiatif/convectif 1D monte jusqu'à 8km.
Mais 2km en 3D.


Quelles sorties pour le rayonnement pur ?
Espoir : calculer DTLW + DTSW sur une atmosphère initialement uniforme en latitude.
Mais ne marche pas à cause du fait que l'atmophère absorbe peu.
Tout passe par la surface ...
Comment faire un calcul qui rende compte de ca ?
c'est à dire du forcage de la circulation en latitude ?
PPeut être pas possible

Rayonnement + diffusion OK avec day_step=240 iphysiq=10 (pas de temps de une heure)
pour 39 niveaux

Rayonnement + diffusion + thermiques en L39 :
plante avec pas de temps de 1heure
OK avec day_step=720 iphysiq=5 -> pas de temps de 10min


Introduction de la conservation de l'énergie.
Sinon erreurs de qq W/m2
