Changes between Version 4 and Version 5 of Documentation/UserGuide/ImposeVegetation


Ignore:
Timestamp:
2020-02-28T10:37:10+01:00 (4 years ago)
Author:
bastrikov
Comment:

--

Legend:

Unmodified
Added
Removed
Modified

  • Documentation/UserGuide/ImposeVegetation

    v4 v5  
    11= How to prescribe (or impose in ORCHIDEE terminology) the vegetation in a pixel? = 
    22 
    3 Many new ORCHIDEE users will impose their preferred vegetation to a single pixel. This has the advantage that simulations go faster, output files are smaller, data-model comparison is more straightforward and it is easier to understand the effect of specific settings and parameters. To do so you will to make changes to ../config/ORCHIDEE_OL/COMP/sechiba.card and ../config/ORCHIDEE_OL/PARAM/run.def  
     3Author: S.Luyssaert[[BR]] 
     4Last revision: 2020/02/28, V. Bastrikov 
    45 
    5 First, state that you want to impose the vegetation 
     6Objective: Many new ORCHIDEE users will impose their preferred vegetation to a single pixel. This has the advantage that simulations go faster, output files are smaller, data-model comparison is more straightforward and it is easier to understand the effect of specific settings and parameters. 
     7 
     8To prescribe the vegetation you need to make changes to COMP/sechiba.card and PARAM/run.def. 
     9 
     10First, in COMP/sechiba.card set that you don't want to update the land cover map (because you're not using one): 
     11 
    612{{{ 
    7 # IMPOSE_VEG=y/n 
    8 IMPOSE_VEG=y 
    9 }}} 
    10  
    11 Make sure you don't overwrite your own settings in the next time step, hence, switch of land cover change 
    12 {{{ 
    13 # Set LAND_COVER_CHANGE=y if VEGET_UPDATE > 0 
    14 LAND_COVER_CHANGE=n 
    15 }}} 
    16  
    17  
    18 Obviously, you don't want to update the land cover map because you're not using one 
    19 {{{ 
    20 # VEGET_UPDATE=0Y no change in PFTmap. PFTmap should be set only in InitialStateFiles/List.  
    21 # VEGET_UPDATE=1Y change PFTmap every year in december : PFTmap should be added in BoundaryFiles/List and removed from InitialStateFiles/List. 
    22 # This parameter is only used if IMPOSE_VEG=n. 
     13# VEGET_UPDATE=0Y : no change in vegetation map. PFTmap.nc should be set only in InitialStateFiles/List. 
     14# VEGET_UPDATE=1Y : change vegetation map every year. PFTmap.nc should be set only in BoundaryFiles/List. 
     15# Note that the model changed behaviour: 
     16# - since rev 2717 in the trunk, the map is read in January. In BoundaryFiles/List, copy the map for current year, use the variable ${year}. 
     17# - for older version, the map is read in December. In BoundaryFiles/List, copy the map for next coming year using the variable ${year_p1}. 
    2318VEGET_UPDATE=0Y 
    2419}}} 
    2520 
    26 The land cover map should not be specified in Initial state list of the model 
     21The land cover map should not be specified in either "initial state files" or "boundary files" blocks: 
     22 
    2723{{{ 
    2824[InitialStateFiles] 
    29 List=       (${R_INIT}/SRF/${config_UserChoices_TagName}/routing.nc, .),     \ 
    30             (${R_INIT}/SRF/${config_UserChoices_TagName}/soils_param.nc, .), \ 
    31             (${R_INIT}/SRF/${config_UserChoices_TagName}/cartepente2d_15min.nc, .), \ 
    32             (${R_INIT}/SRF/${config_UserChoices_TagName}/floodplains.nc, .) 
     25List=       (${R_IN}/SRF/ROUTING/routing.nc, .),     \ 
     26            (${R_IN}/SRF/SOIL/soils_param4.nc, soils_param.nc), \ 
     27            (${R_IN}/SRF/SOIL/soil_bulk_and_ph.nc, .), \ 
     28            (${R_IN}/SRF/cartepente2d_15min4.nc, cartepente2d_15min.nc), \ 
     29            (${R_IN}/SRF/ROUTING/floodplains4.nc, floodplains.nc), \ 
     30            (${R_IN}/SRF/reftemp4.nc, reftemp.nc), \ 
     31            (${R_IN}/SRF/albedo/alb_bg_modisopt_2D_ESA_v3.nc, alb_bg.nc) 
     32 
     33[BoundaryFiles] 
     34List=   () 
    3335}}} 
    3436 
    35 In this example there are no boundary files 
    36 {{{      
    37 [BoundaryFiles] 
    38 List=   () 
     37Second, in PARAM/run.def state that you want to impose the vegetation: 
     38 
     39{{{ 
     40# Prescribed vegetation (default n) 
     41IMPOSE_VEG = y 
    3942}}} 
    4043 
     44And specify the vegetation distribution you need using the keywords SECHIBA_VEG and SECHIBA_VEGMAX: 
    4145 
    42 Because there is no initial state for the model it needs to be specified in the run.def.  
    43  
    44 For consistency reasons add to the run.def 
    45 {{{ 
    46 impose_veg=y 
    47 }}} 
    48  
    49 It is essential to specify  
    5046{{{ 
    5147SECHIBA_VEG__01 = x1 
    5248... 
    5349SECHIBA_VEG__nvm = xnvm 
     50 
     51SECHIBA_VEGMAX__01 = x1 
     52... 
     53SECHIBA_VEGMAX__nvm = xnvm 
    5454}}} 
    5555 
    56 Where the sum of x1 to xnvm should equal 1 
    57  
    58  
    59 Although some branches no longer make use of it, it is better to specify SECHIBA_VEGMAX as well. Just copy the values of SECHIBA_VEG to avoid problems 
    60  
     56Note that: 
     57 - the number of keywords (from 01 to nvm) should match the number of PFTs used in the ORCHIDEE simulation 
     58 - the sum of x1 to xnvm should equal 1 
     59 - the values set for SECHIBA_VEG and SECHIBA_VEGMAX should be identical.