problems using the seawater equation of state

[francispoulin]

We have a minimum working example that tries to use a nonlinear equation of state, but when we define the model, we get an error saying, the following,

nested task error: MethodError: no method matching ρ′(::Float64, ::Float64, ::Float64, ::typeof(SeawaterPolynomials.TEOS10.TEOS10EquationOfState))
    Closest candidates are:
      ρ′(::Any, ::Any, ::Any, ::Any, ::Any, ::Any, ::Any) at ~/Software/Oceananigans.jl/src/BuoyancyModels/nonlinear_equation_of_state.jl:27
      ρ′(::Any, ::Any, ::Any, ::SeawaterPolynomials.BoussinesqEquationOfState{<:SeawaterPolynomials.SecondOrderSeawaterPolynomials.SecondOrderSeawaterPolynomial}) at ~/.julia/packages/SeawaterPolynomials/JiuSX/src/SecondOrderSeawaterPolynomials.jl:64
      ρ′(::Any, ::Any, ::Any, ::SeawaterPolynomials.BoussinesqEquationOfState{<:SeawaterPolynomials.TEOS10.TEOS10SeawaterPolynomial}) at ~/.julia/packages/SeawaterPolynomials/JiuSX/src/TEOS10.jl:186

We are following the documentation shown here.

Can someone help us fix this?

The code is copied below.

using Oceananigans
using SeawaterPolynomials

grid = RectilinearGrid(CPU();
                       size = (10, 10), 
                          x = ( 0, 1),
                          z = (-1, 0),
                   topology = (Periodic, Flat, Bounded))
                   
eos = SeawaterPolynomials.TEOS10EquationOfState
buoyancy = SeawaterBuoyancy(equation_of_state=eos)

model = NonhydrostaticModel(; grid, buoyancy, tracers = (:T, :S),)

[simone-silvestri]

you have to instantiate the equation of state

eos = SeawaterPolynomials.TEOS10EquationOfState()

[francispoulin]

Thank you @simone-silvestri ! That worked beautifully.