Interception from modified Rutter model

- The `interception` output (interception loss by evaporation) from the modified Rutter interception model should be used in `SBM`, also as input for the computation of total actual evapotranspiration `actevap`. - Update of docs and changelog.
Deltares · Sep 1, 2023 · c7c615e · c7c615e
1 parent d02f900
commit c7c615e
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 3 deletions.
diff --git a/docs/src/changelog.md b/docs/src/changelog.md
@@ -14,6 +14,10 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
   amounts. And the first drainage amount `dd`, controlled by a change over time in canopy
   storage capacity `cmax`, should not be subtracted from precipitation to compute net
   interception.
+- The `netinterception` (net interception) computed by the modified Rutter interception
+  model was stored as `interception` in `SBM`, while this should be the `interception`
+  (interception loss by evaporation) output of the modified Rutter interception model. The
+  `interception` of `SBM` is used to compute the total actual evapotranspiration `actevap`.
 
 ## v0.7.1 - 2023-06-30
 

diff --git a/docs/src/model_docs/params_vertical.md b/docs/src/model_docs/params_vertical.md
@@ -76,7 +76,7 @@ specific_leaf = "Sl"
 | `pottrans_soil` | interception subtracted from potential evaporation) | mm Δt``^{-1}`` | - |
 | `transpiration` | transpiration | mm Δt``^{-1}`` | - |
 | `ae_ustore` | actual evaporation from unsaturated store | mm Δt``^{-1}`` | - |
-| `interception` | interception | mm Δt``^{-1}`` | - |
+| `interception` | interception loss by evaporation | mm Δt``^{-1}`` | - |
 | `soilevap` | total soil evaporation from unsaturated and saturated store | mm Δt``^{-1}`` | - |
 | `soilevapsat` | soil evaporation from saturated store | mm Δt``^{-1}`` | - |
 | `actcapflux` | actual capillary rise | mm Δt``^{-1}`` | - |

diff --git a/src/sbm.jl b/src/sbm.jl
@@ -84,7 +84,7 @@
     transpiration::Vector{T}
     # Actual evaporation from unsaturated store [mm Δt⁻¹]
     ae_ustore::Vector{T}
-    # Interception [mm Δt⁻¹]
+    # Interception loss by evaporation [mm Δt⁻¹]
     interception::Vector{T}
     # Soil evaporation from unsaturated and saturated store [mm Δt⁻¹]
     soilevap::Vector{T}
@@ -650,7 +650,6 @@ function update_until_snow(sbm::SBM, config)
                     cmax,
                 )
             pottrans_soil = max(0.0, leftover)  # now in mm
-            interception = netinterception
         end
 
         if modelsnow