Changeset 7326 for branches/ORCHIDEE_2_2/ORCHIDEE/src_stomate/stomate_io.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_stomate/stomate_io.f90 (modified) (6 diffs)

branches/ORCHIDEE_2_2/ORCHIDEE/src_stomate/stomate_io.f90

@@ -75 +75 @@
        &  turnover_longterm, gpp_week, biomass, resp_maint_part, &
        &  leaf_age, leaf_frac, senescence, when_growthinit, age, &
-       &  resp_hetero, resp_maint, resp_growth, co2_fire, co2_to_bm_dgvm, co2_to_bm_radia, &
+       &  resp_hetero, resp_maint, resp_growth, co2_fire, co2_to_bm_dgvm, &
        &  veget_lastlight, everywhere, need_adjacent, RIP_time, &
        &  time_hum_min, hum_min_dormance, &
@@ -249 +249 @@
     ! biomass uptaken (gC/(m**2 of total ground)/day)
     REAL(r_std),DIMENSION(npts,nvm),INTENT(out) :: co2_to_bm_dgvm
-    ! biomass uptaken (gC/(m**2 of total ground)/dt_sechiba)
-    REAL(r_std),DIMENSION(npts,nvm),INTENT(out) :: co2_to_bm_radia
     ! vegetation fractions (on ground) after last light competition
     REAL(r_std),DIMENSION(npts,nvm),INTENT(out) :: veget_lastlight
@@ -941 +939 @@
     IF (ALL(co2_to_bm_dgvm(:,:) == val_exp)) co2_to_bm_dgvm(:,:) = zero
 
-    co2_to_bm_radia(:,:) = val_exp
-    var_name = 'co2_to_bm_radia'
-    CALL restget_p (rest_id_stomate, var_name, nbp_glo, nvm     , 1, itime, &
-         &                .TRUE., co2_to_bm_radia, 'gather', nbp_glo, index_g)
-    IF (ALL(co2_to_bm_radia(:,:) == val_exp)) co2_to_bm_radia(:,:) = zero
     !-
     ! 14 vegetation distribution after last light competition
@@ -1347 +1340 @@
        &  turnover_longterm, gpp_week, biomass, resp_maint_part, &
        &  leaf_age, leaf_frac, senescence, when_growthinit, age, &
-       &  resp_hetero, resp_maint, resp_growth, co2_fire, co2_to_bm_dgvm, co2_to_bm_radia, &
+       &  resp_hetero, resp_maint, resp_growth, co2_fire, co2_to_bm_dgvm, &
        &  veget_lastlight, everywhere, need_adjacent, RIP_time, &
        &  time_hum_min, hum_min_dormance, &
@@ -1513 +1506 @@
     ! biomass uptaken (gC/(m**2 of total ground)/day)
     REAL(r_std),DIMENSION(npts,nvm),INTENT(in) :: co2_to_bm_dgvm
-    ! biomass uptaken (gC/(m**2 of total ground)/dt_sechiba)
-    REAL(r_std),DIMENSION(npts,nvm),INTENT(in) :: co2_to_bm_radia
     ! vegetation fractions (on ground) after last light competition
     REAL(r_std),DIMENSION(npts,nvm),INTENT(in) :: veget_lastlight
@@ -2024 +2015 @@
     CALL restput_p (rest_id_stomate, var_name, nbp_glo, nvm, 1, itime, &
          &                co2_to_bm_dgvm, 'scatter', nbp_glo, index_g)
-    !-
-    var_name = 'co2_to_bm_radia'
-    CALL restput_p (rest_id_stomate, var_name, nbp_glo, nvm, 1, itime, &
-         &                co2_to_bm_radia, 'scatter', nbp_glo, index_g)
     !-
     ! 14 vegetation distribution after last light competition