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.