The surface gradient method for the treatment of source terms in the shallow-water equations