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")