oce (version 0.9-21)

swSigmaTheta: Seawater potential density anomaly

Description

Compute \(\sigma_\theta\), the potential density of seawater, minus 1000 kg/m\(^3\).

Usage

swSigmaTheta(salinity, temperature = NULL, pressure = NULL,
  referencePressure = 0, longitude = 300, latitude = 30,
  eos = getOption("oceEOS", default = "gsw"))

Arguments

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 [\(^\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 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]

referencePressure

The reference pressure, in 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].

Value

Potential density anomaly [kg/m\(^3\)], defined as \(\sigma_\theta=\rho(S,\theta(S,t,p),0\) - 1000 kg/m\(^3\).

Details

If the first argument is an oce object, then salinity, etc., are extracted from it, and used for the calculation instead of any values provided in the other arguments.

References

See citations provided in the swRho documentation.

See Also

Other functions that calculate seawater properties: T68fromT90, T90fromT48, T90fromT68, swAbsoluteSalinity, swAlphaOverBeta, swAlpha, swBeta, swCSTp, swConservativeTemperature, swDepth, swDynamicHeight, swLapseRate, swN2, swPressure, swRho, swRrho, swSCTp, swSTrho, swSigma0, swSigma1, swSigma2, swSigma3, swSigma4, swSigmaT, swSigma, swSoundAbsorption, swSoundSpeed, swSpecificHeat, swSpice, swTFreeze, swTSrho, swThermalConductivity, swTheta, swViscosity, swZ

Examples

# NOT RUN {
swSigmaTheta(35, 13, 1000)             # 26.42514 (gsw)
swSigmaTheta(35, 13, 1000, eos="unesco") # 26.4219

# }