| 98 | | |
| | 98 | === supplement === |
| | 99 | {{{ |
| | 100 | ! Calculate the C that needs to be distributed to the |
| | 101 | ! labile pool. The fraction is proportional to the ratio |
| | 102 | ! between the total allocatable biomass and the unallocated |
| | 103 | ! biomass per tree (b_inc now contains the unallocated |
| | 104 | ! biomass). At the end of the allocation scheme bm_alloc_tot |
| | 105 | ! is substracted from the labile biomass pool to update the |
| | 106 | ! biomass pool (biomass(:,:,ilabile) = biomass(:,:,ilabile) - |
| | 107 | ! bm_alloc_tot(:,:)). At that point, the scheme puts the |
| | 108 | ! unallocated b_inc into the labile pool. What we |
| | 109 | ! want is that the unallocated fraction is removed from |
| | 110 | ! ::bm_alloc_tot such that only the allocated C is removed |
| | 111 | ! from the labile pool. b_inc_tot will be moved back into |
| | 112 | ! the labile pool in 5.2.11 |
| | 113 | |
| | 114 | !DSGdebug_01a: we do not have to add the b_inc_tot to ilabile as |
| | 115 | ! we do not have subtracted it due to DSGdebug_01 |
| | 116 | |
| | 117 | bm_alloc_tot(ipts,j) = bm_alloc_tot(ipts,j) - b_inc_tot |
| | 118 | !DSGdebug_01a biomass(ipts,j,ilabile,icarbon) = & |
| | 119 | !DSGdebug_01a |
| | 120 | }}} |
