Add wood harvest-related Carbon fluxes

[nvuilsce]

Based on what has been developed by Nicolas Viovy, we have included the process of wood harvest and its consequences on Carbon fluxes.
It uses an additional forcing which corresponds to the harvest wood carbon mass per unit area of forest in each grid cell (gC m-2[forest]).

It is developed in a separate routine named woodharvest_main in the module stomate_woodharvest.

It is flaged by the parameter DO_HARVEST_WOOD.
For the moment, it uses the parameter used for land-use VEGET_UPDATE, to set the frequency at which the forcing file for woodharvest is updated.

By default, the name of the forcing file for woodharvest is "woodharvest.nc" and the default name of the variable read in this file is "woodharvest". They can be modified by the respective keywords, WOODHARVEST_FILE and WOODHARVEST_VARIABLE.

The carbon content in the harvested wood will have the same treatment than the one from land-use change (deforestation). Consequently, the additional output variables are:

• "PROD10_HARVEST" : 10 year wood product pool
• "FLUX10_HARVEST" : annual flux for each 10 year wood product pool
• "PROD100_HARVEST" : 100 year wood product pool
• "FLUX100_HARVEST" : annual flux for each 100 year wood product pool
• "CONVFLUX_HARVEST" : annual release right after wood harvest
• "CFLUX_PROD10_HARVEST" : annual release from all 10 year wood product pools
• "CFLUX_PROD100_HARVEST" : annual release from all 100 year wood product pools
• "WOOD_HARVEST" : harvested wood biomass as read in the map
• "WOOD_HARVEST_PFT" : harvested wood biomass used per PFT

[nvuilsce]

This new feature has been comitted in the trunk in revision r4657. Three other minor commits are also related to this development (r4658, r4659, r4660)

[nvuilsce]

script for generating the wood harvest data

[nvuilsce]

Forcing data used for defining the wood harvest has been produced. This data is derived from LUH2v2h. It is computed as the total wood biomass (primf_bioh+primn_bioh+secmf_bioh+secyf_bioh+secnf_bioh) divided by the related land cell area ((primf+primn+secdf+secdn)*area).
The script used for generating this data is in attachment.

[jgipsl]

Changed parmetername DO_HARVEST_WOOD into DO_WOOD_HARVEST to be coherent with module name. Done in trunk [4677].

Script for adapting forcing files are added to SVN here : ​http://forge.ipsl.jussieu.fr/igcmg/browser/TOOLS/CMIP6_FORCING/WOOD_HARVEST/generate_woodharvest.jnl

Forcing files are added to shared repository here .../IGCM/SRF/WOODHARVEST/LUH2v2/historical

[jgipsl]

Bug correction done in the trunk [4847]