Software

Some publicly available projects are listed below. Additional ones can be found on my GitHub page.

h-NUMO

h-NUMO is a multilayer ocean model that utilizes shallow equations and is implemented within the Galerkin Numerical Modeling Environment (GNuME) framework. This framework uses an arbitrary polynomial basis function expansion and offers a choice for both continuous Galerkin and discontinuous Galerkin methods (Abdi and Giraldo, 2016). The work of Giraldo et al. (2013) is related to the current model through this framework. Previously, this framework was used to develop the Non-Hydrostatic Unified Model of the Ocean (Kopera et al., 2023).

Below, you can view animations of the double-gyre circulation using different polynomial approximations (Nop = 2, 4, 6), with the number of elements adjusted to maintain a grid resolution of 10 km across all simulations.

MPAS-Seaice: Sea ice component of E3SM, the climate model of the US Department of Energy

B-grid formulations of the sea ice dynamics have been dominant historically because they have matched the grid type used by ocean models. The divergence of the stress is the most challenging term to discretize in the sea-ice dynamics. For the computation of the divergence of the stress on a spherical mesh, the current formulation of MPAS-Seaice projects the cell coordinates onto planes tangent to the cell centers. We present an alternative approach where the cell coordinates are projected on common planes tangent to the vertex. This approach improves convergence and reduces errors on idealized test cases compared to the state-of-the-art method previous used in MPAS-Seaice model. More details can be found here.