Changeset 3857 for LMDZ6/trunk/libf/phylmd/carbon_cycle_mod.F90

Timestamp:

Feb 26, 2021, 8:37:59 PM (3 years ago)

Author:

oboucher

Message:

Additions for the interactive carbon cycle

File:

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

LMDZ6/trunk/libf/phylmd/carbon_cycle_mod.F90

@@ -39 +39 @@
   LOGICAL, PUBLIC :: carbon_cycle_tr        ! 3D transport of CO2 in the atmosphere, parameter read in conf_phys
 !$OMP THREADPRIVATE(carbon_cycle_tr)
-  LOGICAL, PUBLIC :: carbon_cycle_rad       ! CO2 interactive radiatively 
+  LOGICAL, PUBLIC :: carbon_cycle_rad       ! flag to activate CO2 interactive radiatively 
 !$OMP THREADPRIVATE(carbon_cycle_rad)
-  INTEGER, PUBLIC :: level_coupling_esm ! Level of coupling for the ESM - 0, 1, 2, 3
+  INTEGER, PUBLIC :: level_coupling_esm     ! Level of coupling for the ESM - 0, 1, 2, 3
 !$OMP THREADPRIVATE(level_coupling_esm)
+  LOGICAL, PUBLIC :: read_fco2_ocean_cor    ! flag to read corrective oceanic CO2 flux
+!$OMP THREADPRIVATE(read_fco2_ocean_cor)  
+  REAL, PUBLIC :: var_fco2_ocean_cor        ! corrective oceanic CO2 flux
+!$OMP THREADPRIVATE(var_fco2_ocean_cor)
+  REAL, PUBLIC :: ocean_area_tot            ! total oceanic area to convert flux
+!$OMP THREADPRIVATE(ocean_area_tot)
+  LOGICAL, PUBLIC :: read_fco2_land_cor     ! flag to read corrective land CO2 flux
+!$OMP THREADPRIVATE(read_fco2_land_cor)  
+  REAL, PUBLIC :: var_fco2_land_cor         ! corrective land CO2 flux
+!$OMP THREADPRIVATE(var_fco2_land_cor)
+  REAL, PUBLIC :: land_area_tot             ! total land area to convert flux
+!$OMP THREADPRIVATE(land_area_tot)
+
   REAL, PUBLIC :: RCO2_glo
 !$OMP THREADPRIVATE(RCO2_glo)
@@ -95 +108 @@
   REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_ocean ! Net flux from ocean [kgCO2/m2/s]
 !$OMP THREADPRIVATE(fco2_ocean)
+  REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_ocean_cor ! Net corrective flux from ocean [kgCO2/m2/s]
+!$OMP THREADPRIVATE(fco2_ocean_cor)
+  REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_land_cor  ! Net corrective flux from land [kgCO2/m2/s]
+!$OMP THREADPRIVATE(fco2_land_cor)
 
   REAL, DIMENSION(:,:), ALLOCATABLE :: dtr_add       ! Tracer concentration to be injected 
@@ -252 +269 @@
        IF (.NOT.ALLOCATED(fco2_ocean)) ALLOCATE(fco2_ocean(klon), stat=ierr)
        IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation fco2_ocean',1)
-       fco2_bb(1:klon) = 0.
+       fco2_ocean(1:klon) = 0.
+
+       IF (.NOT.ALLOCATED(fco2_ocean_cor)) ALLOCATE(fco2_ocean_cor(klon), stat=ierr)
+       IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation fco2_ocean_cor',1)
+       fco2_ocean_cor(1:klon) = 0.
+       IF (.NOT.ALLOCATED(fco2_land_cor)) ALLOCATE(fco2_land_cor(klon), stat=ierr)
+       IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation fco2_land_cor',1)
+       fco2_land_cor(1:klon) = 0.
+
     ENDIF