Sylt elevation

## ============================================================================
## perspective plots of surface elevation in a mudflat
## data Syltsurf is part of package OceanView
## this is an example of the help file of Syltsurf
## implemented by Karline Soetaert
## ============================================================================
 library(OceanView)
 par(mfrow = c(3, 2), mar = c(0, 0, 1, 0))

 clim <- range(Syltsurf$elev, na.rm = TRUE)
 for (i in 1:5) 
   persp3D(z = -depth, colvar = Syltsurf$elev[,,i], x = x, y = y, 
     clim = clim, inttype = 2,  d = 2, 
     scale = FALSE, expand = 0.1, colkey = FALSE, shade = 0.5, 
       main = paste(format(Syltsurf$time[i], digits = 3), " hr"))
 par(mar = c(3, 3, 3, 3))
 colkey(clim = clim, clab = c("elevation", "m"))