data(section)
ctd
objects, in
section$data$station[[i]]
, with i
in the range from 1 to
7.8
was evidently
typed as 3
).library(oce) stn08 <- ctdTrim(read.ctd("BED0308.CNV", station="8")) stn09 <- ctdTrim(read.ctd("BED0309.CNV", station="9")) stn10 <- ctdTrim(read.ctd("BED0310.CNV", station="10")) stn01 <- ctdTrim(read.ctd("BED0301.CNV", station="1")) stn11 <- ctdTrim(read.ctd("BED0311.CNV", station="11")) stn12 <- ctdTrim(read.ctd("BED0312.CNV", station="12")) stn12 <- oceEdit(stn12, "latitude", 44 + 39.894 / 60, reason="typo in lat minutes", person="Dan Kelley") stn13 <- ctdTrim(read.ctd("BED0313.CNV", station="13")) stn15 <- ctdTrim(read.ctd("BED0315.CNV", station="15")) stn15 <- oceEdit(stn15, "latitude", stn15$metadata$latitude + 1, reason="typo in lat degrees", person="Dan Kelley") section <- makeSection(stn08,stn09,stn10,stn01,stn11,stn12,stn13,stn15) section <- oceEdit(section, "sectionId", "Halifax-Harbour-2003-oct", person="Dan Kelley") section <- sectionGrid(section, p=seq(0,60,2)) save(section, file="section.rda") The stations cover a 14 km region running from station 308, near the Sackville River, seaward to station 315, at the entrance to the general Harbour, offshore of the aptly-named Point Pleasant Park.
summary.section
,
and a plot can be created with plot.section
.data(section)
data(coastlineHalifax)
plot(section, coastline=coastlineHalifax, grid=TRUE)
Run the code above in your browser using DataLab