A fully implicit parallel ocean model using MUMPS

The formulation, implementation and performance of a new fully implicit parallel model of the ocean circulation is presented. Within this model, steady states can be traced in one of the control parameters. In addition, transient flows can be computed using relatively (compared to traditional ocean models) large time steps such that long integration times can be reached. The discretized equations of the ocean model are solved by the Newton-Raphson technique and the emerging linear systems are solved by a (MPI) version of the MUltifrontal Massively Parallel Solver. The performance of the code on an SGI Origin 2000 platform is presented here using typical results for a sector ocean flow.