Abstract: ABSTRACT We consider near-inertial waves continuously excited by a localized source and their subsequent radiation and evolution on a two-dimensional $\ensuremath{\beta} $-plane. Numerical simulations are used to quantify the wave propagation and the energy flux in a realistically stratified ocean basin. We focus on the dynamics near and poleward of the inertial latitude where the local value of the Coriolis parameter $f$ matches the forcing frequency $\sigma $, contrasting the behaviour of waves under the traditional approximation (TA), where only the component of the Earths rotation aligned with gravity is retained in the dynamics, with that obtained under the non-traditional approach (non-TA) in which the horizontal component of rotation is retained. Under the TA, assuming inviscid linear wave propagation in the WKB limit, all energy radiated from the source eventually propagates toward the equator, with the initially poleward propagation being internally reflected at the inertial latitude. Under the non-TA however, these waves propagate sub-inertially beyond their inertial latitude, exhibiting multiple reflections between internal turning points that lie poleward of the inertial latitude and the bottom. The numerical experiments complement and extend existing theory by relaxing the linearity and WKB approximations, and by illustrating the time development of the steadily forced flow and the spatial patterns of energy flux and flux divergence. The flux divergence of the flow at both the forcing frequency and its first harmonic reveal the spatial patterns of nonlinear energy transfer and highlight the importance of nonlinearity in the vicinity of near-critical bottom reflection at the inertial latitude of the forced waves.