R: Seawater density anomaly

swSigma {oce}

R Documentation

Seawater density anomaly

Compute sigma, the density of seawater, minus 1000 kg/m^3.

swSigma(salinity, temperature = NULL, pressure = NULL,
  longitude = NULL, latitude = NULL, eos = getOption("oceEOS",
  default = "gsw"))

`salinity`	either practical salinity (in which case `temperature` and `pressure` must be provided) or an `oce` object, in which case `salinity`, `temperature` (in the ITS-90 scale; see next item), etc. are inferred from the object.
`temperature`	in-situ temperature [degC], 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 within `ctd` 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, using `T68fromT90`.
`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"` [1,2] or `"gsw"` [3,4].

Density anomaly [kg/m^3], defined as swRho - 1000 kg/m^3.

Dan Kelley

See citations provided in the swRho documentation.

library(oce)
swSigma(35, 13, 1000, longitude=300, latitude=30, eos="gsw") # 30.82374
swSigma(35, T90fromT68(13), 1000, eos="unesco")    # 30.8183

[Package oce version 1.0-1 Index]