surf_sim function

Simulate surface data based on a spatial model