diff --git a/src/Physics/Surface/SEB/seb.jl b/src/Physics/Surface/SEB/seb.jl
index 0db86369721863149c4de5b689b254e3ee27fcc3..186be78f51897ee2366593d909be9d14caf7d363 100644
--- a/src/Physics/Surface/SEB/seb.jl
+++ b/src/Physics/Surface/SEB/seb.jl
@@ -19,7 +19,7 @@ struct Iterative <: SolutionScheme end
 Equations from Westermann 2016, but use nonlinear solver.
 """
 Base.@kwdef struct Numerical{TSolver} <: SolutionScheme
-    solver::TSolver = SimpleNewtonRaphson()
+    solver::TSolver = SimpleNewtonRaphson(autodiff=AutoFiniteDiff())
 end
 Flatten.flattenable(::Type{<:Numerical}, ::Type) = false