septembre 2009.
--------------

fichiers modifies : 

M      Antarctique_general_files/lect-Ant_CISM_gen_dat.f90    (cosmetique)
A      resol_adv_diff_2D-sept2009.f90 (changement de place de certaines declarations)
M      Makefile     (pour compiler avec ifort 10)
M      Ant15_CISM_files/module_choix_antar15CISM.f90 (use dragging_vit_bil_CISM_gen)
M      Fichiers-parametres/CISMA20_param_list.dat
M      Fichiers-parametres/CISMA15_param_list.dat
M      Fichiers-parametres/CISMA25_param_list.dat
A      conserv-mass-adv-diff_sept2009_mod.f90 (idem le precedent "juin"). prescription de hauteur 
A      litho-0.4.f90   (version utilisant les operations tableaux de f90)
M      iso_declar_mod-0.3.f90 (cosmetique)



Quelques explications : 
---------------------------

Modifications pour essayer de faire tourner Grisli sur une plus petite grille.

Au depart pour les grilles < 25 km :

memory fault au premier appel 
call resol_adv_diff_2D_vect(Dmx,Dmy,advmx,advmy,Hp,i_Hp,bilmass,vieuxH,H)

D'aprs Michel c'est une histoire de memory stack.

En effet, en faisant ulimit -s 65000  qui change la mmemoire allouee au stack ca tourne.

Ca peut provenir des tableaux automatiques : les tableaux de travail
qui sont declares dans une routine mais en mettant ceux de
resol_adv_diff_2D_vect dans le module avant contains (du coup ils
deviennent globaux)  ca ne change rien.


Test des changements de compilateur ifort9 -> ifort10
-------------------------------------------------------

- motivation. Le probleme de fonctionnement aux petites grilles vient
du memory stack. A partir de la version 10, il y a un ordre
-heap-array qui permet de gerer differemment cette memoire.


- Probleme additionnel, en version 10 la routine litho-0.3 plante,
modification (litho-0.4) avec une syntaxe f90. Il faut verifier que ça
donne les mêmes resultats.


Tests en 25 km

Pour tourner en ifort9, il faut modifier l'appel des librairies.
export LD_LIBRARY_PATH=/usr/localcat/intel/compiler9/lib:/usr/local/intel/mkl9/lib/32/:$LD_LIBRARY_PATH 

A25tst01 : ifort9  litho-0.3 run 10000 ans depuis un stationnaire.
A25tst02 : ifort9  litho-0.4 run 10000 ans depuis un stationnaire.

Resultat, les differences n'excedent pas 1 m mais plutôt de l'ordre du mm sauf quelques regions cm-dm.
Conclusion : On garde litho-0.4

A25tst03 : ifort10  litho-0.4 run 10000 ans depuis un stationnaire.
avec l'option -heap-array

De differences en moyenne assez faibles mais sur certains ice streams
peuvent être fortes (max quelques m). D'apres Michel S. ce n'est pas
normal compte tenu des options de compil qui devraient donner des
resultats parfaits. Les differences pourraient venir de mkl dont on ne
sait pas comment ils ont ete compiles.