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

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)

Fixed bug where boolean variables might not be set to their default values when not given in the input file, added option for fixed cloud

DOES NOT COMPILE. This version with Janice's adjustment to calculate_density_profiles, which will now differentiate between initializing the strings and subsequently applying compaction. Note that currently the LSFs must be used for all calls to compaction; not yet compatible withs surface removal.