Learn R Programming

oce (version 0.8-3)

geosecs235: GEOSECS station 235 data

Description

GEOSECS station 235 data

Usage

data(geosecs235)

Arguments

source

Data from http://iridl.ldeo.columbia.edu/SOURCES/.GEOSECS/ISTA/237/VALUE/data.cdf

References

D. E. Kelley and Van Scoy, K. A., 1999. A basin-wide estimate of vertical mixing in the upper pycnocline: spreading of bomb tritium in the North Pacific Ocean. Journal of Physical Oceanography, 29: 1759--1771.

Examples

Run this code
library(oce)
data(geosecs235)
## reproduce part of Figure 10 of Kelley and Van Scoy (1999)
plotProfile(geosecs235, xtype="tritium", type='p', ylim=c(600, 0))
trit <- tritium(geosecs235)
p <- pressure(geosecs235)
m <- nls(trit~a*exp(-((p-p0)/D)^2)+a*exp(-((p+p0)/D)^2), start=list(a=10,p0=200,D=100))
pp <- seq(600, 0, -10)
lines(predict(m, newdata=list(p=pp)), pp)

Run the code above in your browser using DataLab