Changes between Version 7 and Version 8 of Documentation/UserGuide/ImposeVegetation

Timestamp:

2020-02-28T12:43:24+01:00 (4 years ago)

Author:

bastrikov

Comment:

--

Documentation/UserGuide/ImposeVegetation

@@ -1 @@
-= How to prescribe (or impose in ORCHIDEE terminology) the vegetation in a pixel? =
+= How to prescribe the land cover? =
 
 Author: S.Luyssaert[[BR]]
 Last revision: 2020/02/28, V. Bastrikov
 
-Objective: 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.
+Objective: There are several ways to prescribe the vegetation land cover in ORCHIDEE. For a site simulation quite often it is needed to manually impose the preferred vegetation for a single pixel. This has the advantage that simulation 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. For a regional or global simulation it is more appropriate to use currently available vegetation maps. The vegetation then can be read from a land cover map either once at the beginning of a simulation, or changed each year to take into account the land cover change, and possibly the land cover use. This page describes different settings in use for prescribing the land cover with libIGCM.
 
-To prescribe the vegetation you need to make changes to COMP/sechiba.card and PARAM/run.def.
+== Impose land cover for a pixel ==
+
+To manually impose the vegetation you need to make changes to COMP/sechiba.card and PARAM/run.def.
 
 First, in COMP/sechiba.card set that you don't want to update the land cover map (because you're not using one):
@@ -33 +35 @@
 [BoundaryFiles]
 List=   ()
+ListNonDel= ()
 }}}
 
@@ -58 +61 @@
  - the sum of x1 to xnvm should equal 1
  - the values for SECHIBA_VEG and SECHIBA_VEGMAX should be identical
+
+== Read land cover from a map ==
+
+To read the vegetation cover from a map, state that you no longer want to impose the vegetation in PARAM/run.def:
+
+{{{
+# Prescribed vegetation (default n)
+IMPOSE_VEG = n
+}}}
+
+Set that you want to use a single map for the whole simulation (VEGET_UPDATE=0Y) in COMP/sechiba.card:
+
+{{{
+# VEGET_UPDATE=0Y : no change in vegetation map. PFTmap.nc should be set only in InitialStateFiles/List.
+# VEGET_UPDATE=1Y : change vegetation map every year. PFTmap.nc should be set only in BoundaryFiles/List.
+# Note that the model changed behaviour:
+# - 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}.
+# - for older version, the map is read in December. In BoundaryFiles/List, copy the map for next coming year using the variable ${year_p1}.
+VEGET_UPDATE=0Y
+}}}
+
+And add the path to the land cover map in the initial state files list:
+
+{{{
+[InitialStateFiles]
+List=       (${R_IN}/SRF/ROUTING/routing.nc, .),     \
+            (${R_IN}/SRF/SOIL/soils_param4.nc, soils_param.nc), \
+            (${R_IN}/SRF/SOIL/soil_bulk_and_ph.nc, .), \
+            (${R_IN}/SRF/cartepente2d_15min4.nc, cartepente2d_15min.nc), \
+            (${R_IN}/SRF/ROUTING/floodplains4.nc, floodplains.nc), \
+            (${R_IN}/SRF/reftemp4.nc, reftemp.nc), \
+            (${R_IN}/SRF/albedo/alb_bg_modisopt_2D_ESA_v3.nc, alb_bg.nc), \
+            (${R_IN}/SRF/PFTMAPS/CMIP6/ESA-LUH2v2/historical/15PFT.v2/PFTmap_1860.nc, PFTmap.nc)
+
+[BoundaryFiles]
+List=   ()
+ListNonDel= ()
+}}}
+
+Note that:
+ - if you start from existing restart files, the vegetation will not be updated, the settings from restart file will be kept
+ - run.def has priority over the parameters set in sechiba.card, therefore make sure the run.def does NOT contain IMPOSE_VEG = y
+
+== Use land cover change ==
+
+To read different land cover maps at the beginning of each simulation year, set VEGET_UPDATE=1Y in COMP/sechiba.card:
+
+{{{
+# VEGET_UPDATE=0Y : no change in vegetation map. PFTmap.nc should be set only in InitialStateFiles/List.
+# VEGET_UPDATE=1Y : change vegetation map every year. PFTmap.nc should be set only in BoundaryFiles/List.
+# Note that the model changed behaviour:
+# - 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}.
+# - for older version, the map is read in December. In BoundaryFiles/List, copy the map for next coming year using the variable ${year_p1}.
+VEGET_UPDATE=1Y
+}}}
+
+Move the path to the land cover map to the boundary files list, putting also the special field "${year}" in the path to be replaced with the actual year number during each simulation year. Depending on your needs, additional land use maps can be provided (like the wood harvest map in the example below):
+
+{{{
+[InitialStateFiles]
+List=       (${R_IN}/SRF/ROUTING/routing.nc, .),     \
+            (${R_IN}/SRF/SOIL/soils_param4.nc, soils_param.nc), \
+            (${R_IN}/SRF/SOIL/soil_bulk_and_ph.nc, .), \
+            (${R_IN}/SRF/cartepente2d_15min4.nc, cartepente2d_15min.nc), \
+            (${R_IN}/SRF/ROUTING/floodplains4.nc, floodplains.nc), \
+            (${R_IN}/SRF/reftemp4.nc, reftemp.nc), \
+            (${R_IN}/SRF/albedo/alb_bg_modisopt_2D_ESA_v3.nc, alb_bg.nc)
+
+[BoundaryFiles]
+List=   (${R_IN}/SRF/PFTMAPS/CMIP6/ESA-LUH2v2/historical/15PFT.v2/PFTmap_${year}.nc, PFTmap.nc), \
+        (${R_IN}/SRF/WOODHARVEST/LUH2v2/historical4/woodharvest_${year}.nc, woodharvest.nc)
+ListNonDel= ()
+}}}
+
+Note that:
+ - if you set VEGET_UPDATE = nY (n is an INTEGER > 1), the land cover will be updated the first year and after that every n years