Changeset 7326 for branches


Ignore:
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

Location:
branches/ORCHIDEE_2_2/ORCHIDEE
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/diffuco.f90

    r7265 r7326  
    236236     & evap_bare_lim, evap_bare_lim_ns, evapot, evapot_corr, snow, flood_frac, flood_res, frac_nobio, snow_nobio, totfrac_nobio, & 
    237237     & swnet, swdown, coszang, ccanopy, humrel, veget, veget_max, lai, qsintveg, qsintmax, assim_param, & 
    238      & vbeta , vbeta1, vbeta2, vbeta3, vbeta3pot, vbeta4, vbeta5, gsmean, rveget, rstruct, cimean, gpp, co2_to_bm, & 
     238     & vbeta , vbeta1, vbeta2, vbeta3, vbeta3pot, vbeta4, vbeta5, gsmean, rveget, rstruct, cimean, gpp, & 
    239239     & lalo, neighbours, resolution, ptnlev1, precip_rain, frac_age, tot_bare_soil, frac_snow_veg, frac_snow_nobio, & 
    240240     & hist_id, hist2_id) 
    … …  
    288288    REAL(r_std),DIMENSION (kjpindex,nvm), INTENT (in)  :: qsintmax         !! Maximum water on vegetation for interception  
    289289                                                                           !! (kg m^{-2}) 
    290     REAL(r_std),DIMENSION (kjpindex,nvm), INTENT (in)  :: co2_to_bm        !! virtual gpp ((gC m^{-2} dt_sechiba^{-1}), total area)  
    291290    REAL(r_std),DIMENSION (kjpindex,nvm,npco2), INTENT (in) :: assim_param !! min+max+opt temps, vcmax, vjmax 
    292291                                                                           !! for photosynthesis (K ??) 
    … …  
    374373         veget, veget_max, lai, qsintveg, qsintmax, vbeta3, vbeta3pot, & 
    375374         rveget, rstruct, cimean, gsmean, gpp, & 
    376          co2_to_bm, vbeta23, hist_id, indexveg, indexlai, index, kjit, cim) 
     375         vbeta23, hist_id, indexveg, indexlai, index, kjit, cim) 
    377376 
    378377    ! 
    … …  
    13821381                                veget, veget_max, lai, qsintveg, qsintmax, vbeta3, vbeta3pot, rveget, rstruct, & 
    13831382                                cimean, gsmean, gpp, & 
    1384                                 co2_to_bm, vbeta23, hist_id, indexveg, indexlai, index, kjit, cim) 
     1383                                vbeta23, hist_id, indexveg, indexlai, index, kjit, cim) 
    13851384 
    13861385    ! 
    … …  
    14251424    REAL(r_std),DIMENSION (kjpindex,nvm), INTENT (in)        :: vbeta23          !! Beta for fraction of wetted foliage that will  
    14261425                                                                                 !! transpire (unitless)  
    1427     REAL(r_std),DIMENSION (kjpindex,nvm), INTENT (in)        :: co2_to_bm        !! virtual gpp ((gC m^{-2} dt_sechiba ^{-1}), total area) 
    14281426    INTEGER(i_std),INTENT (in)                               :: hist_id          !! _History_ file identifier (-)     
    14291427    INTEGER(i_std),DIMENSION (kjpindex*nvm), INTENT (in) :: indexveg       !! Indeces of the points on the 3D map (-)    
    … …  
    23102308      ! 
    23112309   END DO         ! loop over vegetation types 
    2312    ! 
    2313  
    2314    ! Add virtual gpp (co2_to_bm) to the gpp. 
    2315    ! Virtual gpp can be created when introducing new pft or for correction of carbon fluxes  
    2316    ! for instance for adjustment of Ra at end of the day. 
    2317    gpp(:,:) = gpp(:,:) + co2_to_bm(:,:) 
    23182310       
    23192311   IF (printlev>=3) WRITE (numout,*) ' diffuco_trans_co2 done ' 

  • branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/sechiba.f90

    r7206 r7326  
    234234  REAL(r_std), ALLOCATABLE, SAVE, DIMENSION (:,:)  :: co2_flux       !! CO2 flux (gC/m**2 of average ground/one_day) 
    235235!$OMP THREADPRIVATE(co2_flux) 
    236   REAL(r_std), ALLOCATABLE, SAVE, DIMENSION (:,:)  :: co2_to_bm      !! virtual CO2 flux (gC/m**2 of average ground/s) 
    237 !$OMP THREADPRIVATE(co2_to_bm) 
    238236  REAL(r_std), ALLOCATABLE, SAVE, DIMENSION (:)    :: evapot         !! Soil Potential Evaporation 
    239237!$OMP THREADPRIVATE(evapot) 
    … …  
    459457                              frac_nobio,    njsc,         veget_max,      fraclut,           & 
    460458                              nwdfraclut,    tot_bare_soil,totfrac_nobio,  qsintmax,          & 
    461                               co2_to_bm, temp_growth) 
     459                              temp_growth) 
    462460     
    463461    !! 1.4 Initialize diffusion coefficients 
    … …  
    727725         & frac_nobio, snow_nobio, totfrac_nobio, & 
    728726         & swnet, swdown, coszang, ccanopy, humrel, veget, veget_max, lai, qsintveg, qsintmax, assim_param, & 
    729          & vbeta, vbeta1, vbeta2, vbeta3, vbeta3pot, vbeta4, vbeta5, gsmean, rveget, rstruct, cimean, gpp, co2_to_bm, & 
     727         & vbeta, vbeta1, vbeta2, vbeta3, vbeta3pot, vbeta4, vbeta5, gsmean, rveget, rstruct, cimean, gpp, & 
    730728         & lalo, neighbours, resolution, ptnlev1, precip_rain, frac_age, tot_bare_soil, frac_snow_veg, frac_snow_nobio, & 
    731729         & hist_id, hist2_id) 
    … …  
    809807         lai, frac_age, height, veget, frac_nobio, veget_max, totfrac_nobio, qsintmax, & 
    810808         rest_id, hist_id, hist2_id, rest_id_stom, hist_id_stom, hist_id_stom_IPCC, & 
    811          co2_flux, fco2_lu, fco2_wh, fco2_ha, co2_to_bm, temp_growth, tot_bare_soil) 
     809         co2_flux, fco2_lu, fco2_wh, fco2_ha, temp_growth, tot_bare_soil) 
    812810 
    813811 
    … …  
    14001398                            ks,         nvan,      avan,     mcr,        & 
    14011399                            mcs,        mcfc,      mcw,                  & 
    1402                             co2_to_bm,  assim_param, frac_age) 
     1400                            assim_param, frac_age) 
    14031401     
    14041402    IF (printlev_loc>=3) WRITE (numout,*) 'sechiba_finalize done' 
    … …  
    17381736    IF (ier /= 0) CALL ipslerr_p(3,'sechiba_init','Pb in alloc for co2_flux','','') 
    17391737    co2_flux(:,:)=zero 
    1740  
    1741     ALLOCATE (co2_to_bm(kjpindex,nvm),stat=ier) 
    1742     IF (ier /= 0) CALL ipslerr_p(3,'sechiba_init','Pb in alloc for co2_to_bm','','') 
    17431738 
    17441739    ALLOCATE (shumdiag(kjpindex,nslm),stat=ier) 
    … …  
    19821977    IF ( ALLOCATED (stempdiag)) DEALLOCATE (stempdiag) 
    19831978    IF ( ALLOCATED (co2_flux)) DEALLOCATE (co2_flux) 
    1984     IF ( ALLOCATED (co2_to_bm)) DEALLOCATE (co2_to_bm) 
    19851979    IF ( ALLOCATED (shumdiag)) DEALLOCATE (shumdiag) 
    19861980    IF ( ALLOCATED (shumdiag_perma)) DEALLOCATE (shumdiag_perma) 

  • branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/slowproc.f90

    r7325 r7326  
    283283                                  frac_nobio,    njsc,         veget_max,      fraclut,           & 
    284284                                  nwdfraclut,    tot_bare_soil,totfrac_nobio,  qsintmax,          & 
    285                                   co2_to_bm,     temp_growth) 
     285                                  temp_growth) 
    286286 
    287287!! 0.1 Input variables 
    … …  
    302302     
    303303!! 0.2 Output variables  
    304     REAL(r_std), DIMENSION (kjpindex,nvm), INTENT(out)     :: co2_to_bm      !! Virtual gpp per average ground area (gC m^{-2} dt_stomate^{-1}) 
    305304    REAL(r_std),DIMENSION (kjpindex), INTENT (out)         :: temp_growth    !! Growth temperature (°C) - Is equal to t2m_month  
    306305    INTEGER(i_std), DIMENSION(kjpindex), INTENT(out)       :: njsc           !! Index of the dominant soil textural class in the grid cell (1-nscm, unitless) 
    … …  
    365364             contfrac,       totfrac_nobio,          clayfraction, temp_air,          & 
    366365             lai,            veget,                  veget_max,                       & 
    367              co2_to_bm,      deadleaf_cover,         assim_param,  temp_growth ) 
    368     ELSE 
    369        !! ok_stomate is not activated 
    370        !! Define the CO2 fluxes to zero (no carbone cycle) 
    371        co2_to_bm(:,:) = zero 
     366             deadleaf_cover,         assim_param,  temp_growth ) 
    372367    ENDIF 
    373368     
    … …  
    447442       rest_id, hist_id, hist2_id, rest_id_stom, hist_id_stom, hist_id_stom_IPCC, & 
    448443       co2_flux, fco2_lu, fco2_wh, fco2_ha, & 
    449        co2_to_bm, temp_growth, tot_bare_soil) 
     444       temp_growth, tot_bare_soil) 
    450445   
    451446!! INTERFACE DESCRIPTION 
    … …  
    484479    REAL(r_std),DIMENSION (kjpindex), INTENT (out)      :: fco2_wh             !! CO2 Flux to Atmosphere from Wood Harvesting (gC m^{-2} dt_stomate^{-1}) 
    485480    REAL(r_std),DIMENSION (kjpindex), INTENT (out)      :: fco2_ha             !! CO2 Flux to Atmosphere from Crop Harvesting (gC m^{-2} dt_stomate^{-1}) 
    486     REAL(r_std), DIMENSION (kjpindex,nvm), INTENT(out)  :: co2_to_bm           !! virtual gpp flux per average ground area (gC m^{-2} dt_stomate^{-1}) 
    487481    REAL(r_std),DIMENSION (kjpindex), INTENT (out)      :: temp_growth         !! Growth temperature (°C) - Is equal to t2m_month  
    488482    REAL(r_std), DIMENSION (kjpindex), INTENT(out)      :: tot_bare_soil       !! Total evaporating bare soil fraction in the mesh 
    … …  
    595589            rest_id_stom, hist_id_stom, hist_id_stom_IPCC, & 
    596590            co2_flux, fco2_lu, fco2_wh, fco2_ha, & 
    597             resp_maint, resp_hetero, resp_growth, co2_to_bm, temp_growth) 
     591            resp_maint, resp_hetero, resp_growth, temp_growth) 
    598592 
    599593 
    … …  
    664658       fco2_wh(:) = zero 
    665659       fco2_ha(:) = zero 
    666        co2_to_bm(:,:) = zero 
    667660    ENDIF 
    668661 
    … …  
    746739                                frac_nobio, veget_max, reinf_slope,          & 
    747740                                ks,  nvan, avan, mcr, mcs, mcfc, mcw,        & 
    748                                 co2_to_bm,  assim_param, frac_age ) 
     741                                assim_param, frac_age ) 
    749742 
    750743!! 0.1 Input variables 
    … …  
    768761    REAL(r_std),DIMENSION (kjpindex), INTENT(in)         :: mcfc           !! Volumetric water content at field capacity (m^{3} m^{-3}) 
    769762    REAL(r_std),DIMENSION (kjpindex), INTENT(in)         :: mcw            !! Volumetric water content at wilting point (m^{3} m^{-3}) 
    770  
    771     REAL(r_std),DIMENSION (kjpindex,nvm),INTENT(in)      :: co2_to_bm      !! virtual gpp flux between atmosphere and biosphere 
    772763    REAL(r_std),DIMENSION (kjpindex,nvm,npco2),INTENT (in):: assim_param   !! min+max+opt temperatures & vmax for photosynthesis (K, \mumol m^{-2} s^{-1}) 
    773764    REAL(r_std),DIMENSION (kjpindex,nvm,nleafages), INTENT(in):: frac_age  !! Age efficacity from STOMATE for isoprene 
    … …  
    842833    ! 2.2 Write restart variables managed by STOMATE 
    843834    IF ( ok_stomate ) THEN 
    844        CALL stomate_finalize (kjit, kjpindex, indexLand, clayfraction, co2_to_bm, assim_param)  
     835       CALL stomate_finalize (kjit, kjpindex, indexLand, clayfraction, assim_param)  
    845836    ENDIF 
    846837     

  • branches/ORCHIDEE_2_2/ORCHIDEE/src_stomate/stomate.f90

    r6369 r7326  
    646646         contfrac,       totfrac_nobio,     clay,         temp_air,          & 
    647647         lai,            veget,             veget_max,                       & 
    648          co2_to_bm_radia,deadleaf_cover,    assim_param,  temp_growth ) 
     648         deadleaf_cover,    assim_param,  temp_growth ) 
    649649 
    650650 
    … …  
    674674    !! 0.2 Output variables 
    675675 
    676     REAL(r_std),DIMENSION(kjpindex,nvm),INTENT(out) :: co2_to_bm_radia   !! virtual gpp flux between atmosphere and biosphere 
    677676    REAL(r_std),DIMENSION(kjpindex),INTENT(out)     :: deadleaf_cover    !! Fraction of soil covered by dead leaves (unitless) 
    678677    REAL(r_std),DIMENSION(kjpindex,nvm,npco2),INTENT(out) :: assim_param !! min+max+opt temperatures (K) & vmax for photosynthesis   
    … …  
    805804         leaf_age, leaf_frac, & 
    806805         senescence, when_growthinit, age, & 
    807          resp_hetero_d, resp_maint_d, resp_growth_d, co2_fire, co2_to_bm_dgvm, co2_to_bm_radia, & 
     806         resp_hetero_d, resp_maint_d, resp_growth_d, co2_fire, co2_to_bm_dgvm, & 
    808807         veget_lastlight, everywhere, need_adjacent, RIP_time, & 
    809808         time_hum_min, hum_min_dormance, & 
    … …  
    11961195       &  rest_id_stom, hist_id_stom, hist_id_stom_IPCC, & 
    11971196       &  co2_flux_out, fco2_lu_out, fco2_wh_out, fco2_ha_out, & 
    1198        &  resp_maint, resp_hetero, resp_growth, co2_to_bm_radia, temp_growth) 
     1197       &  resp_maint, resp_hetero, resp_growth, temp_growth) 
    11991198     
    12001199    IMPLICIT NONE 
    … …  
    12621261    REAL(r_std),DIMENSION(kjpindex,nvm),INTENT(out) :: resp_hetero       !! Heterotrophic respiration in   
    12631262                                                                         !! @tex $(gC m^{-2} dt_stomate^{-1})$ @endtex   
    1264     REAL(r_std),DIMENSION(kjpindex,nvm),INTENT(out) :: co2_to_bm_radia   !! Virtual gpp created for equilibrium of carbon mass   
    1265                                                                          !! @tex $(gC m^{-2} dt_stomate^{-1})$ @endtex   
    12661263    REAL(r_std),DIMENSION(kjpindex),INTENT(out)     :: temp_growth       !! Growth temperature (°C)   
    12671264                                                                         !! Is equal to t2m_month  
    … …  
    16021599            &             veget_lastlight, everywhere, need_adjacent, RIP_time, & 
    16031600            &             lai, rprof,npp_daily, turnover_daily, turnover_time,& 
    1604             &             control_moist_inst, control_temp_inst, soilcarbon_input_inst, & 
     1601            &             control_moist_inst, control_temp_inst, soilcarbon_input_daily, & 
    16051602            &             co2_to_bm_dgvm, co2_fire, & 
    16061603            &             resp_hetero_d, resp_hetero_litter_d, resp_hetero_soil_d, resp_maint_d, resp_growth_d, & 
    … …  
    17711768       ! - (6) co2 emission from fire 
    17721769       ! co2_to_bm is not added as it is already encounted in gpp 
    1773        nep_daily(:,:)= gpp_daily(:,:)       & 
     1770       nep_daily(:,:)= gpp_daily(:,:) + co2_to_bm_dgvm(:,:)      & 
    17741771                     - resp_maint_d(:,:)  - resp_growth_d(:,:)   & 
    17751772                     - resp_hetero_d(:,:) - co2_fire(:,:)  
    … …  
    17801777       ! This variable will be used for the coupling to LMDZ for ESM configuration. 
    17811778       co2_flux(:,:) = (resp_hetero_d(:,:) + resp_maint_d(:,:) + resp_growth_d(:,:) & 
    1782             + co2_fire(:,:) - gpp_daily(:,:))*veget_cov_max 
     1779            + co2_fire(:,:) - co2_to_bm_dgvm(:,:) - gpp_daily(:,:))*veget_cov_max 
    17831780       
    17841781       IF ( hist_id_stom_IPCC > 0 ) THEN 
    … …  
    18751872    resp_maint(:,ibare_sechiba) = zero 
    18761873    resp_growth(:,:)= resp_growth_d(:,:)*veget_cov_max(:,:)*dt_sechiba/one_day 
    1877     co2_to_bm_radia(:,:)=co2_to_bm_dgvm(:,:)*veget_cov_max(:,:)*dt_sechiba/one_day 
    18781874    resp_hetero(:,:) = resp_hetero_radia(:,:)*veget_cov_max(:,:) 
    18791875     
    … …  
    20172013!_ ================================================================================================================================ 
    20182014 
    2019   SUBROUTINE stomate_finalize (kjit, kjpindex, index, clay, co2_to_bm_radia, assim_param)  
     2015  SUBROUTINE stomate_finalize (kjit, kjpindex, index, clay, assim_param)  
    20202016     
    20212017    IMPLICIT NONE 
    … …  
    20272023    INTEGER(i_std),DIMENSION(kjpindex),INTENT(in)   :: index             !! Indices of the terrestrial pixels only (unitless) 
    20282024    REAL(r_std),DIMENSION(kjpindex),INTENT(in)      :: clay              !! Clay fraction of soil (0-1, unitless) 
    2029     REAL(r_std),DIMENSION(kjpindex,nvm),INTENT(in)  :: co2_to_bm_radia   !! virtual gpp flux between atmosphere and biosphere 
    20302025    REAL(r_std),DIMENSION(kjpindex,nvm,npco2),INTENT(in) :: assim_param  !! min+max+opt temperatures (K) & vmax for photosynthesis   
    20312026 
    … …  
    21322127         leaf_age, leaf_frac, & 
    21332128         senescence, when_growthinit, age, & 
    2134          resp_hetero_d, resp_maint_d, resp_growth_d, co2_fire, co2_to_bm_dgvm, co2_to_bm_radia, & 
     2129         resp_hetero_d, resp_maint_d, resp_growth_d, co2_fire, co2_to_bm_dgvm, & 
    21352130         veget_lastlight, everywhere, need_adjacent, & 
    21362131         RIP_time, & 

  • branches/ORCHIDEE_2_2/ORCHIDEE/src_stomate/stomate_io.f90

    r6319 r7326  
    7575       &  turnover_longterm, gpp_week, biomass, resp_maint_part, & 
    7676       &  leaf_age, leaf_frac, senescence, when_growthinit, age, & 
    77        &  resp_hetero, resp_maint, resp_growth, co2_fire, co2_to_bm_dgvm, co2_to_bm_radia, & 
     77       &  resp_hetero, resp_maint, resp_growth, co2_fire, co2_to_bm_dgvm, & 
    7878       &  veget_lastlight, everywhere, need_adjacent, RIP_time, & 
    7979       &  time_hum_min, hum_min_dormance, & 
    … …  
    249249    ! biomass uptaken (gC/(m**2 of total ground)/day) 
    250250    REAL(r_std),DIMENSION(npts,nvm),INTENT(out) :: co2_to_bm_dgvm 
    251     ! biomass uptaken (gC/(m**2 of total ground)/dt_sechiba) 
    252     REAL(r_std),DIMENSION(npts,nvm),INTENT(out) :: co2_to_bm_radia 
    253251    ! vegetation fractions (on ground) after last light competition 
    254252    REAL(r_std),DIMENSION(npts,nvm),INTENT(out) :: veget_lastlight 
    … …  
    941939    IF (ALL(co2_to_bm_dgvm(:,:) == val_exp)) co2_to_bm_dgvm(:,:) = zero 
    942940 
    943     co2_to_bm_radia(:,:) = val_exp 
    944     var_name = 'co2_to_bm_radia' 
    945     CALL restget_p (rest_id_stomate, var_name, nbp_glo, nvm     , 1, itime, & 
    946          &                .TRUE., co2_to_bm_radia, 'gather', nbp_glo, index_g) 
    947     IF (ALL(co2_to_bm_radia(:,:) == val_exp)) co2_to_bm_radia(:,:) = zero 
    948941    !- 
    949942    ! 14 vegetation distribution after last light competition 
    … …  
    13471340       &  turnover_longterm, gpp_week, biomass, resp_maint_part, & 
    13481341       &  leaf_age, leaf_frac, senescence, when_growthinit, age, & 
    1349        &  resp_hetero, resp_maint, resp_growth, co2_fire, co2_to_bm_dgvm, co2_to_bm_radia, & 
     1342       &  resp_hetero, resp_maint, resp_growth, co2_fire, co2_to_bm_dgvm, & 
    13501343       &  veget_lastlight, everywhere, need_adjacent, RIP_time, & 
    13511344       &  time_hum_min, hum_min_dormance, & 
    … …  
    15131506    ! biomass uptaken (gC/(m**2 of total ground)/day) 
    15141507    REAL(r_std),DIMENSION(npts,nvm),INTENT(in) :: co2_to_bm_dgvm 
    1515     ! biomass uptaken (gC/(m**2 of total ground)/dt_sechiba) 
    1516     REAL(r_std),DIMENSION(npts,nvm),INTENT(in) :: co2_to_bm_radia 
    15171508    ! vegetation fractions (on ground) after last light competition 
    15181509    REAL(r_std),DIMENSION(npts,nvm),INTENT(in) :: veget_lastlight 
    … …  
    20242015    CALL restput_p (rest_id_stomate, var_name, nbp_glo, nvm, 1, itime, & 
    20252016         &                co2_to_bm_dgvm, 'scatter', nbp_glo, index_g) 
    2026     !- 
    2027     var_name = 'co2_to_bm_radia' 
    2028     CALL restput_p (rest_id_stomate, var_name, nbp_glo, nvm, 1, itime, & 
    2029          &                co2_to_bm_radia, 'scatter', nbp_glo, index_g) 
    20302017    !- 
    20312018    ! 14 vegetation distribution after last light competition 

  • branches/ORCHIDEE_2_2/ORCHIDEE/src_stomate/stomate_lpj.f90

    r6369 r7326  
    875875    sum_cVegTree = zero 
    876876 
    877     DO j=2,nvm 
     877    DO j=1,nvm 
    878878 
    879879       tot_litter_carb(:,j) = tot_litter_carb(:,j) + (litter(:,istructural,j,iabove,icarbon) + & 
    … …  
    958958    CALL xios_orchidee_send_field("LAI",lai) 
    959959    CALL xios_orchidee_send_field("VEGET_COV_MAX",veget_cov_max) 
    960     CALL xios_orchidee_send_field("NPP_STOMATE",npp_daily) 
    961     CALL xios_orchidee_send_field("GPP",gpp_daily) 
     960    CALL xios_orchidee_send_field("NPP_STOMATE",npp_daily+co2_to_bm) 
     961    CALL xios_orchidee_send_field("GPP",gpp_daily+co2_to_bm) 
    962962    CALL xios_orchidee_send_field("IND",ind) 
    963963    CALL xios_orchidee_send_field("CN_IND",cn_ind) 
    … …  
    10531053     
    10541054    ! Carbon fluxes transformed from gC/m2/d into kgC/m2/s 
    1055     CALL xios_orchidee_send_field("gpp_ipcc",SUM(gpp_daily*veget_cov_max,dim=2)/1e3/one_day) 
     1055    CALL xios_orchidee_send_field("gpp_ipcc",SUM((gpp_daily+co2_to_bm)*veget_cov_max,dim=2)/1e3/one_day) 
    10561056    CALL xios_orchidee_send_field("ra",SUM((resp_maint+resp_growth)*veget_cov_max,dim=2)/1e3/one_day) 
    10571057    vartmp(:)=zero 
    … …  
    10791079    CALL xios_orchidee_send_field("raTree",vartmp/1e3/one_day) 
    10801080 
    1081     CALL xios_orchidee_send_field("npp_ipcc",SUM(npp_daily*veget_cov_max,dim=2)/1e3/one_day) 
     1081    CALL xios_orchidee_send_field("npp_ipcc",SUM((npp_daily+co2_to_bm)*veget_cov_max,dim=2)/1e3/one_day) 
     1082 
    10821083    vartmp(:)=zero 
    10831084    DO j = 2, nvm 
    10841085       IF ( .NOT. is_tree(j) .AND. natural(j) ) THEN 
    1085           vartmp(:) = vartmp(:) + npp_daily(:,j)*veget_cov_max(:,j) 
     1086          vartmp(:) = vartmp(:) +( npp_daily(:,j)+co2_to_bm(:,j))*veget_cov_max(:,j) 
    10861087       ENDIF 
    10871088    ENDDO 
    … …  
    10901091    DO j = 2, nvm 
    10911092       IF ( (.NOT. is_tree(j)) .AND. (.NOT. natural(j)) ) THEN 
    1092           vartmp(:) = vartmp(:) + npp_daily(:,j)*veget_cov_max(:,j) 
     1093          vartmp(:) = vartmp(:) + (npp_daily(:,j)+co2_to_bm(:,j))*veget_cov_max(:,j) 
    10931094       ENDIF 
    10941095    ENDDO 
    … …  
    10981099    DO j = 2, nvm 
    10991100       IF ( is_tree(j) ) THEN 
    1100           vartmp(:) = vartmp(:) + npp_daily(:,j)*veget_cov_max(:,j) 
     1101          vartmp(:) = vartmp(:) + (npp_daily(:,j)+co2_to_bm(:,j))*veget_cov_max(:,j) 
    11011102       ENDIF 
    11021103    ENDDO 
    … …  
    11961197    CALL xios_orchidee_send_field("flulccatmlut",flulccatmlut) 
    11971198 
    1198    ! co2_to_bm is not added as it is already included in gpp 
    1199     CALL xios_orchidee_send_field("nbp",(SUM((gpp_daily-(resp_maint+resp_growth+resp_hetero)-co2_fire) * & 
     1199    CALL xios_orchidee_send_field("nbp",(SUM((gpp_daily+co2_to_bm-(resp_maint+resp_growth+resp_hetero)-co2_fire) * & 
    12001200          veget_cov_max,dim=2)-cflux_prod_total-cflux_prod_harvest_total-harvest_above)/1e3/one_day) 
    12011201    CALL xios_orchidee_send_field("fVegLitter",SUM((tot_bm_to_litter(:,:,icarbon) + tot_turnover(:,:,icarbon))*& 
    … …  
    12671267               resp_hetero(:,j)*veget_cov_max(:,j)/1e3/one_day 
    12681268          npplut(:,id_psl) = npplut(:,id_psl) + & 
    1269                npp_daily(:,j)*veget_cov_max(:,j)/1e3/one_day 
     1269               (npp_daily(:,j)+co2_to_bm(:,j))*veget_cov_max(:,j)/1e3/one_day 
    12701270       ELSE 
    12711271          clitterlut(:,id_crp) = clitterlut(:,id_crp) + tot_litter_carb(:,j)*veget_cov_max(:,j)/1e3 
    … …  
    12781278               resp_hetero(:,j)*veget_cov_max(:,j)/1e3/one_day 
    12791279          npplut(:,id_crp) = npplut(:,id_crp) + & 
    1280                npp_daily(:,j)*veget_cov_max(:,j)/1e3/one_day 
     1280               (npp_daily(:,j)+co2_to_bm(:,j))*veget_cov_max(:,j)/1e3/one_day 
    12811281       END IF 
    12821282    END DO 
    … …  
    14221422         veget_cov_max, npts*nvm, horipft_index) 
    14231423    CALL histwrite_p (hist_id_stomate, 'NPP', itime, & 
    1424          npp_daily, npts*nvm, horipft_index) 
     1424         npp_daily+co2_to_bm, npts*nvm, horipft_index) 
    14251425    CALL histwrite_p (hist_id_stomate, 'GPP', itime, & 
    1426          gpp_daily, npts*nvm, horipft_index) 
     1426         gpp_daily+co2_to_bm, npts*nvm, horipft_index) 
    14271427    CALL histwrite_p (hist_id_stomate, 'IND', itime, & 
    14281428         ind, npts*nvm, horipft_index) 
    … …  
    15281528       CALL histwrite_p (hist_id_stomate_IPCC, "lai", itime, & 
    15291529            vartmp, npts, hori_index) 
    1530        vartmp(:)=SUM(gpp_daily*veget_cov_max,dim=2)/1e3/one_day 
     1530       vartmp(:)=SUM((gpp_daily+co2_to_bm)*veget_cov_max,dim=2)/1e3/one_day 
    15311531       CALL histwrite_p (hist_id_stomate_IPCC, "gpp", itime, & 
    15321532            vartmp, npts, hori_index) 
    … …  
    15341534       CALL histwrite_p (hist_id_stomate_IPCC, "ra", itime, & 
    15351535            vartmp, npts, hori_index) 
    1536        vartmp(:)=SUM(npp_daily*veget_cov_max,dim=2)/1e3/one_day 
     1536       vartmp(:)=SUM((npp_daily+co2_to_bm)*veget_cov_max,dim=2)/1e3/one_day 
    15371537       CALL histwrite_p (hist_id_stomate_IPCC, "npp", itime, & 
    15381538            vartmp, npts, hori_index) 
    … …  
    15531553            vartmp, npts, hori_index) 
    15541554       ! co2_to_bm is not added as it is already included in gpp 
    1555        vartmp(:)=(SUM((gpp_daily-(resp_maint+resp_growth+resp_hetero)-co2_fire) & 
     1555       vartmp(:)=(SUM((gpp_daily+co2_to_bm-(resp_maint+resp_growth+resp_hetero)-co2_fire) & 
    15561556            &        *veget_cov_max,dim=2)-cflux_prod_total-cflux_prod_harvest_total-harvest_above)/1e3/one_day 
    15571557       CALL histwrite_p (hist_id_stomate_IPCC, "nbp", itime, & 
Note: See TracChangeset for help on using the changeset viewer.