library(oce)
data(levitus)
par(mfrow=c(2,1))
imagep(levitus$longitude, levitus$latitude, levitus$SST, col=oceColorsJet, zlim=c(-2, 30))
imagep(levitus$longitude, levitus$latitude, levitus$SSS, col=oceColorsJet, zlim=c(20, 40))
Run the code above in your browser using DataLab