Updated cumulative TVAI function to estimate the mid-point of layer o…

…f interest
rgknox · Jul 11, 2018 · 516670a · 516670a
1 parent c0dbd7d
commit 516670a
Showing 1 changed file with 2 additions and 1 deletion.
diff --git 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