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.