Add workaround for julia 1.10 and CUDA

As of CUDA 5.1.1, there are some issues with ClimaAtmos and Julia 1.10.
It seems that the issue is with this particular function and the error
is due to added calles to gc_frame functions that cannot be compiled on
a GPU.

For unknown reasons, the code compiles when we remove the keyword
arguments. While the issue is being investigated, we can work around it
by excluding the offending branch

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)