Ignore:
Timestamp:
2019-11-08T13:00:52+01:00 (5 years ago)
Author:
josefine.ghattas
Message:

Improvment for the ESM CO2 configuration:

  • Separate variable fco2_lu into 3 parts: fco2_lu, fco2_wh and fco2_ha
  • Move calculation of co2_flux from dt_sechiba time-step to daily time-step (in the part for do_slow)
  • Removed co2_flux and fco2_lu from stomate_intialize argument list. These variables were never used in the intialization phase.
  • Add co2_flux, and fco2_wh, fco2_ha to restart file
  • Corrected output unit for nee to be consistent with LMDZ and stomate output variables. It is now in kgC/m2/s.
  • Corrected output for znetco2
  • Added fCO2_fWoodharvest and fCO2_fHarvest as new possible tracers in LMDZ (intersurf).
  • Added diagnostic output for fCO2_fWoodharvest and fCO2_fHarvest
  1. Cadule
File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r6145 r6319  
    248248                                  frac_nobio,    njsc,         veget_max,      fraclut,           & 
    249249                                  nwdfraclut,    tot_bare_soil,totfrac_nobio,  qsintmax,          & 
    250                                   co2_flux,      co2_to_bm,    fco2_lu,      temp_growth) 
     250                                  co2_to_bm,     temp_growth) 
    251251 
    252252!! 0.1 Input variables 
    … …  
    267267     
    268268!! 0.2 Output variables  
    269     REAL(r_std), DIMENSION (kjpindex,nvm), INTENT(out)     :: co2_flux       !! CO2 flux per average ground area (gC m^{-2} dt_stomate^{-1}) 
    270269    REAL(r_std), DIMENSION (kjpindex,nvm), INTENT(out)     :: co2_to_bm      !! Virtual gpp per average ground area (gC m^{-2} dt_stomate^{-1}) 
    271     REAL(r_std),DIMENSION (kjpindex), INTENT (out)         :: fco2_lu        !! CO2 flux from land-use (without forest management) (gC m^{-2} dt_stomate^{-1}) 
    272270    REAL(r_std),DIMENSION (kjpindex), INTENT (out)         :: temp_growth    !! Growth temperature (°C) - Is equal to t2m_month  
    273271    INTEGER(i_std), DIMENSION(kjpindex), INTENT(out)       :: njsc           !! Index of the dominant soil textural class in the grid cell (1-nscm, unitless) 
    … …  
    322320             contfrac,       totfrac_nobio,          clayfraction, temp_air,          & 
    323321             lai,            veget,                  veget_max,                       & 
    324              co2_flux,       co2_to_bm,              fco2_lu,                         & 
    325              deadleaf_cover, assim_param,            temp_growth ) 
     322             co2_to_bm,      deadleaf_cover,         assim_param,  temp_growth ) 
    326323    ELSE 
    327324       !! ok_stomate is not activated 
    328325       !! Define the CO2 fluxes to zero (no carbone cycle) 
    329        co2_flux(:,:) = zero 
    330326       co2_to_bm(:,:) = zero 
    331327    ENDIF 
    … …  
    383379!! RECENT CHANGE(S): None 
    384380!! 
    385 !! MAIN OUTPUT VARIABLE(S):  ::co2_flux, ::fco2_lu, ::lai, ::height, ::veget, ::frac_nobio,   
     381!! MAIN OUTPUT VARIABLE(S):  ::co2_flux, ::fco2_lu,::fco2_wh, ::fco2_ha, ::lai, ::height, ::veget, ::frac_nobio,   
    386382!! ::veget_max, ::woodharvest, ::totfrac_nobio, ::soiltype, ::assim_param, ::deadleaf_cover, ::qsintmax, 
    387383!! and resp_maint, resp_hetero, resp_growth, npp that are calculated and stored 
    … …  
    405401       lai, frac_age, height, veget, frac_nobio, veget_max, totfrac_nobio, qsintmax, & 
    406402       rest_id, hist_id, hist2_id, rest_id_stom, hist_id_stom, hist_id_stom_IPCC, & 
    407        co2_flux, fco2_lu, co2_to_bm, temp_growth, tot_bare_soil) 
     403       co2_flux, fco2_lu, fco2_wh, fco2_ha, & 
     404       co2_to_bm, temp_growth, tot_bare_soil) 
    408405   
    409406!! INTERFACE DESCRIPTION 
    … …  
    439436!! 0.2 Output variables  
    440437    REAL(r_std), DIMENSION (kjpindex,nvm), INTENT(out)  :: co2_flux            !! CO2 flux per average ground area (gC m^{-2} dt_stomate^{-1}) 
     438    REAL(r_std),DIMENSION (kjpindex), INTENT (out)      :: fco2_lu             !! CO2 flux from land-use (without forest management) (gC m^{-2} dt_stomate^{-1}) 
     439    REAL(r_std),DIMENSION (kjpindex), INTENT (out)      :: fco2_wh             !! CO2 Flux to Atmosphere from Wood Harvesting (gC m^{-2} dt_stomate^{-1}) 
     440    REAL(r_std),DIMENSION (kjpindex), INTENT (out)      :: fco2_ha             !! CO2 Flux to Atmosphere from Crop Harvesting (gC m^{-2} dt_stomate^{-1}) 
    441441    REAL(r_std), DIMENSION (kjpindex,nvm), INTENT(out)  :: co2_to_bm           !! virtual gpp flux per average ground area (gC m^{-2} dt_stomate^{-1}) 
    442     REAL(r_std),DIMENSION (kjpindex), INTENT (out)      :: fco2_lu             !! CO2 flux from land-use (without forest management) (gC m^{-2} dt_stomate^{-1}) 
    443442    REAL(r_std),DIMENSION (kjpindex), INTENT (out)      :: temp_growth         !! Growth temperature (°C) - Is equal to t2m_month  
    444443    REAL(r_std), DIMENSION (kjpindex), INTENT(out)      :: tot_bare_soil       !! Total evaporating bare soil fraction in the mesh 
    … …  
    550549            veget_max_new, woodharvest, totfrac_nobio_new, fraclut, & 
    551550            rest_id_stom, hist_id_stom, hist_id_stom_IPCC, & 
    552             co2_flux, fco2_lu, resp_maint,resp_hetero,resp_growth,co2_to_bm,temp_growth) 
     551            co2_flux, fco2_lu, fco2_wh, fco2_ha, & 
     552            resp_maint, resp_hetero, resp_growth, co2_to_bm, temp_growth) 
     553 
    553554 
    554555       !! 4.2 Output the respiration terms and the net primary 
    … …  
    615616       !! Define the CO2 flux from the grid point to zero (no carbone cycle) 
    616617       co2_flux(:,:) = zero 
     618       fco2_lu(:) = zero 
     619       fco2_wh(:) = zero 
     620       fco2_ha(:) = zero 
    617621       co2_to_bm(:,:) = zero 
    618622    ENDIF 
Note: See TracChangeset for help on using the changeset viewer.