5. Computational options

5.1. Barotropic mode solver

As part of the implicit solution of the barotropic mode, a two-dimensional elliptic equation for the surface pressure is solved. Two solver methods are available, both iterative: preconditioned conjugate gradient and a Chronopouos-Gear form of pcg requiring fewer global reductions. Convergence of the iterative solvers is governed by the two parameters convergenceCriterion and maxIterations as shown in the table below. The convergence criterion convergenceCriterion should be chosen small enough such that the pressure balance (printed as part of the model global diagnostics) agrees to 3-4 digits. The parameter maxIterations must be chosen large enough to allow the solver to converge (typically a few hundred), but small enough so that the code will terminate in a reasonable time if the solver is unable to converge.

Todo

Add link to namelist group &solvers for Barotropic solver namelist

Occasionally, when benchmarking the code, it is useful to fix the number of iterations to give a consistent iteration count between runs. In this case, convergenceCriterion is set to exactly zero and the solver will iterate maxIterations and continue with the simulation without terminating. This feature should only be used for benchmarking and not for actual ocean simulations.

The convergenceCheckFreq provides a means to improve performance by checking for convergence every convergenceCheckFreq iterations, thus eliminating an extra global sum on most iterations. Another means for improving performance is to supply a preconditioner to improve convergence. The preconditioner must be computed off-line and must be in the form of a nine point stencil operator. The preconditioner is then supplied in a file named preconditionerFile containing the nine operator weights.

5.2. Advection methods

Currently, advection of momentum is always done by leapfrog centered advection with periodic ‘mixing’ steps. For tracer advection, three options are available. The first is standard leapfrog centered advection; the second is a 3rd-order upwinding [Leonard] which, although not monotone, will improve monotonicity at a somewhat increased computational cost, and the third, available only in CESM1 POP2, is a one-dimensional Lax-Wendroff with one-dimensional flux limiters.

Todo

Add link to namelist group &advect for Advection namelist

5.3. Pressure gradient options

The pressure-averaging technique was explained in a previous section. Because it increases the timestep, it should always be enabled. The option to turn it off is provided only to permit comparisons with and without pressure-averaging or between POP and other codes that do not incorporate this technique.

The pressure gradient term includes a density factor which is assumed to be a constant reference density in Boussinesq models. The depth-dependent pressure effects on this density can be corrected for using simple depth-dependent factors.

Todo

Add link to namelist group &advect for Pressure averaging namelist