| 33 | |
| 34 | The are still some issues with the multi-layer energy budget, and it is currently only possible to run for one PFT. Thus, we suggest you use the single layer energy budget. This can, however, be done by two different methods that gives the exact same results: |
| 35 | |
| 36 | A: use the energy scheme as found in the original enerbil.f90 |
| 37 | |
| 38 | B: use the multi-layer energy scheme for a single canopy layer. This collapses the multi-layer energy scheme to the original enerbil code, but if you wish to take water stress into account in your simulation, you need to choose this method. |
| 39 | |
| 40 | Several flags exist to control these choices. For simplicity, a flag has been made to automatically control several of the flags related to the multi-layering, called '''multi_layer_control'''. Multi_layer_control has 4 possibilities |
| 41 | |
| 42 | 1 – single layer, using the multi layer energy scheme (i.e. choice B above) |
| 43 | |
| 44 | 2 – multi-layer energy budget |
| 45 | |
| 46 | 3 – user specific as set in the run.def |
| 47 | |
| 48 | 4 – single layer, using the original enerbil scheme (i.e. choice A above) |
| 49 | |
| 50 | The default setting of multi_layer_control in ORHIDEE-CN-CAN is 1. This is well tested and save to use. More details of the flags controlling the multi-layer can be found below and within the model code. |
| 51 | |
| 52 | The flags controlled by the multi_layer_control are: |
| 53 | |
| 54 | '''ok_hydrol_arch:''' |
| 55 | Flag that activates the hydraulic architecture routine (true/false). The trunk version of ORCHIDEE (false) uses soil water as a proxy for water stress and applies the stress to Vcmax. When set to true the hydraulic architecture of the vegetation is accounted for to calculate the amount of water that can be transported through the plant given the soil and leaf potential and the conductivities of the roots, wood and leaves. Water supply through the plant is compared against the atmospheric demand for water. If the supply is smaller then the demand, the plant experiences water stress and the stomata will be closed (water stress is now on gs rather than Vcmax). Note that whether stomatal regulation is used or not is controled by a separate flag: ok_gs_feedback. |
| 56 | ''' |
| 57 | ok_gs_feedback:''' |
| 58 | Flag that activates water stress on stomata (true/false). This flag is for debugging only! It allows developers to calculate GPP without any water stress. If the model is used in production mode and ok_hydrol_arch is true this flag should be true as well. |
| 59 | |
| 60 | '''ok_mleb:''' |
| 61 | Flag that activates the multilayer energy budget (true/false). The model uses 10 (default) canopy layers to calculate the albedo, transmittance, absorbance and GPP. These canopy layers can be combined with 10 (default) layers below and 9 layers above the canopy to calculate the energy budget (ok_mleb=y). If set to no, this flag will make the model use 10 layers for the canopy albedo, transmittance, absorbance and GPP and just a single layer for the energy budget. Be aware that if you wish to run with hydraulic architechture ok_mleb needs to be se to true as well. Furthermore, if you wish to run with the original energy scheme (enerbil), set the layers for mleb to 1. |
| 62 | |
| 63 | '''ok_impose_can_structure:''' |
| 64 | This flag is for debugging only! It allows developers to use a prescribed canopy structure rather then the structure calculate by ORCHIDEE. The flag activates the sections of code which directly link the energy budget scheme to the the size and LAI profile of the canopy for the respective PFT and age class that is calculated in stomate, for the albedo. If set to TRUE and the multi-layer budget is activated the model takes LAI profile information and canopy level heights from the run.def. If set to FALSE, and the multi-layer energy budget is used the profile information and canopy levels heights comes from the PGap-based processes for calculation of stand profile information in stomate. |
| 65 | |
| 66 | '''ok_mleb_history_file:''' |
| 67 | Flag that controls the writing of an output file with the multi-layer energy simulations (true/false). Note that this is a large file and writing slows down the code. |
| 68 | |
| 69 | The multi_layer_control set these flags in the following manner: |
| 70 | |
| 71 | 1 – single layer, using the multi layer energy scheme |
| 72 | ok_hydrol_arch = .TRUE.; |
| 73 | ok_gs_feedback = .TRUE.; |
| 74 | ok_mleb = .TRUE.; |
| 75 | ok_impose_can_structure = .TRUE.; |
| 76 | ok_mleb_history_file = .FALSE. |
| 77 | 2 – multi-layer energy budget |
| 78 | ok_hydrol_arch = .TRUE.; |
| 79 | ok_gs_feedback = .TRUE.; |
| 80 | ok_mleb = .TRUE.; |
| 81 | ok_impose_can_structure = .FALSE.; |
| 82 | ok_mleb_history_file = .FALSE. |
| 83 | 3 – user specific as set in the run.def |
| 84 | All read from the run.def |
| 85 | 4 – single layer, using the original enerbil scheme |
| 86 | ok_hydrol_arch = .FALSE.; |
| 87 | ok_gs_feedback = .FALSE.; |
| 88 | ok_mleb = .FALSE.; |
| 89 | ok_impose_can_structure = .FALSE.; |
| 90 | ok_mleb_history_file = .FALSE. |
| 91 | |