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

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

Added flag for handling surface versus surface0, fixed typos, changed code to resize surface arrays rather than deallocate/reallocate entire surface array structure