3D interpolation via gmt surface

See also my answer in Plot 2D profile using file having xyz and diffusion