Changes between Version 36 and Version 37 of IPSLCM5ACo2

Timestamp:

02/21/11 22:29:37 (13 years ago)

Author:

cetlod

Comment:

--

IPSLCM5ACo2

@@ -23 +23 @@
  1. on doit rajouter la lecture d'un fichier d'émissions ff (fossil fuel),
  1. on doit récupérer les flux moyens de PISCES et d'ORCHIDEE,
- 1. Ã  la place de la lecture des concentrations IPCC_CO2 imposées, on va calculer la concentration moyenne mensuelle
+ 1. A la place de la lecture des concentrations IPCC_CO2 imposées, on va calculer la concentration moyenne mensuelle
    (ou annuelle).
  1. on doit conserver ces valeurs de co2 (une ligne par mois)
@@ -63 +63 @@
 }}}
 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). [[BR]]
+On va devoir modifier PISCES pour afficher ce CO2 calculé tous les mois ( A faire). [[BR]]
 MM: récupération de ces affichages faite par lmdz_analyse_opa_out.awk.[[BR]]
     voir /work/cont003/p86manci/IPSLCM5A-C/config/ pour les tests.
@@ -139 +139 @@
 On doit alors récupérer chaque mois le co2_prec=co2 du mois précédent pour appliquer la formule.
 
-La condition initiale pour le co2_prec est un nouveau paramètre dans la lmdz.card : lmdz_UserChoices_co2_init.
+La condition initiale pour le co2_prec est un nouveau paramè&tre dans la lmdz.card : lmdz_UserChoices_co2_init.
 
 === Tests et runs ===
@@ -170 +170 @@
 [[BR]]
 
-Pour la version IPSLCM5A avec un CO2 atmosphèrique moyen ( calculé offline ), '''il ne faut pas activer la clé key_cpl_carbon_cycle''' ( l'océan s'attend à  recevoir dans ce cas une pco2 "spatialisée" ) (CE)
+Pour la version IPSLCM5A avec un CO2 atmosphérique moyen ( calculé offline ), '''il ne faut pas activer la clé key_cpl_carbon_cycle''' ( l'océan s'attend à  recevoir dans ce cas une pco2 "spatialisée" ) (CE)
  * PISCES est modifié pour diagnostiquer le flux total de carbone océanique cumulé sur la durée du run
 {{{
@@ -210 +210 @@
 Test historical interactif hist02 : voir [[BR]]
 /work/cont003/p86manci/IPSLCM5A-C/config/IPSLCM5A_C/Chistorical2 [[BR]]
-Avec le script commité dans IPSLCM5A/Chistorical, on peut restarter correctement du run de control.
+Avec le script commité dans IPSLCM5A/Chistorical, on peut restarter correctement du run de control.
  * voir [[BR]] :
    http://dods.ipsl.jussieu.fr/orchidee/SANORCHIDEE/CpiControl/fegg_plot01_Cpi02Chist02_prod/ [[BR]]
@@ -282 +282 @@
 CMIP5 ... c'est bon donc ?!
 }}}
+
+=== Janvier 2011 : Christian Ethé ===
+
+La formule pour le calcul de la pco2 est la suivante
+{{{
+    co2_ppm = co2_ppm_prec + ( EMI(Fossil) + Flux(ocean) + Total Flux(land) ) / 2.12
+     ppmv      ppmv             GtC/mth        GtC/mth       GtC/mth
+}}}
+
+Dans ORCHIDEE  
+{{{  
+   Total Flux(land) = Flux(land) + EMI(Land-Use) +    HARVEST
+       NBP          =   (NEP)    +  (FLUC)       + (HARVEST_ABOVE ) 
+}}} 
+
+Les craintes de MM sur LAND COVER CHANGE étaient justifiées. Il faut intégrer ces flux dans la routine stomate.f90 comme c'est fait pour le flux NEP :
+{{{
+!         Total ( land) Cumulative fluxes of fluc and harvest
+          net_cflux_prod_monthly_sum=&
+              &  SUM(cflux_prod_monthly(:)*resolution(:,1)*resolution(:,2)*contfrac(:))*1e-15
+          CALL reduce_sum(net_cflux_prod_monthly_sum,net_cflux_prod_monthly_tot)
+          CALL bcast(net_cflux_prod_monthly_tot)
+
+          net_harvest_above_monthly_sum=&
+             &   SUM(harvest_above_monthly(:)*resolution(:,1)*resolution(:,2)*contfrac(:))*1e-15
+          CALL reduce_sum(net_harvest_above_monthly_sum,net_harvest_above_monthly_tot)
+          CALL bcast(net_harvest_above_monthly_tot)
+
+          net_co2_flux_monthly = net_co2_flux_monthly*1e-15
+          CALL reduce_sum(net_co2_flux_monthly,net_co2_flux_monthly_sum)
+          CALL bcast(net_co2_flux_monthly_sum)
+
+          WRITE(numout,9010) 'GLOBAL net_cflux_prod_monthly    (Peta gC/month)  = ',net_cflux_prod_monthly_tot
+          WRITE(numout,9010) 'GLOBAL net_harvest_above_monthly (Peta gC/month)  = ',net_harvest_above_monthly_tot
+          WRITE(numout,9010) 'GLOBAL net_co2_flux_monthly      (Peta gC/month)  = ',net_co2_flux_monthly_sum
+
+!         Calculation of net biospheric production
+          net_biosp_prod_monthly_tot =  &
+             &    ( net_co2_flux_monthly_sum + net_cflux_prod_monthly_tot + net_harvest_above_monthly_tot )
+          WRITE(numout,9010) 'GLOBAL net_biosp_prod_monthly    (Peta gC/month)  = ',net_biosp_prod_monthly_tot
+
+9010  FORMAT(A52,F17.14)
+}}}
+
+=== fichier fuel ===
+
+Le fichier "emissions" preparé par Patricia est la :
+http://dods.ipsl.jussieu.fr/cpipsl/ANDRES/CMIP5_gridcar_CO2_emissions_fossil_fuel_Andres_1751-2007_monthly_SC.txt
+Dans ce fichier, FF est déjà en PgC/month
+