Changes between Version 57 and Version 58 of Branches/MergeOCN/Goll

Timestamp:

2015-06-16T14:23:46+02:00 (9 years ago)

Author:

dgoll

Comment:

--

Branches/MergeOCN/Goll

@@ -689 +689 @@
                                                       dt / leaffall(ivm) * recycle_root(ivm)
 
-                circ_class_biomass(:,ivm,icirc,ilabile,iphosphorus) =                                      &
-                                                  circ_class_biomass(:,ivm,icirc,ilabile,iphosphorus)      &
-                                                  +  circ_class_biomass(:,ivm,icirc,ileaf,iphosphorus) *   & ! add from leaves
-                     !DSGdebug_12                                (un - dt / leaffall(ivm) * p_recycle_leaf)            &
-                                                     ( dt / leaffall(ivm)) * p_recycle_leaf(ivm)           &
-                                                  +  circ_class_biomass(:,ivm,icirc,iroot,iphosphorus) *   & ! add from roots
-                     !DSGdebug_12                                (un - dt / leaffall(ivm) * p_recycle_root)           
-                                                     ( dt / leaffall(ivm)) * p_recycle_root(ivm)
-
-                     !DSGdebug_12 
 
 
@@ -712 +702 @@
                                                                       (un - dt / leaffall(ivm))
 
-                circ_class_biomass(:,ivm,icirc,ileaf,iphosphorus) = circ_class_biomass(:,ivm,icirc,ileaf,iphosphorus) * &
-                                                                      (un - dt / leaffall(ivm))
-                circ_class_biomass(:,ivm,icirc,iroot,iphosphorus) = circ_class_biomass(:,ivm,icirc,iroot,iphosphorus) * &
-                                                                      (un - dt / leaffall(ivm))
+
              !DSGdebug_16   ! add the recycled nutrients to labile pool:      
              !DSGdebug_16   circ_class_biomass(:,ivm,icirc,ilabile,initrogen) =                                      &
@@ -726 +713 @@
              !DSGdebug_16                                         dt / leaffall(ivm) * recycle_root(ivm)
 
-             !DSGdebug_16   circ_class_biomass(:,ivm,icirc,ilabile,iphosphorus) =                                      &
-             !DSGdebug_16                                     circ_class_biomass(:,ivm,icirc,ilabile,iphosphorus)      &
-             !DSGdebug_16                                     +  circ_class_biomass(:,ivm,icirc,ileaf,iphosphorus) *   & ! add from leaves
-             !DSGdebug_16        !DSGdebug_12                                (un - dt / leaffall(ivm) * p_recycle_leaf)            &
-             !DSGdebug_16                                        ( dt / leaffall(ivm)) * p_recycle_leaf(ivm)           &
-             !DSGdebug_16                                     +  circ_class_biomass(:,ivm,icirc,iroot,iphosphorus) *   & ! add from roots
-             !DSGdebug_16        !DSGdebug_12                                (un - dt / leaffall(ivm) * p_recycle_root)           
-             !DSGdebug_16                                        ( dt / leaffall(ivm)) * p_recycle_root(ivm)           
-             !DSGdebug_16        !DSGdebug_12 
+
+}}}
+[...]
+{{{
+       !! 6.1 For deciduous trees: next to leaves, also fruits and fine roots are dropped 
+       !  For deciduous trees: next to leaves, also fruits and fine roots are dropped: fruit ::biomass(:,ivm,ifruit) 
+       !  and fine root ::biomass(:,ivm,iroot) carbon pools are set to zero.
+       IF ( is_tree(ivm) .AND. ( senescence_type(ivm) .NE. 'none' ) ) THEN
+
+          ! Check whether we shed the remaining leaves. The condition depends on biomass(ileaf)
+          ! so first calculate the sheding at the tree level. because biomass(ileaf) is not 
+          ! changed at the tree level the same statement can then be used at the stand level.
+          DO icirc = 1,ncirc
+
+             ! check whether we shed the remaining leaves
+             WHERE ( ( biomass(:,ivm,ileaf,icarbon) .GT. zero ) .AND. senescence(:,ivm) .AND. &
+                  ( biomass(:,ivm,ileaf,icarbon) .LT. (lai_initmin(ivm) / 2.)/sla(ivm) ) )
+
+                ! Tree level (g tree-1)
+                !DSGdebug_16: the recycling has to be done before the source
+                !pools are set to zero, as these pool are part of the equation:
+
+                ! relocate nutrient 
+                circ_class_biomass(:,ivm,icirc,ilabile ,initrogen)  =                           &
+                !DSGdebug_14
+                         circ_class_biomass(:,ivm,icirc,ilabile ,initrogen)                     &
+
+
+                !DSGdebug_14
+                         + circ_class_biomass(:,ivm,icirc,ileaf ,initrogen) * recycle_leaf(ivm) &
+                         + circ_class_biomass(:,ivm,icirc,iroot ,initrogen) * recycle_root(ivm)
+
+
+
+                circ_class_biomass(:,ivm,icirc,ileaf,icarbon)       = zero
+                circ_class_biomass(:,ivm,icirc,iroot,icarbon)       = zero
+                circ_class_biomass(:,ivm,icirc,ifruit,icarbon)      = zero
+                circ_class_biomass(:,ivm,icirc,ileaf ,initrogen)    = zero
+                circ_class_biomass(:,ivm,icirc,iroot ,initrogen)    = zero
+                circ_class_biomass(:,ivm,icirc,ifruit,initrogen)    = zero
+
+
+               !DSGdebug_16 circ_class_biomass(:,ivm,icirc,ilabile ,initrogen)  =                           &
+               !DSGdebug_16 !DSGdebug_14
+               !DSGdebug_16          circ_class_biomass(:,ivm,icirc,ilabile ,initrogen)                     &
+               !DSGdebug_16 !DSGdebug_14
+               !DSGdebug_16          + circ_class_biomass(:,ivm,icirc,ileaf ,initrogen) * recycle_leaf(ivm) &
+               !DSGdebug_16          + circ_class_biomass(:,ivm,icirc,iroot ,initrogen) * recycle_root(ivm)
+
+
 }}}