### Abstract

Following the work of Matthies (2001), it is shown how Gevrey (exponential) regularity of the solution and a classical method used together to prove an exponentially-accurate approximation result for a singular perturbation problem with a small parameter. The model considered is the viscous primitive equations of the ocean, although the method is applicable more generally.