A hybrid characteristics-finite difference method for unsteady flow in pipe networks