| | 149 | |
| | 150 | = '''Mass conservation checks''' = |
| | 151 | |
| | 152 | !DSG mass conservation ======================================== |
| | 153 | mass_before(:,:,:) = SUM(biomass(:,:,:,:),DIM=3) |
| | 154 | |
| | 155 | !! 5. Grow new biomass - respiration, npp and allocation |
| | 156 | |
| | 157 | ! Call the allometry based allocation (based on Sitch et al 2003 and Zaehle et al 2010) |
| | 158 | CALL growth_fun_all (npts, dt_days, veget_max, PFTpresent, & |
| | 159 | senescence, when_growthinit, moiavail_growingseason, t2m_week, & |
| | 160 | gpp_daily, gpp_week, resp_maint_part, resp_maint, & |
| | 161 | resp_growth, npp_daily, biomass, age, & |
| | 162 | leaf_age, leaf_frac, use_reserve, t_photo_stress, & |
| | 163 | lab_fac, lai_target, ind, rue_longterm, & |
| | 164 | circ_class_n, circ_class_biomass, c0_alloc, cn_leaf_season, np_leaf_season, & |
| | 165 | KF, n_uptake_daily, p_uptake_daily) |
| | 166 | |
| | 167 | !DSG mass conservation ============================================ |
| | 168 | mass_change(:,:,icarbon) = npp_daily(:,:)*dt_days |
| | 169 | mass_change(:,:,initrogen) = SUM(n_uptake_daily(:,:,:),DIM=3) |
| | 170 | mass_change(:,:,iphosphorus) = p_uptake_daily(:,:) |
