1 | septembre 2009. |
---|
2 | -------------- |
---|
3 | |
---|
4 | fichiers modifies : |
---|
5 | |
---|
6 | M Antarctique_general_files/lect-Ant_CISM_gen_dat.f90 (cosmetique) |
---|
7 | A resol_adv_diff_2D-sept2009.f90 (changement de place de certaines declarations) |
---|
8 | M Makefile (pour compiler avec ifort 10) |
---|
9 | M Ant15_CISM_files/module_choix_antar15CISM.f90 (use dragging_vit_bil_CISM_gen) |
---|
10 | M Fichiers-parametres/CISMA20_param_list.dat |
---|
11 | M Fichiers-parametres/CISMA15_param_list.dat |
---|
12 | M Fichiers-parametres/CISMA25_param_list.dat |
---|
13 | A conserv-mass-adv-diff_sept2009_mod.f90 (idem le precedent "juin"). prescription de hauteur |
---|
14 | A litho-0.4.f90 (version utilisant les operations tableaux de f90) |
---|
15 | M iso_declar_mod-0.3.f90 (cosmetique) |
---|
16 | |
---|
17 | |
---|
18 | |
---|
19 | Quelques explications : |
---|
20 | --------------------------- |
---|
21 | |
---|
22 | Modifications pour essayer de faire tourner Grisli sur une plus petite grille. |
---|
23 | |
---|
24 | Au depart pour les grilles < 25 km : |
---|
25 | |
---|
26 | memory fault au premier appel |
---|
27 | call resol_adv_diff_2D_vect(Dmx,Dmy,advmx,advmy,Hp,i_Hp,bilmass,vieuxH,H) |
---|
28 | |
---|
29 | D'aprs Michel c'est une histoire de memory stack. |
---|
30 | |
---|
31 | En effet, en faisant ulimit -s 65000 qui change la mmemoire allouee au stack ca tourne. |
---|
32 | |
---|
33 | Ca peut provenir des tableaux automatiques : les tableaux de travail |
---|
34 | qui sont declares dans une routine mais en mettant ceux de |
---|
35 | resol_adv_diff_2D_vect dans le module avant contains (du coup ils |
---|
36 | deviennent globaux) ca ne change rien. |
---|
37 | |
---|
38 | |
---|
39 | Test des changements de compilateur ifort9 -> ifort10 |
---|
40 | ------------------------------------------------------- |
---|
41 | |
---|
42 | - motivation. Le probleme de fonctionnement aux petites grilles vient |
---|
43 | du memory stack. A partir de la version 10, il y a un ordre |
---|
44 | -heap-array qui permet de gerer differemment cette memoire. |
---|
45 | |
---|
46 | |
---|
47 | - Probleme additionnel, en version 10 la routine litho-0.3 plante, |
---|
48 | modification (litho-0.4) avec une syntaxe f90. Il faut verifier que ça |
---|
49 | donne les mêmes resultats. |
---|
50 | |
---|
51 | |
---|
52 | Tests en 25 km |
---|
53 | |
---|
54 | Pour tourner en ifort9, il faut modifier l'appel des librairies. |
---|
55 | export LD_LIBRARY_PATH=/usr/localcat/intel/compiler9/lib:/usr/local/intel/mkl9/lib/32/:$LD_LIBRARY_PATH |
---|
56 | |
---|
57 | A25tst01 : ifort9 litho-0.3 run 10000 ans depuis un stationnaire. |
---|
58 | A25tst02 : ifort9 litho-0.4 run 10000 ans depuis un stationnaire. |
---|
59 | |
---|
60 | Resultat, les differences n'excedent pas 1 m mais plutôt de l'ordre du mm sauf quelques regions cm-dm. |
---|
61 | Conclusion : On garde litho-0.4 |
---|
62 | |
---|
63 | A25tst03 : ifort10 litho-0.4 run 10000 ans depuis un stationnaire. |
---|
64 | avec l'option -heap-array |
---|
65 | |
---|
66 | De differences en moyenne assez faibles mais sur certains ice streams |
---|
67 | peuvent être fortes (max quelques m). D'apres Michel S. ce n'est pas |
---|
68 | normal compte tenu des options de compil qui devraient donner des |
---|
69 | resultats parfaits. Les differences pourraient venir de mkl dont on ne |
---|
70 | sait pas comment ils ont ete compiles. |
---|