Creates a multi-panel summary plot of data measured in a meteorological data set. cast. The panels are controlled by the which argument.

# S4 method for met
plot(x, which = 1:4, mgp, mar, tformat, debug = getOption("oceDebug"))



a met object.


list of desired plot types.

  • which=1 gives a time-series plot of temperature

  • which=2 gives a time-series plot of pressure

  • which=3 gives a time-series plot of the x (eastward) component of velocity

  • which=4 gives a time-series plot of the y (northward) component of velocity

  • which=5 gives a time-series plot of speed

  • which=6 gives a time-series plot of direction (degrees clockwise from north; note that the values returned by met[["direction"]] must be multiplied by 10 to get the direction plotted)


A 3-element numerical vector used with par("mgp") to control the spacing of axis elements. The default is tighter than the R default.


A 4-element numerical vector used with par("mar") to control the plot margins. The default is tighter than the R default.


optional argument passed to oce.plot.ts(), for plot types that call that function. (See strptime() for the format used.)


an integer specifying whether debugging information is to be printed during the processing. This is a general parameter that is used by many oce functions. Generally, setting debug=0 turns off the printing, while higher values suggest that more information be printed. If one function calls another, it usually reduces the value of debug first, so that a user can often obtain deeper debugging by specifying higher debug values.


If more than one panel is drawn, then on exit from plot.met, the value of par will be reset to the value it had before the function call. However, if only one panel is drawn, the adjustments to par made within plot.met are left in place, so that further additions may be made to the plot.


Dan Kelley


plot(met, which = 3:4)

# Wind speed and direction during Hurricane Juan
# Compare with the final figure in a white paper by Chris Fogarty
# (available at
# downloaded 2017-01-02).
t0 <- as.POSIXct("2003-09-29 04:00:00", tz = "UTC")
dt <- 12 * 3600
juan <- subset(met, t0 - dt <= time & time <= t0 + dt)
par(mfrow = c(2, 1))
plot(juan, which = 5)
abline(v = t0)
plot(juan, which = 6)
abline(v = t0)