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

Reverted to using 'include mpif.h' as the use of mpi.mod requires the MPI module be compiled with the same compiler as the DOUAR source, which is not true on geodyncomp

Modified nonlinear iteration convergence criterion to calculate the maximum velocity difference between iterations over the maximum vector length from the previous iteration