Numerical methods for the primitive equations (PEs) of oceanic flow are presented in this paper. First, a two-dimensional Poisson equation with a suitable boundary condition is derived to solve the ...