swRho(salinity, temperature=NULL, pressure=NULL,
longitude, latitude, eos=getOption("oceEOS", default="gsw"))temperature
and pressure must be provided) or an oce object, in
which case salinity, temperature (in the ITS-90 scale; see
eos="gsw"), and is the value contained within
ctd objects (and probably mosteos="gsw"; see
eos="gsw"; see
"unesco" [1,2] or "gsw"
[3,4].sw* functions took temperature to be in
IPTS-68 units. As GSW capabilities were added in early 2015, the assumed unit
of temperature was taken to be ITS-90. This change means that old code
has to be modified, by replacing e.g. swRho(S, T, p) with
swRho(S, T90fromT68(T), p). At typical oceanic values, the difference
between the two scales is a few millidegrees.oce object, then values for
salinity, etc., are extracted from it, and used for the calculation, and the
corresponding arguments to the present function are ignored. If eos="unesco", the density is calculated using the UNESCO equation of
state for seawater [1,2], and if eos="gsw", the GSW formulation [3,4]
is used.
[2] Gill, A.E., 1982. Atmosphere-ocean Dynamics, Academic Press, New York, 662 pp.
[3] IOC, SCOR, and IAPSO (2010). The international thermodynamic equation of seawater-2010: Calculation and use of thermodynamic properties. Technical Report 56, Intergovernmental Oceanographic Commission, Manuals and Guide.
[4] McDougall, T.J. and P.M. Barker, 2011: Getting started with TEOS-10 and the Gibbs Seawater (GSW) Oceanographic Toolbox, 28pp., SCOR/IAPSO WG127, ISBN 978-0-646-55621-5.
swSigma0 (and equivalents at other pressure horizons),
swSigmaT, and
swSigmaTheta.library(oce)
rho <- swRho(35, 13, 1000)Run the code above in your browser using DataLab