From 516670a48901251bfcc968e20f7eed0cd6267c1a Mon Sep 17 00:00:00 2001
From: Ryan Knox <rgknox@gmail.com>
Date: Wed, 11 Jul 2018 16:26:25 -0700
Subject: [PATCH] Updated cumulative TVAI function to estimate the mid-point of
 layer of interest

---
 biogeochem/FatesAllometryMod.F90 | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/biogeochem/FatesAllometryMod.F90 b/biogeochem/FatesAllometryMod.F90
index 12934256..8ced4209 100644
--- a/biogeochem/FatesAllometryMod.F90
+++ b/biogeochem/FatesAllometryMod.F90
@@ -2155,7 +2155,8 @@ function CumulativeLayerTVAI(icanlayer,     &
        call endrun(msg=errMsg(sourcefile, __LINE__))
     end if
 
-    cum_tvai = cum_tvai + tvai0 + tvai
+    ! Apply the cumulative TVAI at the current layer mid-point (ie 0.5*)
+    cum_tvai = cum_tvai + tvai0 + 0.5*tvai
     
     return
   end function CumulativeLayerTVAI