Published at: Journal of Computational and Applied Mathematics, 116595

A compact scheme for the Munk boundary-layer equation in one dimension

Abstract

In this paper, we introduce a two-scale compact finite difference scheme for the equation
(MK-1D)

This equation serves as a model for the nonlinear barotropic equation (NB) governing oceanic flows.
(NB)

where
and
are the streamfunction and the wind stress tensor, respectively. This equation encodes the western boundary layer problem (Ghil et al. 2008) for the potential vorticity
, which corresponds to the sharp contrast between the gyres flow in the oceanic circulation at mid-latitude and the strong western boundary currents. Numerical results for Equation (MK-1D) show that, with this two-scale scheme, high order accuracy is preserved for
and
both in the boundary layer and in the central zone of the domain. The test cases are taken from Chekroun et al. 2020."

Read the full article