Changes between Version 19 and Version 20 of Branches/MergeOCN/Goll

Timestamp:

2015-05-26T10:30:56+02:00 (10 years ago)

Author:

dgoll

Comment:

--

Branches/MergeOCN/Goll

@@ -164 +164 @@
 }}}
 
-== '''4. Negative biomass: stomate_growth_fun_all.f90''' ==
+== '''5. Negative biomass: stomate_growth_fun_all.f90''' ==
 In case the variable "residual" is not zero and nutrient availability reduces bm_alloc_tot to a value smaller than the value of residual, the allocated biomass gets negative.
 This can be fixed by correcting the biomass to be allocated by the residual before calculating nutrient limitation. 
+{{{
+                ! Move the unallocated carbon back into the labile pool
+                biomass(ipts,j,ilabile,icarbon) = &
+                     biomass(ipts,j,ilabile,icarbon) + residual(ipts,j)
+
+             ENDIF
+
+!DSGdebug_05
+             ! correct the biomass to be allocated by the residual
+             ! this has to be done here because we need the actual biomass to be
+             ! allocated to compute the nutrient demands
+             bm_alloc_tot(ipts,j) = bm_alloc_tot(ipts,j) - residual(ipts,j)
+!DSGdebug_05
+
+
+}}}
+{{{
+!DSGdebug_05
+!          bm_alloc(:,j,k,icarbon) = f_alloc(:,j,k) * (bm_alloc_tot(:,j) - residual(:,j))
+           bm_alloc(:,j,k,icarbon) = f_alloc(:,j,k) * bm_alloc_tot(:,j)  ! residual was already removed from bm_alloc_tot
+!DSGdebug_05
+}}}