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.

DOES NOT COMPILE!!! This is the first attempt at a complete compaction implementation. Needs debugging/testing!

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; fixed reference column weight calculation to work without internal LSFs/surfaces