Solution (source code)

= Solution

<Potential-vorticity conservation> makes the evolving PV distribution a set of parcel-carried labels. In the ideal model it can be rearranged into filaments and complicated geometry, but it is not created or destroyed by advection. Apparent homogenization, such as that assumed for the shallow-water experiment, describes a coarse-grained mixed distribution rather than pointwise destruction of materially conserved PV.

Within a balanced regime, <potential-vorticity inversion> turns those labels and the boundary data into a <velocity> field. This makes understanding a flow a coupled geometric problem: determine how PV is rearranged, invert the resulting distribution, and use the recovered <velocity> to find the next rearrangement. Inversion is nonlocal, so a compact PV anomaly influences fluid well outside its support, as the spherical example explicitly demonstrates.

The <quasi-geostrophic potential-vorticity equation> has only one time derivative:
$$
\partial_tQ+J(\psi,Q)=0,\qquad
J(\psi,Q)=\psi_xQ_y-\psi_yQ_x,\qquad\nabla_*^2\psi=Q-f.
$$
\b[An initial PV field and the inversion data determine the balanced evolution; an independent initial acceleration is not needed.] <Pressure>, <buoyancy> and the leading <velocity> are diagnosed from that field at each instant. This reduced initial-data requirement reflects the removal of freely oscillating inertia-gravity waves, which in the full fluid equations require additional initial information. The principles are consequently powerful for slow large-scale flow, but do not recover every unbalanced motion or remove the need to state <boundary conditions> and the balance approximation.