diff --git a/veros/diagnostics/overturning.py b/veros/diagnostics/overturning.py
index da7345e6..2fba05d5 100644
--- a/veros/diagnostics/overturning.py
+++ b/veros/diagnostics/overturning.py
@@ -121,7 +121,7 @@ def output(self, state):
 
         ovt_vs = self.variables
 
-        mean_variables = ("trans", "vsf_iso", "vsf_depth")
+        mean_variables = ("trans", "vsf_iso", "vsf_depth", "bolus_iso", "bolus_depth")
 
         if ovt_vs.nitts > 0:
             for var in mean_variables: