This function is essentially a combination of predict.surface and
plot.surface. It may not always give a great rendition but is easy to use
for checking the fitted surface. The default of extrap=F is designed to
discourage looking at the estimated surface outside the range of
the observations.
NOTE: that any Z covariates will b edropped and only the spatial part of
the model will be evaluated.