Changes between Version 32 and Version 33 of IPSLCM5ACo2

Timestamp:

01/01/11 01:36:09 (13 years ago)

Author:

mmaipsl

Comment:

--

IPSLCM5ACo2

@@ -215 +215 @@
    pour la comparaison de cette simulation Chist02 avec le control CO2 interactif Cpi02.
 
+=== Alerte pour le LAND COVER CHANGE ===
+Le flux de LAND COVER CHANGE est calculé une fois par an dans ORCHIDEE. Si c'est un flux annuel, alors
+on doit peut-être le divisé par 12 dans l'affichage du modèle.
+
+Voir les deux courriels :
+{{{
+Bonjour Nicolas (Laurent en cc)
+
+J'ai ajouté dans la version 1.9.5.1 (cf
+http://forge.ipsl.jussieu.fr/orchidee/viewvc/ORCHIDEE/src_stomate/stomate_lpj.f90?view=markup&revision=1.29&root=orchidee&pathrev=orchidee_1_9_5_1
+), pour le cycle du carbone interactif que je
+mets en place, l'affichage du flux de carbone total du land use chaque
+mois. Tu peux voir le lien
+https://forge.ipsl.jussieu.fr/igcmg/wiki/IPSLCM5ACo2
+pour les infos sur ce couplé carbone.
+
+Les trois flux convflux(:) + cflux_prod10(:) + cflux_prod100(:) sont
+restartés et sont donc constants toute l'année.
+Je pense que ce sont des flux par an et comme je récupère tous les
+flux par mois, il me semble que je dois (dans le
+code) les divisés par 12... Tu peux voir la sortie co2 dans ce fichier :
+mercure.ccc.cea.fr
+/work/cont003/p86manci/IPSLCM5A-C/config/IPSLCM5A_C/Chistorical2/co2.log
+}}}
+
+{{{
+Bonsoir
+
+je vous écris après mes séances d'hosto .. sachant que je n'aurai
+probablement qu'une heure ou deux
+lundi pour résoudre ce problème essentiel pour CMIP5 !
+Voici les formules de stomate_lpj concernant le flux de LAND COVER
+CHANGE :
+    cflux_prod_total(:) = convflux(:) + cflux_prod10(:) + cflux_prod100(:)
+    IF ( EndOfMonth ) THEN
+      
+cflux_prod_total_sum=SUM(cflux_prod_total(:)*resolution(:,1)*resolution(:,2)*contfrac(:))*1e-15
+       CALL reduce_sum(cflux_prod_total_sum,cflux_prod_sum)
+       CALL bcast(cflux_prod_sum)
+       WRITE(numout,*) 'net cflux_prod_total_sum (Peta gC/month)  =
+',cflux_prod_total_sum
+       WRITE(numout,*) 'GLOBAL net cflux_prod_total_sum (Peta
+gC/month) (Peta gC/month)  = ',cflux_prod_sum
+
+On affiche donc toute l'année ce même flux... (les flux de LAND COVER
+CHANGE n'étant modifiés qu'une fois par an)
+
+[...] puis pour l'écriture du flux flux CMIP5 :
+       vartmp(:)=cflux_prod_total/1e3/one_day*contfrac
+       CALL histwrite (hist_id_stomate_IPCC, "fLuc", itime, &
+         vartmp, npts, hori_index)
+
+Avec la définition dans intersurf de cette sortie qui lui donne les
+unités kg C m-2 s-1 ...
+    ! CO2 Flux to Atmosphere from Land Use Change
+    CALL histdef (hist_id_stom_IPCC, &
+         &               TRIM("fLuc"), &
+         &               TRIM("CO2 Flux to Atmosphere from Land Use
+Change"), &
+         &               TRIM("kg C m-2 s-1"), iim,jjm, hist_hori_id, &
+         &               1,1,1, -99,32, ave(1), dt, hist_dt)
+
+Ais-je bien affiché le flux moyen mensuel ? par seconde ? (plutôt par
+jour, non ?)
+L'autre flux NEP a les mêmes unités et le même calcul pour la sortie
+CMIP5 ... c'est bon donc ?!
+}}}