| 131 | This piece of code leads that less biomass is allocated to leaf and root than was subtracted earlier from labile and rescarb using variable bm_use; |
| 132 | I guess this mass leak is related to the in code comment |
| 133 | {{{ |
| 134 | ! +++CHECK+++ |
| 135 | ! Cl_init + Cr_init can exceed bm_use. bm_use should be used in these equations |
| 136 | circ_class_biomass(i,j,l,ileaf,icarbon) = circ_class_biomass(i,j,l,ileaf,icarbon) + & |
| 137 | Cl_init * ( KF(i,j) * Cs_tree(l) / height(i,j,l) * circ_class_n(i,j,l) ) / & |
| 138 | SUM( KF(i,j) * Cs_tree(:) / height(i,j,:) * circ_class_n(i,j,:) ) |
| 139 | |
| 140 | circ_class_biomass(i,j,l,iroot,icarbon) = circ_class_biomass(i,j,l,iroot,icarbon) + & |
| 141 | Cr_init * ( KF(i,j) * Cs_tree(l) / height(i,j,l) * circ_class_n(i,j,l) ) / & |
| 142 | SUM( KF(i,j) * Cs_tree(:) / height(i,j,:) * circ_class_n(i,j,:) ) |
| 143 | !++++++++++++ |
| 144 | |
| 145 | |
| 146 | |
| 147 | }}} |
| 148 | |