Parallel ocean general circulation modeling (Q687781)

We have developed a global ocean model on the massively parallel CM-2 Connection Machine based on the Bryan-Cox-Semtner ocean general circulation model. This paper discusses the Connection Machine implementation of the model and its performance, as well as major changes that were made in the algorithm and numerical methods. In order to increase the parallel efficiency, we reformulated the barotropic equations to solve for the surface-pressure field rather than the volume- transport streamfunction. We have also developed a new parallelizable preconditioning method for the solution of the elliptic equation for the surface-pressure field. In order to use a standard conjugate gradient algorithm, an approximate-factorization technique is used to split off the Coriolis terms and obtain a symmetric operator. We then apply a symmetrical local approximate-inverse operator as a preconditioning matrix, which is very effective in accelerating convergence, to a solution.