| plot.TS {oce} | R Documentation |
Plot temperature-salinity diagram for seawater (CTD) data.
plot.TS(x,
rho.levels = 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 = 21,
rotate.rho.labels = FALSE,
connect.points = FALSE,
xlab,
ylab,
Slim,
Tlim,
mgp=getOption("oce.mgp"),
mar=c(mgp[1]+5/4,mgp[1]+1,mgp[1],mgp[1]),
lwd.rho=par("lwd"), lty.rho=par("lty"),
...)
x |
A cdt object, e.g. as read by read.ctd. |
rho.levels |
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 of symbols on graph. |
pch |
code for symbols on graph. |
rotate.rho.labels |
if TRUE, labels in right-hand margin are written vertically |
connect.points |
if TRUE, the points are connected with line segments, in sequence |
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. |
Creates a temperature-salinity plot for a CTD cast, with labeled isopycnals.
None.
Dan Kelley
summary.ctd summarizes the information, while
read.ctd scans it from a file.
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
plot.TS(profile)
title("Mid-Atlantic at 36N")