deleting declared but unused variables, removing dummy arguments,
changing variables to double precision.
Not touching: uninitialized and more complicated warnings.

To eliminate warnings, declared but unneeded variables were removed from
define_bc and define_bc_rot_subduction

New boundary conditions:
 - A cylindrical or ellipsoidical indenter that is rotating downward set
   in horizontal flow field
 - Temperature conditions at the base can be set to a gradient or
   different for indenter (uniform or gradient)

New Surface:
 - Top of a cylinder with axis along y or of an ellipsoid rising above plane

Input files edited to accomodate new BCs and new surface

Example input file describing new options provided

Makefile for Taito updated (August 2014)

Added option for a pre-run isostatic compensation stage with a limit on the rate of compensation based on the geometry and a defined max velocity