search

CliMA / SeawaterPolynomials.jl

Polynomials for efficiently computing the density of seawater
https://clima.github.io/SeawaterPolynomials.jl/dev
MIT License
14 stars 3 forks source link

Density perturbation is actually total density for TEOS10 #22

Closed glwagner closed 2 years ago

glwagner commented 2 years ago

For example:

julia> using SeawaterPolynomials

julia> teos10 = SeawaterPolynomials.TEOS10.TEOS10EquationOfState()
SeawaterPolynomials.BoussinesqEquationOfState{SeawaterPolynomials.TEOS10.TEOS10SeawaterPolynomial{Float64}, Int64}(SeawaterPolynomials.TEOS10.TEOS10SeawaterPolynomial{Float64}(), 1020)

julia> SeawaterPolynomials.ρ′(0, 0, 0, teos10)
999.8434119170641

should be

julia> SeawaterPolynomials.ρ′(0, 0, 0, teos10) - teos10.reference_density
-20.156588082935855

cc @simone-silvestri