library(oce)
data(ctd)
sigthe <- sw.sigma.theta(ctd$data$salinity,
ctd$data$temperature,
ctd$data$pressure);
ctd.new <- ctd.add.column(ctd,
sigthe, "sigmatheta",
"sigma-theta", "density, sigma-theta",
"kg/m^3");
Run the code above in your browser using DataLab