pPOM: A nested, scalable, parallel and Fortran 90 implementation of the Princeton Ocean Model

In this work we describe the development of a parallel implementation of the Princeton Ocean Model (POM2k) with a nested-domain feature. Parallelization has been handled using the Run-Time System Library (RSL) and the Fortran Loop Index Converter (FLIC), avoiding a direct use of the MPI library. Modularity and flexibility have been added through advanced Fortran 90 features, such as modules, dynamic memory allocation, pointers and recursion. The “seamount problem”, either in a nested and non-nested configuration, is used as a test bed for showing results scalability.