Changes between Version 50 and Version 51 of Branches/MergeOCN/Goll

Timestamp:

2015-06-15T13:11:51+02:00 (9 years ago)

Author:

dgoll

Comment:

--

Branches/MergeOCN/Goll

@@ -501 +501 @@
 }}}
 === 11. Mass conservation issue: stomate_turnover.f90 ===
-The mass of biomass and circ_class_biomass starts to diverge after leaf fall in stomate_turnover.f90. I fixed this by removing repetitions in a big; I couldn't locate the precise error in the original code.
-{{{
+The mass of biomass and circ_class_biomass starts to diverge after leaf fall in stomate_turnover.f90.
+I rewrote the code in a more fail-safe format and thereby fixed the bug. I don't know where the error is in the original piece of code.
+
+{{{
+       !! 4.3 Loop over leaf age classes
+       DO ilage = 1, nleafages
+             
+          turnover_rate(:) = zero
+
+          WHERE ( leaf_age(:,ivm,ilage) .GT. leaf_age_crit(:,ivm) / deux)
+
+             turnover_rate(:) =  MIN( 0.99_r_std, dt / ( leaf_age_crit(:,ivm) * &
+                  ( leaf_age_crit(:,ivm) / leaf_age(:,ivm,ilage) )**quatre ) )
+             
+          ENDWHERE
+                                 
+          IF (is_tree(ivm)) THEN
+
 !DSGdebug_13             ! Stand level turnover (gC m-2)
 !DSGdebug_13             ! Leaves
@@ -549 +565 @@
 !DSGdebug_13
 
-
-}}}
-by 
-{{{
-   IF (is_tree(ivm)) THEN
              !DSGsimplify
              DO ielem = 1,nelements
@@ -617 +628 @@
              ENDDO
              !DSGsimplify
-
+!DSGdebug_13
 }}}