Modified cloud injection/removal to not inject points above the free surface and to remove points above the free surface

Added first pass at coloring new sediment for stratigraphy; Fixed typo for sedimentation time stepping

First attempt to add a stratigraphy surface to allow proper definition of deposited sediment during sedimentation

With the option to combine progradation and aggradation. This version applies both using the original position of the sedimentation profile, and the total elapsed time. Sedimentation profiles will 'fill in' around one another.

Sedimentation allows for a) 2D-equivalent sinusoidal aggradation profile, and b) progradation profile with Gaussian curve in y superimposed on half-Gaussian in y.  Not yet compiled.

With debug output commented out, updates to written output in case of Cholesky factorization error, and the option to instruct douar to continue in the case of finding a bad value

Added changes (mschmiddunser 12.14) to allow output every nth time step. The also saves the most recent output, which is overwritten at each time step

Creating a new branch for recoding sediment compaction to use the volume of expelled fluid to calculate the compaction velocity

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)