Large-scale ice sheet modeling

We focus on improving the efficiency of ice sheet models coupled to the ocean. Constructing accurate models of the ice sheets on Greenland and Antarctica is crucial in order to reliably estimate sea level rise. Due to low efficiency of todays models it is currently not possible to run accurate enough simulations for long enough time spans.  We aim to run large scale simulations of Greenland and Antarctica at higher resolution and longer time spans than previously achieved with high order models. To do so, we relax a unnecessarily strict time step constraint and improve on the poor performance of linear solvers and preconditioners used. Within SESSI we aim to make sure that our numerical methods scale well and compare the performance of two ice sheet models based on the software Elmer and FEniCS. Furthermore, we hope to improve on load balancing when the model accuracy is varied within an ice sheet.