Changed the divergence term to div((H+h)*u) in the continuity equation.
Updated the model equations in the manual.