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.

Check_delaunay got stuck in an infine loop, always flipping the same
edge. This checks the flipped edges compared to the last time and aborts
loop if it is stuck.

triangle array was redundant

If node has too many neighbours in check_delaunay, a new subroutine
retrinagulate_surface tries to redo the triangulation before aborting.

Using the same algorithm as in move_cloud.f90

Only longest edge in one triangle is split.
If two edges of the same length exist in one triangle, the one with the
lower edge index is split.
If an edge is split no other edge in the adjacent triangles may be
split.

find_leaf crashes if particle lies outside domain too far.

For no advection through z=0-plane

Temperature is initialized by assuming all surfaces are more or less
flat and sp01 defines the base elevation. Temperature is then calculated
as a steady state solution with T=0 at bottom and T=1 at ztemp.

Thermochron Ages are calculated in a self-scheduling manner. Write
output moved to subroutine. Output only in *.vtk

Self-Scheduling Model for multiple processors.

Upon restarting, the restart file is only read by Master process and
then broadcasted to other processes.

Also doing some source code formatting.

Tracking output can be written in regular intervals instead of every
timestep.
Output directory for post-track can be specified in input parameter
file.
Warning added for tracking cloud to be activated after spinup.

Post processor for tracking output files. It produces vtk files for
regular spced output timesteps of either all or only surfaced tracks.
It works with regularly spaced input files.
Output variables are: Speed, vertical speed, temperature, pressure,
time, step number and particle/track ID. Variables and positions are
written either in DOUAR units or natural units.