# NOT RUN {
library(oce)
data(secchi, package="ocedata")
mapPlot(coastlineWorld, longitudelim=c(-5, 20), latitudelim=c(50, 66),
grid=5, fill='gray', proj="lambert",
parameters=c(lat0=50, lat1=65))
col <- rev(oceColorsJet(100))[rescale(secchi$depth, xlow=0, xhigh=20, rlow=1, rhigh=100)]
mapPoints(secchi$longitude, secchi$latitude, pch=20, cex=1, col=col)
mapPolygon(coastlineWorld, col='gray')
# }
Run the code above in your browser using DataLab