“…These considerations show that (21) specifies, at leading order, the background flow, provided that the stream function ψ(ϕ, θ, z) solves the vorticity equation (22) subject to the boundary condition (19) and to the last two constraints in (20); Equations (16)-(17) then determine the associated pressure field, taking also into account the boundary condition represented by the first constraint in (19). Let us note that if we ignore the z-dependence (and thus, implicitly, consider an inviscid setting in which the wind forcing plays no role), then (22) simplifies to the ocean gyre model derived recently in [24] as a shallow-water asymptotic solution of Euler's equation in rotating spherical coordinates (with the stipulation that θ stands in [24] for the polar angle, and not for the angle of latitude) and further investigated in [25][26][27][28][29] in the context of the Antarctic Circumpolar Current -the largest ocean current on Earth, flowing clockwise from west to east around Antarctica (see [30,32,37]) so that, due to the lack of any landmass connecting with Antarctica, it keeps the warm ocean waters from lower latitudes away from Antarctica and thus maintains the huge ice sheets encountered near the South Pole (see the discussion in [13,31]).…”