Plot a summary diagram for argo data.

# S4 method for argo
  which = 1,
  coastline = c("best", "coastlineWorld", "coastlineWorldMedium", "coastlineWorldFine",
  cex = 1,
  pch = 1,
  type = "p",
  col = 1,
  fill = FALSE,
  projection = NULL,
  mgp = getOption("oceMgp"),
  mar = c(mgp[1] + 1.5, mgp[1] + 1.5, 1.5, 1.5),
  debug = getOption("oceDebug"),



an argo object.


list of desired plot types, one of the following. Note that oce.pmatch() is used to try to complete partial character matches, and that an error will occur if the match is not complete (e.g. "salinity" matches to both "salinity ts" and "salinity profile".).

  • which=1, which="trajectory" or which="map" gives a plot of the argo trajectory, with the coastline, if one is provided.

  • which=2 or "salinity ts" gives a time series of salinity at the indicated level(s)

  • which=3 or "temperature ts" gives a time series of temperature at the indicated level(s)

  • which=4 or "TS" gives a TS diagram at the indicated level(s)

  • which=5 or "salinity profile" gives a salinity profile

  • which=6 or "temperature profile" gives a temperature profile

  • which=7 or "sigma0 profile" gives a sigma0 profile

  • which=8 or "spice profile" gives a spiciness profile, referenced to the surface. (This is the same as using which=9.)

  • which=9 or "spiciness0 profile" gives a profile of spiciness referenced to a pressure of 0 dbar, i.e. the surface. (This is the same as using which=8.)

  • which=10 or "spiciness1 profile" gives a profile of spiciness referenced to a pressure of 1000 dbar.

  • which=11 or "spiciness2 profile" gives a profile of spiciness referenced to a pressure of 2000 dbar.


depth pseudo-level to plot, for which=2 and higher. May be an integer, in which case it refers to an index of depth (1 being the top) or it may be the string "all" which means to plot all data.


character string giving the coastline to be used in an Argo-location map, or "best" to pick the one with highest resolution, or "none" to avoid drawing the coastline.


size of plotting symbols to be used if type="p".


type of plotting symbols to be used if type="p".


plot type, either "l" or "p".


optional list of colors for plotting.


either a logical, indicating whether to fill the land with light-gray, or a color name. Owing to problems with some projections, the default is not to fill.


character value indicating the projection to be used in trajectory maps. If this is NULL, no projection is used, although the plot aspect ratio will be set to yield zero shape distortion at the mean float latitude. If projection="automatic", then one of two projections is used: stereopolar (i.e. "+proj=stere +lon_0=X" where X is the mean longitude), or Mercator (i.e. "+proj=merc") otherwise. Otherwise, projection must be a character string specifying a projection in the notation used by oceProject() and mapPlot().


a 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.


value to be used with par("mar").


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


debugging flag.


optional arguments passed to plotting functions.




Dan Kelley


tc <- cut(argo[["time"]], "year")
# Example 1: plot map, which reveals float trajectory.
plot(argo, pch = as.integer(tc))
year <- substr(levels(tc), 1, 4)
contour(topoWorld[["longitude"]], topoWorld[["latitude"]],
    add = TRUE
legend("bottomleft", pch = seq_along(year), legend = year, bg = "white", cex = 3 / 4)

# Example 2: plot map, TS, T(z) and S(z). Note the use
# of handleFlags(), to skip over questionable data.
plot(handleFlags(argo), which = c(1, 4, 6, 5))