A numerically stable coupled-mode formulation for acoustic propagation in range-dependent waveguides