Changeset 3453 for LMDZ6/trunk/libf/phylmd/tracco2i_mod.F90

Timestamp:

Jan 28, 2019, 7:31:11 PM (6 years ago)

Author:

oboucher

Message:

Adding some diagnostics for type_trac=co2i

File:

• LMDZ6/trunk/libf/phylmd/tracco2i_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/tracco2i_mod.F90

@@ -12 +12 @@
     USE infotrac_phy
     USE geometry_mod, ONLY: cell_area
-    USE carbon_cycle_mod, ONLY: id_CO2, nbcf_in, fields_in, cfname_in, fco2_ocn_day, fco2_ff, fco2_bb
+    USE carbon_cycle_mod, ONLY: id_CO2, nbcf_in, fields_in, cfname_in
+    USE carbon_cycle_mod, ONLY: fco2_ocn_day, fco2_ff, fco2_bb, fco2_land, fco2_ocean
     USE carbon_cycle_mod, ONLY: carbon_cycle_tr, carbon_cycle_rad, RCO2_glo, RCO2_tot
     USE mod_grid_phy_lmdz
@@ -50 +51 @@
     INTEGER                        :: it, k, i, nb
     REAL, DIMENSION(klon,klev)     :: m_air     ! mass of air in every grid box [kg]
-    REAL, DIMENSION(klon)          :: co2land   ! surface land CO2 emissions [kg CO2/m2/s]
-    REAL, DIMENSION(klon)          :: co2ocean  ! surface ocean CO2 emissions [kg CO2/m2/s]
     REAL, DIMENSION(klon_glo,klev) :: co2_glo   ! variable temporaire sur la grille global
     REAL, DIMENSION(klon_glo,klev) :: m_air_glo ! variable temporaire sur la grille global
@@ -85 +84 @@
 
 !--retrieving land and ocean CO2 flux
-    co2land(:)=0.0
-    co2ocean(:)=0.0
+    fco2_land(:)=0.0
+    fco2_ocean(:)=0.0
     DO nb=1, nbcf_in
 !--fCO2_nep comes in unit of kg C m-2 s-1
 !--converting to kg CO2 m-2 s-1
-      IF (cfname_in(nb) == "fCO2_nep" )   co2land(:)=fields_in(:,nb)*RMCO2/RMC*pctsrf(:,is_ter)
+      IF (cfname_in(nb) == "fCO2_nep" )   fco2_land(:)=fields_in(:,nb)*RMCO2/RMC*pctsrf(:,is_ter)
 !--fCO2_fgco2 comes in unit of mol C02 m-2 s-1
 !--converting to kg CO2 m-2 s-1 + change sign
-      IF (cfname_in(nb) == "fCO2_fgco2" ) co2ocean(:)=-1.*fco2_ocn_day(:)*RMCO2/1.e3*(pctsrf(:,is_oce)+pctsrf(:,is_sic))
+      IF (cfname_in(nb) == "fCO2_fgco2" ) fco2_ocean(:)=-1.*fco2_ocn_day(:)*RMCO2/1.e3*(pctsrf(:,is_oce)+pctsrf(:,is_sic))
     ENDDO
 
-!--preparing the net anthropogenic flux at the surface for mixing layer 
+!--preparing the net anthropogenic flux at the surface for mixing layer
 !--unit kg CO2 / m2 / s
-    source(:,id_CO2)=fco2_ff(:)+fco2_bb(:)+co2land(:)+co2ocean(:)
+    source(:,id_CO2)=fco2_ff(:)+fco2_bb(:)+fco2_land(:)+fco2_ocean(:)
 
 !--computing global mean CO2 for radiation