Changes between Version 1 and Version 2 of Branches/MergeOCN/Goll

Timestamp:

2015-05-07T09:02:31+02:00 (9 years ago)

Author:

dgoll

Comment:

--

Branches/MergeOCN/Goll

@@ -1 +1 @@
 = Daniel's page =
+CNP-Dev version based on MERGE-OCN revision 2567
+
+
+
+
+== '''1. Mass conservation issue: stomate_growth_fun_all.f90''' ==
+
+The carbon being allocated to biomass pools must no be substracted before nutrient limitation of allocation is computed.
+This can be fixed by:
+{{{
+          ! The calculation of ::resp_maint is solely based on the demand i.e.
+          ! given the biomass and the condition of the plant, how much should be
+          ! respired. It is not sure that this demand can be satisfied i.e. the 
+          ! calculated maintenance respiration may exceed the available carbon
+
+          !DSGdebug_01
+          deficit = zero
+          !DSGdebug_01
+             
+          IF ( bm_alloc_tot(ipts,j) - resp_maint(ipts,j) .LT. zero ) THEN
+}}}
+[...]
+{{{
+                ! Not enough carbon to pay the deficit, the individual 
+                ! is going to die at the end of this day
+                bm_alloc_tot(ipts,j) = bm_alloc_tot(ipts,j) + &
+                     biomass(ipts,j,icarbres,icarbon) 
+                biomass(ipts,j,icarbres,icarbon) = zero
+
+                ! Truncate the maintenance respiration to the available carbon
+                resp_maint(ipts,j) = bm_alloc_tot(ipts,j)
+
+                !DSGdebug_01
+                deficit = zero
+                !DSGdebug_01
+             ENDIF
+
+}}}
+[...]
+{{{
+          ! Final ::resp_maint is know
+          bm_alloc_tot(ipts,j) = bm_alloc_tot(ipts,j) - resp_maint(ipts,j)
+
+          !DSGdebug_01
+          biomass(ipts,j,ilabile,icarbon) = biomass(ipts,j,ilabile,icarbon) - &
+                                            (resp_maint(ipts,j) + deficit)
+          !DSGdebug_01
+
+          !! 3.11 Growth respiration
+          !  Calculate total growth respiration and update allocatable carbon
+          !  Growth respiration is a tax on productivity, not actual allocation
+          !  Total growth respiration has be calculated before the allocation 
+          !  takes place because the allocation itself is not linear. After 
+          !  the allocation has been calculated, growth respiration can be 
+          !  calculated for each biomass component separatly. The unit of
+          !  resp_growth is gC m-2 dt-1
+          resp_growth(ipts,j)  = frac_growthresp(j) * MAX(zero, bm_alloc_tot(ipts,j))
+          bm_alloc_tot(ipts,j) = bm_alloc_tot(ipts,j) - resp_growth(ipts,j)
+
+          !DSGdebug_01
+          biomass(ipts,j,ilabile,icarbon) = biomass(ipts,j,ilabile,icarbon) - &
+                                            resp_growth(ipts,j)
+          !DSGdebug_01
+
+
+}}}
+[...]
+{{{
+    !=======================================================
+    ! Block from OCN but I did not find it in DOFOCO ???
+    !
+    ! 5.1 retrieve allocated biomass from labile pool (nitrogen, or new allocation)
+    !
+
+    !DSGdebug_01
+    biomass(:,:,ilabile,icarbon)   = biomass(:,:,ilabile,icarbon)   - bm_alloc_tot(:,:)
+    !DSGdebug_01
+    biomass(:,:,ilabile,initrogen) = biomass(:,:,ilabile,initrogen) - n_alloc_tot(:,:)
+
+}}}