A particle-in-cell (PIC) numerical method developed for the study of shallow-water dynamics, when the moving fluid layer is laterally confined by the intersection of its top and bottom surfaces, is described. The effect of ambient rotation is included for application to geophysical fluids, particularly open-ocean buoyant vortices in which the underlying density interface outcrops to the surface around the rim of the vortex. Extensions to include the dynamical effect of a second moving layer (baroclinicity) and the presence of a lateral rigid boundary (sidewall) are also described. Although the method was developed for oceanographic investigations, applications to other fluid mechanics problems would be straightforward.