Adds computation of velocity boundary conditions from information on tectonic plates
The merge adds functionality for compute plate velocities to the TerraNeo module. These plate velocities can be used to set velocity boundary conditions for the Stokes part of a mantle convection model.