read.section {oce} | R Documentation |
Read a file that contains a series of ctd
profiles that make up an oceanographic section.
read.section(file, section.id="", debug=FALSE, log.action)
file |
a file containing a set of CTD observations. At present, only the exchange BOT format is accepted (see Details). |
section.id |
optional string indicating the name for the section. If not provided, the section ID is determined by examination of the file header. |
debug |
logical. If TRUE , print some information that might be helpful during debugging. |
log.action |
if provided, the action item to be stored in the log. (Typically only provided for internal calls; the default that it provides is better for normal calls by a user.) |
Only exchange BOT comma-separated value format is
permitted at this time, but other formats may be added later. It
should also be noted that the parsing scheme was developed after
inspection of the A03 data set (see Examples). This may cause problems
if the format is not universal. For example, the header must name the
salinity column "CTDSAL
"; if not, salinity values will not be
read from the file.
An object of class
"section"
, which is a list
containing
data |
a list containing just one item,
station , which is a list of ctd objects, one per station. |
metadata |
a list containing:
header , the header from the data file;
section.id , a name for the section;
station.id , a vector of station IDs, one per station;
latitude , a vector of station latitudes, decimal and positive in northern hemisphere;
and
longitude , a vector of station latitudes, decimal and positive in eastern hemisphere.
|
processing.log |
A processing log, in the standard oce format. |
Dan Kelley
Several repository sites provide section data. An example
that is perhaps likely to exist for years is
http://cchdo.ucsd.edu, but a search on "WOCE bottle data"
should turn up other sites, if this one ceases to exist. Only the
so-called exchange BOT data format can be processed by
read.section() at this time.
make.section
can be used to bind together ctd
objects into a section.
library(oce) data(a03) Gulf.Stream <- subset(a03, indices=124:102) Gulf.Stream.gridded <- section.grid(Gulf.Stream, p=seq(0, 5000, 100)) data(coastline.world) plot(Gulf.Stream.gridded, coastline=coastline.world, map.xlim=c(-80,-60))