POLCA: le couplé Predicate










                                        Mise à jour: 4 fevrier 2002
 
 
 

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).