oce (version 0.8-10)

section-class: Class to store hydrographic section data

Description

Class to store hydrographic section data, with standard slots metadata, data and processingLog.

Methods

The examples given below use s to represent the name of a section.

Extracting values{A list of stations is retrieved by s[["section"]]. Individual stations are retrieved by providing a station number as a second argument in the index, e.g. the first station is s[["section", 1]] (which is a ctd-class object). Aggregated values of the quantities measured at each level of the CTD profiles contained within the section may be accessed as e.g. section[["salinity"]]. This works for any quantity whose name is present in the constituent profiles.

Since it is often useful to pair such quantities with locations, section[["latitude"]] and section[["longitude"]] return vectors with values repeated for each level in each CTD (see the pairs() call in the example section). If just one latitude or longitude is desired per station, e.g. section[["latitude", "bystation"]] may be used. Station-by-station values of dynamic height are provided by e.g. section[["dynamic height"]].}

Assigning values{(This feature has yet to be added.)} Overview of contents{The show method (e.g. show(section)) displays information about the object.}

See Also

Sections can be read with read.section or created with read.section or created from CTD objects by using makeSection or by combining an existing section with a CTD object with the + operator. Sections may be sorted with sectionSort, smoothed with sectionSmooth, and gridded with sectionGrid. Gridded sections may be plotted with plot.section. Statistical summaries are provided by summary.section, while overviews are provided by show.section. Two section datasets are provided in the package: a03, containing data along WOCE line A03, and section, containing data in Halifax Harbour.

Examples

library(oce)
data(a03)
plot(a03[['station', 1]])
pairs(cbind(z=-a03[["pressure"]],T=a03[["temperature"]],S=a03[["salinity"]]))
## T profiles for first few stations in a03, at common scale
par(mfrow=c(3,3))
Tlim <- range(a03[["temperature"]])
ylim <- rev(range(a03[["pressure"]]))
for (stn in a03[["station",1:9]])
    plotProfile(stn, xtype='temperature', ylim=ylim, Tlim=Tlim)