Reverting changes and chosing edge splitting with random addition.

Extending new BC setup for rot_sub_init, editing the input file and
removing debugging output.

The transition zone is shifted so it lies above the indenter. This
ensures that, when using same numbers for h in BC definitions and in
surface parameters, the transition zone occurs in weaker overriding
plate.

x=1-face will be set to indenter velocities, if indenter is large
enough.
Strictflow will only set faces that don't cut the indenter.

Transition zone for BC only depends on distance from center.
define_bc_rot_sub_init.f90 still needs updating to the same setup.

New approach to the transition zone so that no conflicts at corners can occur.
Indenter peak can be freely chosen.

Allowing asymmetrical shapes, variable position of subduction and
reviving cylindrical option.

Age points are sent to external Python routine for interpolating onto
the current surface, read in again and a VTK surface file created.

Indenter can be rotated around the z-axis for angled influx.

Adjusting the code to call new functions.
Deleting unneeded features from age calculation modules.
Updating the example input file.

Received from Willi Kappler on 2015-02-24.

The list of all points to be removed is passed to remove_point and only
points that will not be removed in the future are accepted as nearest
neighbours that will take up the connections from the point to be
removed.
The remove_list is sorted for faster checking.