baroclinic instability in a pizza box

I describe the domain as a "pizza box" because it is a shallow box. Two simulations are shown: a low-resolution simulation with 41×41×11 grid points and a high-resolution simulation with 161×161×21 grid points.

In distance, it is much more shallow: 2000×2000×1. So we can think of the domain as a box 20,000 km on a side and 10 km deep, draped over the north pole ... but we keep the geometry Cartesian ... an f-plane.

Description of the plots

Upper left: horizontal cross-section at the surface

Lower left: vertical cross-section at mid y

Right: time trace

Hi-res: 161x161x21

The high-resolution simulation.

If you want to see it: the low-resolution simulation