removing superflous variables and function parameters

Rotating indenter, nonmoving bulk and a spherical shell without
temperature or velocity constraints in between

Temperature is set at the left side and unset in the indenter footprint,
use to "pull" temperature from ther into indenter.

Also addition to elliptical dome surface: spherical shells of given
offset can be created.

Squashed commit of the following:

commit 01af6100baccbd389dca45e170c578cf642d73e7
Author: Matthias Schmiddunser <matthias.schmiddunser@uni-tuebingen.de>
Date:   Tue Nov 11 11:42:49 2014 +0100

    Adding option not to set temperature at bottom

    Use for restart after sstemp

commit 3f0d2ab899ffba26030e56353b2c5a12beb550a4
Author: Matthias Schmiddunser <matthias.schmiddunser@uni-tuebingen.de>
Date:   Fri Nov 7 16:18:42 2014 +0100

    Adding option domainfix and using old vex scaling

    domainfix will set velocity in entire domain 2 dz above indenter to 0
    returning to apply vex to osolve%z, since scaling the input parameters
    leads to different and incorrect shape of the indenter

commit 00a4fcd59cef94f7efc227f96297fa4d705bfce7
Author: Matthias Schmiddunser <matthias.schmiddunser@uni-tuebingen.de>
Date:   Thu Nov 6 16:40:12 2014 +0100

    scaling input parameters with vex

    moving vex out of actual code and scaling input parameters at the start

commit 44171a3be6b1bc42eb4d80af3453bd6b3fa2e3cf
Author: Matthias Schmiddunser <matthias.schmiddunser@uni-tuebingen.de>
Date:   Thu Nov 6 14:09:36 2014 +0100

    Changing Surface Calculation

    Trying a supposedly better spherical shell implementation

commit 1ee86ae001542ac5b37c7f41190b2d0cb913b8bd
Author: Matthias Schmiddunser <matthias.schmiddunser@uni-tuebingen.de>
Date:   Thu Nov 6 11:12:51 2014 +0100

    Fixed geotherm calculation

commit 8ef7e3c3fe5e6631e4629cce003e4242c408b5bd
Author: Matthias Schmiddunser <matthias.schmiddunser@uni-tuebingen.de>
Date:   Wed Nov 5 11:07:45 2014 +0100

    Added vex for geotherm calculation

commit 11e0f37f90b125014e254b3be9bab925b71590ef
Author: Matthias Schmiddunser <matthias.schmiddunser@uni-tuebingen.de>
Date:   Wed Nov 5 10:44:29 2014 +0200

    Fixed typos, removed pi

    Pi defined by use constants, set comparison variables in if clauses to
    double precision

commit 15bb1a8e8fa9944d694f0f7fee2354ef8b3109d7
Author: Matthias Schmiddunser <matthias.schmiddunser@uni-tuebingen.de>
Date:   Wed Nov 5 09:33:18 2014 +0100

    fixed minor typos

commit 9e9f1c7b9285f9497e4f1558179427520801eccf
Author: Matthias Schmiddunser <matthias.schmiddunser@uni-tuebingen.de>
Date:   Tue Nov 4 10:30:49 2014 +0100

    Adding rot_sub_tinit for Temperature Initialisation

    New set of boundary conditions similar to rot_subduction, but with
    kfixt=1 at the left boundary and kfixt=0 within indenter footprint.
    Strictflow=True by default.

commit aa2fc2c179ca5b378995879c097bd92ce332b2c7
Author: Matthias Schmiddunser <matthias.schmiddunser@uni-tuebingen.de>
Date:   Thu Oct 30 16:07:01 2014 +0100

    Added spherical shells description

commit 0c7c2532df3f0dd44c57e52a2c32a725b4cecc30
Author: Matthias Schmiddunser <matthias.schmiddunser@uni-tuebingen.de>
Date:   Thu Oct 30 16:02:47 2014 +0100

    Added spherical shells for indenter surface

    parallel surfaces of given distance to originally defined plane

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

This sets velocity and temp for nodes inside the indenter by setting
osolve%kfix... Option to set either velocity (rotation with vin),
temperature (uniform)  or both.

The input template has been updated accordingly.

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

Makefile.inc.esd.argand for compiling on Todd Ehler's work group's
(ESD's) Argand (Linux machine).

Also tidied up input file template input_schmiddu.txt.

5th temperature option added:
Radial Gradient from (0,y0,0) to indenter edge
Updated input template

using strictflow=1 (bc_pamam12) will set the boundary conditions
at the sides and back to (vback,0,0) and prevent vertical movement.

module version updated to version after compilation

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 variable calc_comp_velo, which is used to determined when compaction velocity should (or should not) be calculated

Fixed a few compaction problems: Deallocating certain osolve fields after compaction; Storing densityp=density in initialize_compaction; Several small bugs related to calculating the compaction velocity

Fixed problem with integer values being typed as doubles, set initialization value for cur_lsf, added absolute values for density comparisons, pruned old comments

Updated vermajor to 3 (v0.3), added options for restarting with compaction output, modified code to not read/write output for certain fields when they are not used (isostasy, compaction, nest)

Updated vermajor to 3 (v0.3), tweaked allocation of density string arrays, added options for restarting with compaction output, modified code to not read/write output for certain fields when they are not used (isostasy, compaction, nest)

Added option for viscosity/strain rate strengthening/hardening, added option for sstemp calculation on restart