oce (version 0.1-80)

plot.pt: Plot pT (pressure-temperature) recorder data

Description

Plot pT temperature-depth recorder data

Usage

## S3 method for class 'pt':
plot(x,
     which=1:4, title=deparse(substitute(x)), adorn=NULL,
     tlim, plim, Tlim, 
     xlab, ylab,
     draw.time.range=getOption("oce.draw.time.range"),
     mgp=getOption("oce.mgp"),
     mar=c(mgp[1], mgp[1]+1, 1, 1.5),
     ...)

Arguments

x

pt object, typically result of read.pt.

which

if a subset of the plots is required, specify a subset of the numbers 1:4. These are codes for the panels, which will be drawn from top-left to bottom-right (see Details).

title

character string to be used in the text-summary panel (which=2)

adorn

list of expressions to be executed for the panels in turn, e.g. to adorn the plots. If the number matches the number of panels, then the strings are applied to the appropriate panels, as they are drawn from top-left to bottom-right. If only

tlim

optional limits for time axis. If not provided, the value will be inferred from the data.

plim

optional limits for pressure axis. If not provided, the value will be inferred from the data. (It is helpful to specify this, if the auto-scaled value will be inappropriate, e.g. if more lines are to be added later.)

Tlim

optional limits for temperature axis. If not provided, the value will be inferred from the data. (It is helpful to specify this, if the auto-scaled value will be inappropriate, e.g. if more lines are to be added later.)

xlab

optional label for x axis.

ylab

optional label for y axis.

draw.time.range

boolean that applies to panels with time as the horizontal axis, indicating whether to draw the time range in the top-left margin of the plot.

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

...

optional arguments passed to plotting functions.

Details

Several plots are available: which=1 for a timeseries plot of temperature, which=2 for textual information about the dataset, which=3 for a timeseries plot of pressure, and which=4 for a scatterplot of the temperature profile, . You may supply a list of which values, to get multi-panel plots.

See Also

See read.pt for more information on dealing with the pt data type, including examples.

Examples

data(pt)
plot(pt, which=c(1,3), 
    adorn=expression({abline(v=as.POSIXct('2008-07-04 20:00:00',tz='UTC'),col='red')}))