sw.beta {oce} | R Documentation |
Compute beta using McDougall's (1987) algorithm. NOTE: this conflicts with base::beta, so be sure to prefix as sw.beta().
b <- sw.beta(S, t, p, is.theta=FALSE);
S |
salinity [PSU] |
t |
Temperature or potential temperature [degC], according to is.theta |
p |
pressure [dbar] |
is.theta |
Set TRUE if t is theta or FALSE if t is in-situ |
McDougall (1987). Note the use of potential temperature, not temperature.
Value in 1/psu.
Dan Kelley Dan.Kelley@Dal.Ca
McDougall, T.J. 1987. "Neutral Surfaces" Journal of Physical Oceanography vol 17 pages 1950-1964
N/A.
b <- sw.beta(40, 10, 4000, is.theta=TRUE) # 0.72088e-3 (from p1964 of McDougall 1987)