Learn R Programming

ocedata (version 0.1.6)

papa: OWS Papa hydrographic record during 2010

Description

This dataset holds the temperature and salinity measured at Ocean Weather Station Papa (50N, 145W) in the year 2010. Time and vertical coordinate are in vectors papa$t and papa$z, while temperature, salinity and \(\sigma_\theta\) are in the matrices papa$temperature, papa$salinity and papa$sigmaTheta. This is an excerpt from a larger dataset that spans 2007-06-08 to 2012-10-24, and it has been trimmed to just those z values that are common to the larger dataset, namely -1, -10, -20, -45, -80, -100, -120, -150, and -200m. Time is sampled daily, from 2011-01-01 12:00:00 to 2011-12-31 12:00:00 UTC.

Usage

data(papa)

Arguments

References

The larger data set is described at http://www.pmel.noaa.gov/stnP/

Examples

Run this code
# NOT RUN {
library(oce)
data(papa, package="ocedata")

par(mfrow=c(3,1))
## Page 1: images
imagep(papa$t, papa$z, papa$temperature, col=oceColorsJet, filledContour=TRUE)
imagep(papa$t, papa$z, papa$salinity, col=oceColorsJet, filledContour=TRUE)
imagep(papa$t, papa$z, papa$sigmaTheta, col=oceColorsJet, filledContour=TRUE)
## Page 2: lineplots
xlim <- c(min(papa$t)-3*28*86400, max(papa$t)) # extra for legend
plot(papa$t, papa$temperature[,1], type='l', xlim=xlim,
     xlab="", ylab="Temperature", ylim=range(papa$temperature))
for (i in 2:dim(papa$temperature)[2])
    lines(papa$t, papa$temperature[,i], col=i, lwd=2)
legend("topleft", title="z [m]", lwd=2, legend=papa$z, col=1:length(papa$z), bg="white")
plot(papa$t, papa$salinity[,1], type='l', xlim=xlim,
     xlab="", ylab="Salinity", ylim=range(papa$salinity), lwd=2)
for (i in 2:dim(papa$salinity)[2])
    lines(papa$t, papa$salinity[,i], col=i, lwd=2)
legend("topleft", title="z [m]", lwd=2, legend=papa$z, col=1:length(papa$z), bg="white")
plot(papa$t, papa$sigmaTheta[,1], type='l', xlim=xlim,
     xlab="", ylab=expression(sigma[theta]), ylim=range(papa$sigmaTheta), lwd=2)
for (i in 2:dim(papa$sigmaTheta)[2])
    lines(papa$t, papa$sigmaTheta[,i], col=i, lwd=2)
legend("topleft", title="z [m]", lwd=2, legend=papa$z, col=1:length(papa$z), bg="white")
# }

Run the code above in your browser using DataLab