diff --git a/src/solvers/dgsem_p4est/dg_2d_manifold_in_3d_cartesian.jl b/src/solvers/dgsem_p4est/dg_2d_manifold_in_3d_cartesian.jl
index 59b02c0f..169fa09c 100644
--- a/src/solvers/dgsem_p4est/dg_2d_manifold_in_3d_cartesian.jl
+++ b/src/solvers/dgsem_p4est/dg_2d_manifold_in_3d_cartesian.jl
@@ -8,7 +8,7 @@ function Trixi.rhs!(du, u, t,
                     dg::DG, cache) where {Source}
 
     # Reset du
-    Trixi.@trixi_timeit Trixi.timer() "reset ∂u/∂t" Trixi.reset_du!(du, dg, cache)
+    Trixi.@trixi_timeit Trixi.timer() "reset ∂u/∂t" Trixi.set_zero!(du, dg, cache)
 
     # Calculate volume integral
     Trixi.@trixi_timeit Trixi.timer() "volume integral" begin
diff --git a/src/solvers/dgsem_p4est/dg_2d_manifold_in_3d_covariant.jl b/src/solvers/dgsem_p4est/dg_2d_manifold_in_3d_covariant.jl
index f9b03b9c..bb78e40e 100644
--- a/src/solvers/dgsem_p4est/dg_2d_manifold_in_3d_covariant.jl
+++ b/src/solvers/dgsem_p4est/dg_2d_manifold_in_3d_covariant.jl
@@ -11,7 +11,7 @@ function Trixi.rhs!(du, u, t,
                     boundary_conditions, source_terms::Source,
                     dg::DG, cache) where {Source}
     # Reset du
-    Trixi.@trixi_timeit Trixi.timer() "reset ∂u/∂t" Trixi.reset_du!(du, dg, cache)
+    Trixi.@trixi_timeit Trixi.timer() "reset ∂u/∂t" Trixi.set_zero!(du, dg, cache)
 
     # Calculate volume integral
     Trixi.@trixi_timeit Trixi.timer() "volume integral" begin