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
temperatureandpressuremust be provided) or anoceobject, 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 withinctdobjects (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