library(oce)
data(a03)
plot(a03[['station', 1]])
pairs(cbind(z=-a03[["pressure"]],T=a03[["temperature"]],S=a03[["salinity"]]))
## T profiles for first few stations in a03, at common scale
par(mfrow=c(3,3))
Tlim <- range(a03[["temperature"]])
ylim <- rev(range(a03[["pressure"]]))
for (stn in a03[["station",1:9]])
plotProfile(stn, xtype='temperature', ylim=ylim, Tlim=Tlim)
Run the code above in your browser using DataLab