Added nest type for using nested model parameters. Reorganized things to have debug level 1 only write to screen; 2 to files.

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 to use an initial e2d value for the first plasticity calculation and output the dilatation rate

added storage/output of pre-isostasy velocity field; output of viscosity, is_plastic, elemental pressure