diff --git a/src/SurfaceFluxes.jl b/src/SurfaceFluxes.jl
index 3cd12e5e..8f90cb1f 100644
--- a/src/SurfaceFluxes.jl
+++ b/src/SurfaceFluxes.jl
@@ -291,7 +291,12 @@ function surface_conditions(
     noniterative_stable_sol::Bool = true,
 ) where {FT}
     uft = SFP.universal_func_type(param_set)
-    L_MO = obukhov_length(param_set, sc, uft, scheme; tol, tol_neutral, maxiter, soltype, noniterative_stable_sol)
+    # FIXME: Workaround for julia 1.10 and GPUs.
+    if sc isa Coefficients || sc isa FluxesAndFrictionVelocity
+        L_MO = obukhov_length(param_set, sc, uft, scheme)
+    else
+        L_MO = obukhov_length(param_set, sc, uft, scheme; tol, tol_neutral, maxiter, soltype, noniterative_stable_sol)
+    end
     ustar = compute_ustar(param_set, L_MO, sc, uft, scheme)
     Cd = momentum_exchange_coefficient(param_set, L_MO, sc, uft, scheme, tol_neutral)
     Ch = heat_exchange_coefficient(param_set, L_MO, sc, uft, scheme, tol_neutral)