Plot a profile, showing variation of some quantity (or quantities) with
pressure, using the oceanographic convention of putting lower pressures
nearer the top of the plot. This works for any
oce object that has a
pressure column in its
The colors (
col.salinity, etc.) are only used if two profiles appear
on a plot.
plotProfile( x, xtype = "salinity+temperature", ytype = "pressure", eos = getOption("oceEOS", default = "gsw"), lty = 1, xlab = NULL, ylab = NULL, col = "black", col.salinity = "darkgreen", col.temperature = "red", col.rho = "blue", col.N2 = "brown", col.dpdt = "darkgreen", col.time = "darkgreen", pt.bg = "transparent", grid = TRUE, col.grid = "lightgray", lty.grid = "dotted", Slim, Clim, Tlim, densitylim, N2lim, Rrholim, dpdtlim, timelim, plim, xlim, ylim, lwd = par("lwd"), xaxs = "r", yaxs = "r", cex = 1, pch = 1, useSmoothScatter = FALSE, df, keepNA = FALSE, type = "l", mgp = getOption("oceMgp"), mar, add = FALSE, inset = FALSE, debug = getOption("oceDebug", 0), ... )
a ctd object.
Item(s) plotted on the x axis, either a vector of length equal
to that of
pressure in the
or a text code from the list below.
"salinity" Profile of salinity.
"conductivity" Profile of conductivity.
"temperature" Profile of in-situ temperature.
"theta" Profile of potential temperature.
"density" Profile of density (expressed as \(\sigma_\theta\)).
"index" Index of sample (useful for working with
"salinity+temperature" Profile of salinity and temperature within a single axis frame.
"N2" Profile of square of buoyancy frequency \(N^2\),
an optional argument setting of
df=length(x[["pressure"]])/4 to do
"density+N2" Profile of sigma0 and
the square of buoyancy frequency within a single axis frame.
"density+dpdt" Profile of sigma0 and dP/dt for the
sensor. The latter is useful in indicating problems with the deployment.
It is calculated by first differencing pressure and then using a smoothing
spline on the result (to avoid grid-point wiggles that result because the
SBE software only writes 3 decimal places in pressure). Note that dP/dt may
be off by a scale factor; this should not be a problem if there is a
time column in the
data slot, or a
sample.rate in the
"sigma4" Profile of potential density referenced
to 0dbar (i.e. the surface), 1000dbar, 2000dbar, 3000dbar, and 4000dbar.
"spice" Profile of spice.
"Rrho" Profile of Rrho, defined in the diffusive sense.
"RrhoSF" Profile of Rrho, defined in the salt-finger sense.
variable to use on y axis. The valid choices are:
"pressure" (the default),
equation of state to be used, either
line type for the profile.
optional label for x axis (at top of plot). If not
provided, the value of
xtype is used (along with the units
of the variable) to construct a label. For user-supplied
the unit (if any) must be specified as in e.g.
optional label for y axis. See
xlab for a note on
ylab="" prevents labelling the axis.
color for a general profile.
color for salinity profile (see “Details”).
color for temperature (see “Details”).
color for density (see “Details”).
color for square of buoyancy frequency (see “Details”).
color for dP/dt.
color for delta-time.
inside color for symbols with
pch in 21:25
logical, set to
TRUE to get a grid.
color for grid.
line type for grid.
Optional limit for S axis
Optional limit for conductivity axis
Optional limit for T axis
Optional limit for density axis
Optional limit for N2 axis
Optional limit for Rrho axis
Optional limit for dp/dt axis
Optional limit for delta-time axis
Optional limit for pressure axis, ignored unless
ytype=="pressure", in which case it takes precedence over
Optional limit for x axis, which can apply to any plot type.
This is ignored if the plotted x variable is something for which a limit
may be specified with an argument, e.g.
xlim is ignored for a
salinity profile, because
Slim ought to be given in such a case.
Optional limit for y axis, which can apply to any plot type,
although is overridden by
lwd value for data line
xaxs to use
yaxs to use
size to be used for plot symbols (see
code for plotting symbol (see
optional argument, passed to
swN2() if provided, and
if a plot using \(N^2\) is requested.
type of plot to draw, using the same scheme as
Four-element numerical value to be used to set the plot
margins, with a call to par
(mar) prior to the plot.
If this is not supplied, a reasonable default will be set up.
A logical value that controls whether to add to an existing plot. (It
makes sense to use
add=TRUE in the
panel argument of a
coplot(), for example.)
A logical value indicating whether to draw an inset plot.
Setting this to
TRUE will prevent the present function from adjusting
the margins, which is
necessary because margin adjustment is the basis for the method used by
a flag that turns on debugging. Set to 1 to get a moderate amount of debugging information, or to 2 to get more.
optional arguments passed to other functions. A common example
is to set
df, for use in
Other functions that plot oce data:
Other things related to ctd data: