Turbulent diffusion in shallow water equations

Add diffusion terms to the shallow water equations. One or more turbulence models will follow later. For now a rudimentary implementation is chosen which is only accurate on orthogonal grids. Possible non-orthogonal corrections will be added later on.