oce (version 0.9-13)

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 -200. Time is sampled daily, from 2011-01-01 12:00:00 to 2011-12-31 12:00:00 UTC.

Usage

data(papa)

source

Salinity is from http://www.pmel.noaa.gov/stnP/data/daily/s50n145w_dy.ascii and temperature is from http://www.pmel.noaa.gov/stnP/data/daily/t50n145w_dy.ascii. Although there is a data file for $\sigma_\theta$ on the NOAA website, here it was calculated with swSigmaTheta.

References

The larger data set is described at http://www.pmel.noaa.gov/stnP/, and related datasets can be found at http://www.pmel.noaa.gov/stnP/data.html in various formats.

Examples

library(oce)
data(papa)
par(mfrow=c(3,1))
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)
plot(papa$t, papa$temperature[,1], type='l',ylim=range(papa$temperature), lwd=2)
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',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',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")