Wide-angle mode parabolic equations for the modelling of horizontal refraction in underwater acoustics and their numerical solution on unbounded domains