Maps of υ_R and υ_θ. Top panels: comparison of the velocity fields υ_R and υ_θ in the form of residuals between the simulation and the theoretical predictions in Eqs. (12) and (34): (υ_R,sim − υ_R,th)/υ_θ,circ (top left panel) and (υ_θ,sim − υ_θ,th)/υ_θ,circ (top right panel), respectively. Bottom panels: similar to top panels but subtracting a circular Keplerian velocity profile (υ_θ,circ) from the υ_R,sim and υ_θ,sim in order to highlight the impact on the residuals of not accounting for the eccentric nature of the disc and of pressure corrections in the radial/azimuthal motion; (υ_R,sim − υ_R,circ)/υ_θ,circ (bottom left panel), where obviously υ_R,circ = 0, and (υ_θ,sim − υ_θ,circ)/υ_θ,circ (bottom right panel). In the bottom right panel it can be clearly seen that the simulation azimuthal velocity υ_θ,sim at large radii is sub-Keplerian, as υ_θ,circ does not account for the pressure support term.