Changes between Version 39 and Version 40 of Branches/MergeOCN/Goll

Timestamp:

2015-06-05T14:55:12+02:00 (9 years ago)

Author:

dgoll

Comment:

--

Branches/MergeOCN/Goll

@@ -384 +384 @@
 }}}
 
+== Mass conservation issue: stomate_turnover.f90 ===
+The fraction of nitrogen recycled is calculated wrong when updating circ_class_biomass (*,ilabile,introgen). This can be fixed by:
+{{{
+                !nitrogen
+                circ_class_biomass(:,ivm,icirc,ileaf,initrogen) = circ_class_biomass(:,ivm,icirc,ileaf,initrogen) * &
+                     (un - leaf_frac(:,ivm,ilage) * turnover_rate(:))
+                circ_class_biomass(:,ivm,icirc,iroot,initrogen) = circ_class_biomass(:,ivm,icirc,iroot,initrogen) * &
+                     (un - leaf_frac(:,ivm,ilage) * turnover_rate(:))
+                circ_class_biomass(:,ivm,icirc,ifruit,initrogen) = circ_class_biomass(:,ivm,icirc,ifruit,initrogen) * &
+                     (un - leaf_frac(:,ivm,ilage) * turnover_rate(:))
+
+                ! add the recycled nutrients to labile pool:      
+                circ_class_biomass(:,ivm,icirc,ilabile,initrogen) =                                        &
+                                       circ_class_biomass(:,ivm,icirc,ilabile,initrogen)                   &
+                                        +  circ_class_biomass(:,ivm,icirc,ileaf,initrogen) *               & ! add from leaves
+                                        !DSGdebug_12  (un - leaf_frac(:,ivm,ilage) * turnover_rate(:)  * recycle_leaf) &
+                                          (  leaf_frac(:,ivm,ilage) * turnover_rate(:)  * recycle_leaf) &
+                                        +  circ_class_biomass(:,ivm,icirc,iroot,initrogen) *               & ! add from roots
+                                        !DSGdebug_12 (un - leaf_frac(:,ivm,ilage) * turnover_rate(:)  * recycle_root) 
+                                          (  leaf_frac(:,ivm,ilage) * turnover_rate(:)  * recycle_root)
+
+}}}
+
+
+{{{
+                !nitrogen
+                circ_class_biomass(:,ivm,icirc,ileaf,initrogen) = circ_class_biomass(:,ivm,icirc,ileaf,initrogen) * &
+                     (un - leaf_frac(:,ivm,ilage) * turnover_rate(:))
+                circ_class_biomass(:,ivm,icirc,iroot,initrogen) = circ_class_biomass(:,ivm,icirc,iroot,initrogen) * &
+                     (un - leaf_frac(:,ivm,ilage) * turnover_rate(:))
+                circ_class_biomass(:,ivm,icirc,ifruit,initrogen) = circ_class_biomass(:,ivm,icirc,ifruit,initrogen) * &
+                     (un - leaf_frac(:,ivm,ilage) * turnover_rate(:))
+
+                ! add the recycled nutrients to labile pool:      
+                circ_class_biomass(:,ivm,icirc,ilabile,initrogen) =                                        &
+                                       circ_class_biomass(:,ivm,icirc,ilabile,initrogen)                   &
+                                        +  circ_class_biomass(:,ivm,icirc,ileaf,initrogen) *               & ! add from leaves
+                                        !DSGdebug_12  (un - leaf_frac(:,ivm,ilage) * turnover_rate(:)  * recycle_leaf) &
+                                          (  leaf_frac(:,ivm,ilage) * turnover_rate(:)  * recycle_leaf) &
+                                        +  circ_class_biomass(:,ivm,icirc,iroot,initrogen) *               & ! add from roots
+                                        !DSGdebug_12 (un - leaf_frac(:,ivm,ilage) * turnover_rate(:)  * recycle_root) 
+                                          (  leaf_frac(:,ivm,ilage) * turnover_rate(:)  * recycle_root)
+}}}
+
+
 = '''Mass conservation checks''' =
 Mass closure given by: