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