Compute \(\sigma_\theta\), the density of seawater, minus 1000 kg/m\(^3\).
Usage
swSigma(
salinity,
temperature = NULL,
pressure = NULL,
longitude = NULL,
latitude = NULL,
eos = getOption("oceEOS", default = "gsw")
)
Arguments
- salinity
either practical salinity (in which case
temperature
andpressure
must be provided) or anoce
object, in which casesalinity
,temperature
(in the ITS-90 scale; see next item), etc. are inferred from the object, ignoring the other parameters, if they are supplied.- temperature
in-situ temperature (\(^\circ\)C), defined on the ITS-90 scale. This scale is used by GSW-style calculation (as requested by setting
eos="gsw"
), and is the value contained withinctd
objects (and probably most other objects created with data acquired in the past decade or two). Since the UNESCO-style calculation is based on IPTS-68, the temperature is converted within the present function, usingT68fromT90()
.- pressure
pressure (dbar)
- longitude
longitude of observation (only used if
eos="gsw"
; see “Details”).- latitude
latitude of observation (only used if
eos="gsw"
; see “Details”).- eos
equation of state, either
"unesco"
(references 1 and 2) or"gsw"
(references 3 and 4).
Value
Density anomaly (kg/m\(^3\)), as computed with swRho()
, minus-
1000 kg/m\(^3\).
References
See citations provided in the swRho()
documentation.
See also
Other functions that calculate seawater properties:
T68fromT90()
,
T90fromT48()
,
T90fromT68()
,
computableWaterProperties()
,
locationForGsw()
,
swAbsoluteSalinity()
,
swAlpha()
,
swAlphaOverBeta()
,
swBeta()
,
swCSTp()
,
swConservativeTemperature()
,
swDepth()
,
swDynamicHeight()
,
swLapseRate()
,
swN2()
,
swPressure()
,
swRho()
,
swRrho()
,
swSCTp()
,
swSR()
,
swSTrho()
,
swSigma0()
,
swSigma1()
,
swSigma2()
,
swSigma3()
,
swSigma4()
,
swSigmaT()
,
swSigmaTheta()
,
swSoundAbsorption()
,
swSoundSpeed()
,
swSpecificHeat()
,
swSpice()
,
swSpiciness0()
,
swSpiciness1()
,
swSpiciness2()
,
swSstar()
,
swTFreeze()
,
swTSrho()
,
swThermalConductivity()
,
swTheta()
,
swViscosity()
,
swZ()
Examples
library(oce)
swSigma(35, 13, 1000, longitude = 300, latitude = 30, eos = "gsw") # 30.82374
#> [1] 30.82374
swSigma(35, T90fromT68(13), 1000, eos = "unesco") # 30.8183
#> [1] 30.8183