R/papa.R

In ocedata: Oceanographic Data Sets for 'oce' Package

#' OWS Papa hydrographic record during 2010
#'
#' `papa` 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
#' \eqn{\sigma_\theta}{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.
#'
#' Salinity was inferred from
#' `https://www.pmel.noaa.gov/stnP/data/daily/s50n145w_dy.ascii` and
#' temperature from
#' `https://www.pmel.noaa.gov/stnP/data/daily/t50n145w_dy.ascii`, as downloaded
#' in early 2015; PMEL changed its website subsequent to that date, and as of
#' early 2016 there seems to be no way to link directly to the data; instead,
#' one must use a graphical interface provided at
#' `https://www.pmel.noaa.gov/ocs/data/disdel/`, which includes filling out a
#' form identifying the user and purpose.
#'
#' Although there is a data file for \eqn{\sigma_\theta}{sigma-theta} on the
#' NOAA website, here it was calculated with `swSigmaTheta()` in the
#' \CRANpkg{oce} package.
#'
#' The larger data set is described at `http://www.pmel.noaa.gov/stnP/`.
#'
#' @name papa
#'
#' @docType data
#'
#' @usage data(papa, package="ocedata")
#'
#' @examples
#' data(papa, package="ocedata")
#' 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")
NULL