interpolation is a two-stage process: first combining the well-data
for each 0.2 block as a barycentric approach on a course grid with
0.2 grid spacing, using a weighted combination of the data from the
nearest three available wells or boundary points. In the case of the
sediment thickness, the thickness was assumed to be zero on the
boundary and the barycentric interpolation included at most two
boundary values. In the second stage, the course grid was refined to
0.01 resolution using a Delaunay triangulation.