Compute the interpolation of the gradient only in the surfaces

[Leguark]

1. In principle first you should compute a model normally.
2. then compute the vertices of the layer of interest
3. put the coordinates of the vertices into the function of compute gradients