See Reproducibility repository On Affordable High-Order Entropy-Conservative/Stable and Well-Balanced Methods for Nonconservative Hyperbolic Systems
to generate

This can be generalized in such a way that the surface level is chosen, as well as longitude and latitude, with a user function for a given quantity. In the vortex shedding test cases for example it is interesting to compute the relative vorticity and zonal wind perturbations.
As far as I know, at the moment there's no way to do such plots if not with PyPlot. Ideally I would like to use CairoMakie (https://docs.makie.org/dev/reference/plots/tricontourf) but I cannot reproduce the same results.
See Reproducibility repository On Affordable High-Order Entropy-Conservative/Stable and Well-Balanced Methods for Nonconservative Hyperbolic Systems
to generate

This can be generalized in such a way that the surface level is chosen, as well as longitude and latitude, with a user function for a given quantity. In the vortex shedding test cases for example it is interesting to compute the relative vorticity and zonal wind perturbations.
As far as I know, at the moment there's no way to do such plots if not with PyPlot. Ideally I would like to use CairoMakie (https://docs.makie.org/dev/reference/plots/tricontourf) but I cannot reproduce the same results.