oce (version 0.7-1)

plotTS: Plot temperature-salinity diagram for seawater (CTD) data

Description

Plot temperature-salinity diagram for seawater (CTD) data.

Usage

plotTS(x,
     inSitu=FALSE,
     referencePressure=0,
     rhoLevels = 6,
     grid = TRUE,
     col.grid = "lightgray",
     lty.grid = "dotted",
     rho1000 = FALSE,
     col = par("col"),
     col.rho = "darkgray",
     cex.rho = 0.9*par("cex"),
     cex = par("cex"), pch = 1,
     rotateRhoLabels = FALSE,
     connectPoints = FALSE,
     useSmoothScatter = FALSE,
     xlab,
     ylab,
     Slim,
     Tlim,
     mgp=getOption("oceMgp"),
     mar=c(mgp[1]+1,mgp[1]+1,mgp[1],mgp[1]),
     lwd.rho=par("lwd"), lty.rho=par("lty"),
     ...)

Arguments

x

An object containing salinity and temperature data, typically a cdt object or section object.

inSitu

A boolean indicating whether to use in-situ temperature or (the default) potential temperature, calculated with reference pressure given by referencePressure.

referencePressure

reference pressure, to be used in calculating potential temperature, if inSitu is FALSE.

rhoLevels

Either a list of density levels for which to draw isopycnal lines, or a suggestion for the number of levels. In the latter case, pretty() is used to select levels.

grid

a flag that can be set to TRUE to get a grid.

col.grid

colour for grid.

lty.grid

line type for grid.

rho1000

if TRUE, label isopycnals as e.g. 1024; if FALSE, label as e.g. 24

col

colour for symbols.

col.rho

colour for isopycnal lines.

cex.rho

size of isopycnal labels.

cex

size to be used for plot symbols (see par)

pch

code for plotting symbol (see par).

rotateRhoLabels

if TRUE, labels in right-hand margin are written vertically

connectPoints

if TRUE, the points are connected with line segments, in sequence

useSmoothScatter

if TRUE, use smoothScatter to plot the points.

xlab

optional label for the x axis, with default "Salinity [PSU]".

ylab

optional label for the y axis, with default "Temperature [C]".

Slim

optional limits for salinity axis, otherwise inferred from data.

Tlim

optional limits for temperature axis, otherwise inferred from data.

mgp

3-element numerical vector to use for par(mgp), and also for par(mar), computed from this. The default is tighter than the R default, in order to use more space for the data and less for the axes.

mar

value to be used with par("mar").

lwd.rho

line width for density lines

lty.rho

line type for density lines

...

optional arguments passed to plotting functions.

Value

• None.

Details

Creates a temperature-salinity plot for a CTD cast, with labeled isopycnals.

See Also

summary.ctd summarizes the information, while read.ctd scans it from a file.

Examples

library(oce)
# Plot TS diagram for a mid-Atlantic station at 36N
# See demo(TS) for a section-spanning TS diagram
data(a03)
profile <- a03$data$station[[60]] # 124 stations in total
plotTS(profile)
title("Mid-Atlantic at 36N")