IPSLCM5A et CO2 interactif

Mise en place d'une première version

réunion du 18/11/2010 : présent Laurent Bopp, Christian Ethé, Martial Mancip.

Deux versions de couplé avec CO2 interactif pour CMIP5 sont prévues :

• une première version avec la version IPSLCM5A avec un CO2 atmosphérique moyen (pas de transport ATM, flux océaniques et continentaux moyennés spatialement et temporellement) : comme IPSLCM4_LOOP
• une seconde version "IPSLCM5_ESM" avec les flux de carbone spatialisés et transportés en haute fréquence (jusqu'à la demi-heure pour la surface).

Nous établissons ici la liste des modifications à apporter au couplé IPSLCM5A (modèle et scripts) :

1. on doit rajouter la lecture d'un fichier d'émissions ff (fossil fuel),
2. on doit récupérer les flux moyens de PISCES et d'ORCHIDEE,
3. à la place de la lecture des concentrations IPCC_CO2 imposées, on va calculer la concentration moyenne mensuelle (ou annuelle).
4. on doit conserver ces valeurs de co2 (une ligne par mois)

   date |    co2 mensuel mois actuel | co2 mensuel mois précédent | ff | PISCES flux | ORCHIDEE NEP | ORCHIDEE LUC 
   

1. fichier fuel
2. modifications des modèles

Les flux d'émission de CO2 dans PISCES sont écrites dans ocean.output et si on rajoute la sortie tco2flux dans iodef.xml, il sera sauvé.

        WRITE(numout,*) 'Total Flux of Carbon out of the ocean :'
        WRITE(numout,*) '-------------------- : ',t_oce_co2_flx
        WRITE(numout,*) '(GtC/yr)'

Pour l'instant ce flux de CO2 global n'est sauvé qu'en moyenne mensuelle. L'affichage du flux est dans le ocean.output du mois de décembre. On va devoir modifier PISCES pour afficher ce CO2 calculé tous les mois (à faire).

Les flux d'émission de CO2 (NEP) dans ORCHIDEE sont écrits par processus MPI dans les sorties out_orchidee_000? :

 net_co2_flux_monthly (Peta gC/month)  =   

Dans ORCHIDEE, il faut modifier la sortie du flux de co2 cumulé mensuel et imprimer le flux de LAND USE co2_fluc (moyenné par mois).

On devra donc obligatoirement modifier ces deux modèles pour créer cette version du couplé CO2.

On ne peut pas passer par les fichiers NetCDF sur les machines SX. On doit donc passer par les fichiers textes.

3. Le calcul de la concentration mose fait dans