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

DOUAR with cloud mat modifications from November/December 2012. Needs testing, may not work properly.