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