On the convergence order of difference schemes for ocean dynamics equations

An algorithm for solving a finite-difference scheme that approximates the large-scale ocean dynamics equations with second order in spatial variables is proposed. Some results of numerical experiments performed to estimate the convergence order of the scheme are discussed. It is shown that the obtained numerical estimate is consistent with the previously obtained theoretical result. Moreover, it is shown that this estimate cannot be improved by convergence order.