Changeset 3421

Timestamp:

Dec 7, 2018, 10:08:07 PM (5 years ago)

Author:

oboucher

Message:

Further changes for the interactive CO2 cycle
Pulling out the definition of emission terms from tracco2i into carbon cycle module
Declaration in phyetat0

Location:

LMDZ6/trunk/libf/phylmd

Files:

• carbon_cycle_mod.F90 (modified) (1 diff)
• phyetat0.F90 (modified) (1 diff)
• tracco2i_mod.F90 (modified) (8 diffs)

LMDZ6/trunk/libf/phylmd/carbon_cycle_mod.F90

@@ -74 +74 @@
   REAL, DIMENSION(:), ALLOCATABLE :: fco2_lu_day     ! Emission from land use change for 1 day (cumulated) [gC/m2/d]
 !$OMP THREADPRIVATE(fco2_lu_day)
+  REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_ff ! Emission from fossil fuel [kgCO2/m2/s]
+!$OMP THREADPRIVATE(fco2_ff)
+  REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_bb ! Emission from biomass burning [kgCO2/m2/s]
+!$OMP THREADPRIVATE(fco2_bb)
 
   REAL, DIMENSION(:,:), ALLOCATABLE :: dtr_add       ! Tracer concentration to be injected 

LMDZ6/trunk/libf/phylmd/phyetat0.F90

@@ -435 +435 @@
      END DO
      CALL traclmdz_from_restart(trs)
-
+  ENDIF
+
+!--OB now this is for co2i
+  IF (type_trac == 'co2i') THEN
      IF (carbon_cycle_cpl) THEN
         ALLOCATE(co2_send(klon), stat=ierr)

LMDZ6/trunk/libf/phylmd/tracco2i_mod.F90

@@ -12 +12 @@
     USE infotrac
     USE geometry_mod, ONLY : cell_area
+    USE carbon_cycle_mod, ONLY : nbcf_in, fields_in, cfname_in, fco2_ocn_day, fco2_ff, fco2_bb
     USE mod_grid_phy_lmdz
     USE mod_phys_lmdz_mpi_data, ONLY :  is_mpi_root
@@ -44 +45 @@
 !----------------
 
-    INTEGER, PARAMETER :: id_CO2=1                           !--temporaire OB=> PC to be changed 
-    REAL, PARAMETER    :: MCO2=44.011   !--g/mol
-    INTEGER                                 :: it, k, i
-    REAL, DIMENSION(klon,klev)              :: m_air          ! mass of air in every grid box [kg]
-    REAL, DIMENSION(klon)                   :: co2ff          ! surface fossil-fuel CO2 emissions [kg CO2/m2/s]
-    REAL, DIMENSION(klon)                   :: co2bb          ! surface biomass burning CO2 emissions [kg CO2/m2/s]
+    INTEGER, PARAMETER :: id_CO2=1              !--temporaire OB -- to be changed 
+    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
@@ -68 +68 @@
     IF (debutphy) THEN
       IF (MAXVAL(tr_seri(:,:,id_CO2)).LT.1.e-15) THEN 
-        tr_seri(:,:,id_CO2)=280.e-6/RMD*MCO2
+        tr_seri(:,:,id_CO2)=280.e-6/RMD*RMCO2
       ENDIF 
     ENDIF
 
-!--calculate mass of air in every grid box in kg 
+!--calculate mass of air in every grid box in kg air
     DO i=1, klon
     DO k=1, klev
@@ -81 +81 @@
 !--call CO2 emission routine 
 !--co2bb is zero for now
-    CALL co2_emissions(debutphy,co2ff,co2bb)
+!--unit kg CO2 m-2 s-1 
+    CALL co2_emissions(debutphy)
+
+!--retrieving land and ocean CO2 flux
+!--fCO2_nep comes in unit of g CO2 m-2 dt_stomate-1 
+!--this needs to be changed in ORCHIDEE
+    co2land(:)=0.0
+    co2ocean(:)=0.0
+    DO nb=1, nbcf_in
+      IF (cfname_in(nb) == "fCO2_nep" )   co2land(:)=fields_in(:,nb)*RMCO2/RMC/86400./1000.
+      !!IF (cfname_in(nb) == "fCO2_fgco2" ) co2ocean(:)=fco2_ocn_day(:) !--for now
+    ENDDO
 
 !--preparing the net anthropogenic flux at the surface for mixing layer 
 !--unit kg CO2 / m2 / s
-    source(:,id_CO2)=co2ff(:)+co2bb(:)
+    source(:,id_CO2)=fco2_ff(:)+fco2_bb(:)+co2land(:)+co2ocean(:)
 
 !--computing global mean CO2 for radiation
+!--every timestep for now but enough every month
 !    IF (debutphy.OR.mth_cur.NE.mth_pre) THEN
       CALL gather(tr_seri(:,:,id_CO2),co2_glo)
       CALL gather(m_air,m_air_glo)
 !$OMP MASTER
+!--conversion from kg CO2/kg air into ppm
        IF (is_mpi_root) THEN
-         RCO2_glo=SUM(co2_glo*m_air_glo)/SUM(m_air_glo)*1.e6*RMD/MCO2
+         RCO2_glo=SUM(co2_glo*m_air_glo)/SUM(m_air_glo)*1.e6*RMD/RMCO2
        ENDIF
-       PRINT *,'in tracco2i: global CO2 in ppm =', RCO2_glo
+       PRINT *,'toto in tracco2i: global CO2 in ppm =', RCO2_glo
 !$OMP END MASTER
        CALL bcast(RCO2_glo)
@@ -103 +116 @@
   END SUBROUTINE tracco2i
 
-  SUBROUTINE co2_emissions(debutphy,co2ff,co2bb)
+  SUBROUTINE co2_emissions(debutphy)
 
     USE dimphy
@@ -116 +129 @@
     USE netcdf, ONLY: nf90_get_var, nf90_noerr, nf90_nowrite
 
+    USE carbon_cycle_mod, ONLY : fco2_ff, fco2_bb
+
     IMPLICIT NONE
 
     INCLUDE "YOMCST.h"
     LOGICAL,INTENT(IN) :: debutphy
-    REAL,DIMENSION(klon),INTENT(out) :: co2ff            !  fossil-fuel CO2 emissions
-    REAL,DIMENSION(klon),INTENT(out) :: co2bb            !  biomass burning CO2 emissions
 
 ! For NetCDF:
@@ -130 +143 @@
     REAL, POINTER:: vector(:), time(:)
     REAL,ALLOCATABLE       :: flx_co2ff_glo(:,:) !  fossil-fuel CO2 
+    REAL,ALLOCATABLE       :: flx_co2bb_glo(:,:) !  biomass-burning CO2 
     REAL,ALLOCATABLE, SAVE :: flx_co2ff(:,:)     !  fossil-fuel CO2 
-    REAL,ALLOCATABLE       :: flx_co2bb_glo(:,:) !  biomass-burning CO2 
     REAL,ALLOCATABLE, SAVE :: flx_co2bb(:,:)     !  biomass-burning CO2 
-!$OMP THREADPRIVATE(flx_co2ff, flx_co2bb)
+!$OMP THREADPRIVATE(flx_co2ff,flx_co2bb)
 
 !! may be controlled via the .def later on
@@ -234 +247 @@
     PRINT *,'probleme avec le mois dans co2_ini =', mth_cur
   ENDIF
-  co2ff(:) = flx_co2ff(:,mth_cur)
-  co2bb(:) = flx_co2bb(:,mth_cur)
+  IF (.NOT.ALLOCATED(fco2_ff)) ALLOCATE(fco2_ff(klon))
+  IF (.NOT.ALLOCATED(fco2_bb)) ALLOCATE(fco2_bb(klon))
+  fco2_ff(:) = flx_co2ff(:,mth_cur)
+  fco2_bb(:) = flx_co2bb(:,mth_cur)
 
   END SUBROUTINE co2_emissions