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

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.

This has the retriangultion and flip-loop-check. It does not work with
remove_surf_pts=T.

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.

Conflicts:
	src/update_gitversion.sh

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.