Changeset 7326 for branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/slowproc.f90

Timestamp:

2021-10-20T18:39:22+02:00 (3 years ago)

Author:

josefine.ghattas

Message:

Corrected bug on carbon balance closure. See ticket #785
Integration in branch 2_2 done by P. Cadule

File:

• branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/slowproc.f90 (modified) (10 diffs)

branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/slowproc.f90

@@ -283 +283 @@
                                   frac_nobio,    njsc,         veget_max,      fraclut,           &
                                   nwdfraclut,    tot_bare_soil,totfrac_nobio,  qsintmax,          &
-                                  co2_to_bm,     temp_growth)
+                                  temp_growth)
 
 !! 0.1 Input variables
@@ -302 +302 @@
     
 !! 0.2 Output variables 
-    REAL(r_std), DIMENSION (kjpindex,nvm), INTENT(out)     :: co2_to_bm      !! Virtual gpp per average ground area (gC m^{-2} dt_stomate^{-1})
     REAL(r_std),DIMENSION (kjpindex), INTENT (out)         :: temp_growth    !! Growth temperature (°C) - Is equal to t2m_month 
     INTEGER(i_std), DIMENSION(kjpindex), INTENT(out)       :: njsc           !! Index of the dominant soil textural class in the grid cell (1-nscm, unitless)
@@ -365 +364 @@
              contfrac,       totfrac_nobio,          clayfraction, temp_air,          &
              lai,            veget,                  veget_max,                       &
-             co2_to_bm,      deadleaf_cover,         assim_param,  temp_growth )
-    ELSE
-       !! ok_stomate is not activated
-       !! Define the CO2 fluxes to zero (no carbone cycle)
-       co2_to_bm(:,:) = zero
+             deadleaf_cover,         assim_param,  temp_growth )
     ENDIF
     
@@ -447 +442 @@
        rest_id, hist_id, hist2_id, rest_id_stom, hist_id_stom, hist_id_stom_IPCC, &
        co2_flux, fco2_lu, fco2_wh, fco2_ha, &
-       co2_to_bm, temp_growth, tot_bare_soil)
+       temp_growth, tot_bare_soil)
   
 !! INTERFACE DESCRIPTION
@@ -484 +479 @@
     REAL(r_std),DIMENSION (kjpindex), INTENT (out)      :: fco2_wh             !! CO2 Flux to Atmosphere from Wood Harvesting (gC m^{-2} dt_stomate^{-1})
     REAL(r_std),DIMENSION (kjpindex), INTENT (out)      :: fco2_ha             !! CO2 Flux to Atmosphere from Crop Harvesting (gC m^{-2} dt_stomate^{-1})
-    REAL(r_std), DIMENSION (kjpindex,nvm), INTENT(out)  :: co2_to_bm           !! virtual gpp flux per average ground area (gC m^{-2} dt_stomate^{-1})
     REAL(r_std),DIMENSION (kjpindex), INTENT (out)      :: temp_growth         !! Growth temperature (°C) - Is equal to t2m_month 
     REAL(r_std), DIMENSION (kjpindex), INTENT(out)      :: tot_bare_soil       !! Total evaporating bare soil fraction in the mesh
@@ -595 +589 @@
             rest_id_stom, hist_id_stom, hist_id_stom_IPCC, &
             co2_flux, fco2_lu, fco2_wh, fco2_ha, &
-            resp_maint, resp_hetero, resp_growth, co2_to_bm, temp_growth)
+            resp_maint, resp_hetero, resp_growth, temp_growth)
 
 
@@ -664 +658 @@
        fco2_wh(:) = zero
        fco2_ha(:) = zero
-       co2_to_bm(:,:) = zero
     ENDIF
 
@@ -746 +739 @@
                                 frac_nobio, veget_max, reinf_slope,          &
                                 ks,  nvan, avan, mcr, mcs, mcfc, mcw,        &
-                                co2_to_bm,  assim_param, frac_age )
+                                assim_param, frac_age )
 
 !! 0.1 Input variables
@@ -768 +761 @@
     REAL(r_std),DIMENSION (kjpindex), INTENT(in)         :: mcfc           !! Volumetric water content at field capacity (m^{3} m^{-3})
     REAL(r_std),DIMENSION (kjpindex), INTENT(in)         :: mcw            !! Volumetric water content at wilting point (m^{3} m^{-3})
-
-    REAL(r_std),DIMENSION (kjpindex,nvm),INTENT(in)      :: co2_to_bm      !! virtual gpp flux between atmosphere and biosphere
     REAL(r_std),DIMENSION (kjpindex,nvm,npco2),INTENT (in):: assim_param   !! min+max+opt temperatures & vmax for photosynthesis (K, \mumol m^{-2} s^{-1})
     REAL(r_std),DIMENSION (kjpindex,nvm,nleafages), INTENT(in):: frac_age  !! Age efficacity from STOMATE for isoprene
@@ -842 +833 @@
     ! 2.2 Write restart variables managed by STOMATE
     IF ( ok_stomate ) THEN
-       CALL stomate_finalize (kjit, kjpindex, indexLand, clayfraction, co2_to_bm, assim_param) 
+       CALL stomate_finalize (kjit, kjpindex, indexLand, clayfraction, assim_param) 
     ENDIF