Mise en place sur kami du modele couple Predicate
Arpege Climat V3 Cycle 22a+
1.Compilation du modele
Les repertoires auxquels renvoient ces notes peuvent etre trouves en
local, via cvs, avec CVSROOT=/home/evian/ROOTCVS et sur kami (sous /u/cf/cglo/cglo765/PREDICATE/modipsl/modeles).
Arpege se compose de librairies (ARPEGE/src/library)
-
master[numero_du_cycle].x.a, contenant les .o de toutes les routines du
cycle de base (ici, le cycle 22). Les sources de cette librairie sont sous:
ARPEGE/src/master.
-
dummy[numero_du_cycle].a
-
xrd[numero_du_cycle].x.a
Viennent s'y ajouter des modules propres au cycle de base (module22.m.tar),
des routines de couplage Oasis (couple22.tar) et de .h (header22.tar).
La compilation de ses diverses routines et librairies permet d'obtenir
l'executable du cycle 22.
Pour obtenir celui du 22a, on se contente de compiler en plus les routines
du repertoire ARPEGE/src/22a. Les .o ainsi obtenus remplacent les .o de
la librairie "master" (ils sont sauves dans un tar: lib22a.x.tar).
Pour obtenir celui du 22b, on copie les routines de ARPEGE/src/22a
dans un repertoire de travail. On y copie ensuite les routines du repertoire
ARPEGE/src/22a+ (ce qui peut ecraser certaines des routines de la 22a).
On compile cette combinaison de sources pour avoir la librairie de la 22a+
lib22a+.x.tar.
L'ensemble de ces operations est effectue dans le script ARPEGE/script/compile22a+.
On peut chosir de recompiler l'ensemble des sources (22a et 22a+) ou bien
seulement celles de la 22a+, ou bien juste certaines routines, a preciser
par l'utilisateur.
2. Les adaptations au couple Predicate
-> modifications du code: cvs
-> modifications d'Arpextract: adaptation au netcdf
La deuxieme modification d'importatnce dans le post-traitement concerne
le choix des variables qui feront l'objet de moyennes mensuelles.
Dans la liste
complete, signalons la presence des variables de sol profond et des reservoirs.
En outre, les champs:
> NEBUL.CONVEC
> NEBUL.HAUTE
> NEBUL.BASSE
> NEBUL.MOYENN
ne sont pas disponibles dans nos simulations, la clef LNEBPAR de la
namelist n'ayant pas ete mise a .T..
3. Les problemes resolus
-
L'executable arpege (ou un executable) ne se lance pas (Message d'erreur:
EXEC
2 failed: Device busy) : Oasis a ete compile pour tourner en parallele
(-Wl,-P). Or kami, ne peut pas lancer deux executable parallele sur le
meme processuer (a l'aide la commande c "exec" contenu dans l'executable
de lancement "lance".
-
La date du restart n'est pas correcte: on utilise le petit utilitaire modif_date_restart_kami.
-
les semaphores et segments de memoire partagee ne sont pas liberes correctement
(sur kami, lancer sur le processeur principal le script rmtool pour effacer
les segments et egalement
"rsh cbar1 -n rmtool" sur le processeur interactif pour effacer les semaphores,
comme on fait "rsh vpp700 -n rmtool" au centre europeen).
-
La comparaison du champ
de temperature tel qu'il est passe par Oasis
avec la sortie Arpege SISUTESU
fait apparaitre un probleme. Des structures de temperatures continentales
apparaissent au milieu des oceans. En fait, la routine updclpdm lisait
les donnees d'orographie avec un decalage de 20 points, et replacait les
temperatures continales sur le champ de couplage en reportant ce decalage.
-
Le seuillage entre -70 et -10 W/m2/K du champ dQ/dT (derivee du flux non
solaire par rapport a la temperature) dans la routine wrcpldm.F est indispensable
pour eviter de multiplier le flux non solaire par une valeur trop importante
ou trop petite lors de la re-repartition du flux entre surface englacee
et surface d'eau libre.
-
Calendrier: le calcul du nombre de jours dans l'annee est fait dans la
routine UPDCAL. Ce calcul tient compte des annees multiples de 4 (bissextiles)
et de 100 (non bissextiles) mais PAS des annees multiples de 400. Bref,
on devrait avoir un probleme au passage de l'an 2000 !
-
Bcond: afin de pouvoir nous servir de notre run Predicate comme run de
controle d'un futur run de scenario (CO2), les fichiers de conditions aux
limites doivent etre crees avec des concentrations d'aerosols sulfates
de l'ere pre-industrielle. On ne dispose que des donnees de 1950. On modifie
donc les bcond a l'aide du petit programme ini_cli_aer.22a
utilisant la namelist conf923n.22
et les fichiers de Michel Deque situes sur kami dans le repertoire /workgc/mrga/mrga561.
Le fichier
ASCII qui contient les valeurs des aerosols sulfates anthropiques est
rempli de zero.En outre, on modifie la namelist initiale pour initialiser
les aerosols des gaz a effet de serre aux valeurs pre-industrielles suivant
les valeurs donnees par le Hadley Center pour leur simulation HTCM22
(changements climatiques anthropiques). On modifie egalement, dans cette
meme namelist, la valeur de la constante solaire (pour prendre celle de
1850).
|
|