ctd
profiles that make up an oceanographic section.read.section(file, sectionId="", flags,
ship="", scientist="", institute="",
missingValue=-999,
debug=getOption("oceDebug"),
processingLog)
TRUE
, print some information that
might be helpful during debugging.class
"section"
, which is a list
containingstation
, which is a list of ctd
objects, one per station.header
, the header from the data file;
sectionId
, a name for the section;
stationId
, 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.oce
format.CTDSAL
"; if not, salinity values will not be
read from the file."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.section-class
explains the
structure of section objects, and also outlines the other functions dealing
with them.library(oce)
data(a03)
GS <- subset.oce(a03, indices=124:102)
GSg <- sectionGrid(GS, p=seq(0, 5000, 100))
data(coastlineWorld)
plot(GSg, coastline=coastlineWorld, map.xlim=c(-80,-60))
# Plot stations individually (good for quality-control)
pdf("a03_stations.pdf")
for (s in seq_along(na03@data$station)) {
plot(na03@data$station[[s]])
}
dev.off()
Run the code above in your browser using DataLab