shearStress returns the bed shear stress normalized by the water density. This normalisation is unfavourable, since the name of the function (and also the documentation) does not reflect it. Also the name
shearStress itself is confusing, since the bed shear stress is not a classical shear stress tensor. The name
bedShearStress is more clear, since it points to the context of shallow water equations.
Modify the function
shearStressto return the bed shear stress and not the bed shear stress normalized by the water density.
Return the positive shear stress