Pour faire tourner le modèle en mode grounding line retreat
----------------------------------------------------------------
Pour l'instant, on modifie la définition de time_gl dans precribe dans routine 
init_prescribe_retreat.



Pour faire fondre les ice_shelves: début a t_break dans conservation de la masse 
----------------------------------


Expériences AGU
-----------------
run 16  +100 m (?, sans doute dans la correction beta)
run 17  +200 m
run 18 correction de dS en lisant le beta du run 17
run 25 GL
run 26 melting shelves
Run 27 gl + shelves (plante)
Run 28 idem 27 (fusion +recul)
Run 29 idem
run 30 fusion seule
run 32 break, n'a pas l'air d'avoir marché

run 45 : idem run 25 retreat gl
run 47 : idem run 27 retreat gl + melting ice shelves
run 48 : idem run 27  retreat gl 500 m /an

Warnings : a corriger éventuellement
-------------------------------------
- les modes grounding line retreat, ... ne sont pas dans le fichier param 
- Dans icethick, des scothc sur hdot (10 m/an) et sur les vitesses (3000 m/an)
- noeuf flottants sf*0.4
- bug dans beta
- pourquoi casser les ice shelves donne un effet inverse a celui de les ramolir (effet pente ?)


Idées en vrac
---------------
- routine "outcrops" Quand la vitesse est trop grande dans les régions de montagne, répartir directement aux noeuds en dessous.

- regarder dTb/dH entre plusieurs cartes de spinup



Calendrier
--------------
- Couper step
- refaire base de donnee avec nouveau CISM.
  * Regarder flottaison and co.
  * Regarder le probleme des fleuves de glace
- refaire vitesses de bilan avec donnees Rignot.