| 98 | |
| 99 | |
| 100 | |
| 101 | == '''2. Mass conservation issue: stomate_phenology.f90''' == |
| 102 | The nutrient demand must be calculated AFTER the final Cl_init and Cr_init are known. |
| 103 | This can be fixed by: |
| 104 | |
| 105 | {{{ |
| 106 | ! The biomass available to use is set to be the minimum of the biomass of |
| 107 | ! the labile pool (if carbon not taken from the atmosphere), and |
| 108 | ! the wanted biomass. |
| 109 | bm_use(i) = MIN( biomass(i,j,ilabile,icarbon) + biomass(i,j,icarbres,icarbon), & |
| 110 | bm_wanted(i) ) |
| 111 | |
| 112 | !DSGdebug_02 ! the nutrients need to support the biomass: |
| 113 | !DSGdebug_02 bm_wanted_n(i) = (Cl_init + Cr_init*fcn_root(j))/cn_leaf_prescribed(j) |
| 114 | }}} |
| 115 | [...] |
| 116 | {{{ |
| 117 | ! In case nitrogen or phosphorus is not sufficiently available |
| 118 | ! downregulate new leaf biomass to respect leaf stoichiometry; |
| 119 | ! DSG: this violates the ratio used to calculate the |
| 120 | ! leave-root-sapwood relationships: is this OK? |
| 121 | |
| 122 | !DSGdebug_02: moved after Cl_init and Cr_init are updated |
| 123 | ! the nutrients need to support the biomass: |
| 124 | bm_wanted_n(i) = (Cl_init + Cr_init*fcn_root(j))/cn_leaf_prescribed(j) |
| 125 | !DSGdebug_02 |
| 126 | }}} |
| 127 | |