source: trunk/mesoscale/NOTES.txt @ 91

Last change on this file since 91 was 91, checked in by aslmd, 14 years ago

mars: test outliers [dans initracer.F, commente] LMD_MM_MARS: modifications mineures [retrocompatibilite ancienne physique, callphys.def test pour nouvelle physique] PLOT: generalisation de la routine map_uvt pour pouvoir tracer des figures en projection polaire complete

File size: 6.4 KB
Line 
1
2concat.e puis localtime.e puis
3ncwa -O -v mtot,icetot -a longitude -d longitude,-179.0,179.0 diagfi.nc yeye.nc
4ncwa -O -v mtot -a longitude -d longitude,-180.0,180.0 concat_LT.nc mawd.nc
5
6        A FAIRE:::: mettre des flags precompilo dans les meso_
7        les reporter dans makegcm
8
9changer le renormalisation dans aeropacity ????
10on ne laisse pas aerosol comme le lifting veut qu'il soit !
11tenter des taux de soulevement pour que taudust_tmp soit les obs
12en prescivant une dust bomb fixe d opacite, on aura au moins la structure verticale
13
14        tester traceurs radiativement actifs avec la nouvelle physique ?????
15
16A FAIRE: PB LES sur iDATAPLEX (les points HFX nuls) (pas de soucis sur ciclad)
17METTRE SUR LE svn LA BASE d'ETATS INITIAUX ????
18
19more than 4 procs w/ nest ??? y reflechir
20        -----------------------------------------------------------------------
21        -- si possible comment determiner taille ?
22        nproc doit diviser e_we-1 (1er nest)
23        grid_ratio doit diviser e_we-1 +4 (1er nest)
24        soit e_we=ye+1
25        grid_ratio divise ye+4 et nproc divise ye
26        soit nproc=8, ye=8*i
27        ainsi il existe j tel que 8i + 4 = 3j ou encore 4*[2i+1] = 3j
28        verifie par exemple si 2i+1 est multiple de 3
29        il suffit donc de trouver un multiple impair de 3 et de deduire i
30        par exemple 2i+1=33 >>>> i=16
31        >>>> e_we = 129 pour le 1er nest (et ajouter 4 pour les suivants)
32        ------------------------------------------------------------------------
33
34        ne pas utiliser le FASTCASE avec traceurs (instabilites en haut)
35            ces instabilites sont cependant reglees si on augmente radt a 10 par exemple
36
37        pour le cycle de l'eau c'est OK de regler caps=F dans le mesoscale
38        sauf si on commence a devoiler la calotte permanente nord
39        ---> corrige, scenario caps specifique au mesoscale
40
41        NE SERAIT-CE PAS MIEUX DE TOUT TRANSMETTRE AUX BORNES ???
42        tous les traceurs, pas seulement vapor
43
44
45        - attention il faut les trois MARS sinon il s arrete sans message clair
46        - attention a ne pas lancer le modele s il est deja lance
47        - important que pd_scalar soit a T ... le mettre par defaut ????
48
49
50ROUTINES a AJOUTER sont dans COMMON_GCM
51- passer aux nouveaux makegcm [en commun avec Ehouarn si on veut le nouveau
52  readtesassim qui est en F90]
53- il faut tester le nest pour verifier les lignes trop longues
54
55        (ok) lier gr_fi_dyn qui est dans dyn3d
56        (ok) regler le pb du nouveau readtesassim (ou alors le lier tout simplement ou
57          l'appeler meso_readtesassim)
58        (ok) regler le pb meso_dustlift (le lier dans makemeso comme point precedent)
59             (car le souci c que dustlift est appele dans vdifc)
60
61        RESTE a ADAPTER le LES a la NOUVELLE PHYSIQUE
62        il y a normalement peu a faire
63        reste a faire egalement le -DNEWPHYS pour le LES
64
65        attention pb d'affichage des valeurs dans le fichier texte avec LES ???
66        bien que les valeurs du fichier soient tout a fait raisonnables
67        ... n'est-ce pas un effet de bord cache ????
68
69
70        apres fusion, le LES est reconnu par module_lmd_driver lorsque diff_opt=2 km_opt=2
71
72
73        -attention PB si on ne sort pas HFX et USTM (note dans le Registry)
74        -il faut run.def nouvelle physique [c est meme ce qui est utilise par runmeso]
75        - IL FAUT SE PENCHER SUR LE FAIT QU'ON INDIQUE q2val=0 dans lmd_driver ....
76
77-----------------------
78ATTENTION NOUVELLE PHYSIQUE
79Oui, c'est quelque chose qu'il faut probablement changer partout
80maintenant que la version de pgf90 à changé (sur les machines du LMD).
81Avec cette nouvelle version (7.1-6), le '-fast' est plus agressif
82qu'avant (et inclue entre autre les horribles '-Mvect=sse -Mscalarsse'
83qui dégradent la précision de certains calculs pour accélérer le code);
84je préconise de ne plus s'en servir. Bon d'accord, je n'ai pas fait une
85étude approfondie de l'impact de '-fast', mais j'ai vu qu'avec,
86j'obtenais des résultats différents lorsque je changeais simplement
87l'ordre des traceurs...
88
89Aymeric Spiga wrote:
90> je détecte ces changements d'option de compilation ; ont-ils de
91> l'importance ?
92>
93> Aymeric
94>
95> < #   set optim90=" -fast"
96> < #   set optimtru90=" -fast -c -Mfree "
97> < #   set optim90=" -O2 -Munroll=c:1 -Mnoframe -Mcache_align"
98> < #   set optimtru90=" -O2 -Munroll=c:1 -Mnoframe -Mcache_align"
99> <    set optim90=" -O2 -Munroll -Mcache_align"
100> <    set optimtru90=" -O2 -Munroll -Mcache_align"
101> ---
102>   
103>>    set optim90=" -fast"
104>>    set optimtru90=" -fast -c -Mfree "
105------------------------------
106
107
108        - attention a cp et R, normaliser une bonne fois pour toutes
109        - il manque sur le SVN les cas idealises
110- il manque sur le SVN les scripts MPI
111        - il faut recompiler les librairies NETCDF
112        - mettre la nouvelle physique
113        - mettre les DEF du meso-echelle
114
115        - modele ok sur auric
116- modele pas ok sur ciclad avec pgf2010, erreur inedite un seul module manquant
117        - modele LES OK sur ciclad
118        - modele LES ok sur auric
119
120        24/01/2011
121        tests g95 en 64bits natif sur systeme Linux
122        -- modifications de makemeso, tests
123        -- tout est OK sauf les libraires NETCDF, probleme d'underscore
124        -- OK avec libraires maison compilees avec g95 standard sur flores [et tourne OK]
125
126
127
128        mpi_64_pgf7_ncdf4_mpi1.2.txt
129        - probleme lors de la compilation de solve_em : LINUX runs out of memory [huchard]
130        - IL FAUT COMPILER SUR auric
131        nougaro est lent a la compilation, utiliser surtout auric
132
133
134
135
136______________________________________________________
137
138
139PB MPI
140/donnees/aslmd/MODELES/MPI/mpich2-1.2.1p1_PGF7/lib/libmpich.a(simple_pmi.o):
141In function `PMI_Init':
142simple_pmi.c:(.text+0x15c0): warning: Using 'gethostbyname' in statically
143linked applications requires at runtime the shared libraries from the glibc
144version used for linking
145/donnees/aslmd/MODELES/MPI/mpich2-1.2.1p1_PGF7/lib/libmpich.a(simple_pmi.o):
146In function `PMI_Init':
147simple_pmi.c:(.text+0x15c0): warning: Using 'gethostbyname' in statically
148linked applications requires at runtime the shared libraries from the glibc
149version used for linking
150/donnees/aslmd/MODELES/MPI/mpich2-1.2.1p1_PGF7/lib/libmpich.a(simple_pmi.o):
151In function `PMI_Init':
152simple_pmi.c:(.text+0x15c0): warning: Using 'gethostbyname' in statically
153linked applications requires at runtime the shared libraries from the glibc
154version used for linking
155/donnees/aslmd/MODELES/MPI/mpich2-1.2.1p1_PGF7/lib/libmpich.a(simple_pmi.o):
156In function `PMI_Init':
157simple_pmi.c:(.text+0x15c0): warning: Using 'gethostbyname' in statically
158linked applications requires at runtime the shared libraries from the glibc
159version used for linking
160
161
162POSSIBLE mars.sed
163
164s+ *../frame/module_internal_header_util.o ../frame/pack_utils.o
165-L../external/esmf_time_f90 -lesmf_time+& -L../mars_lmd/libo -llmd
166-Mmpi=mpich2+g
167
Note: See TracBrowser for help on using the repository browser.