Add a throw for a timeloop error where the precision is exceeded
Renamed the get_dn functions properly
Remove write statements from diffusion
Check if highest grid point exceeds domain top to avoid senseless
simulations
Increased the precision of the model: WARNING THIS BREAKS ALL OLD
RESTART FILES FOR TIME