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