Journal cover Journal topic
The Cryosphere An interactive open-access journal of the European Geosciences Union
The Cryosphere, 2, 67-76, 2008
© Author(s) 2008. This work is distributed under
the Creative Commons Attribution 3.0 License.
13 Jun 2008
The ISMIP-HOM benchmark experiments performed using the Finite-Element code Elmer
O. Gagliardini1 and T. Zwinger2 1LGGE, CNRS, UJF-Grenoble I, BP 96, 38402 Saint-Martin d'Hères Cedex, France
2CSC-Scientific Computing Ltd., Keilaranta 14, P.O. Box 405, 02101 Espoo, Finland
Abstract. The aim of this paper is to describe in detail how the benchmark tests ISMIP-HOM (Ice Sheet Model Intercomparison Project-Higher-Order ice-sheet Model) have been performed using the open source finite element (FE) code Elmer ( The ISMIP-HOM setup consists of five diagnostic and one prognostic experiments, for both 2-D and 3-D geometries. For all the tests, the full-Stokes equations are solved. Some technical points concerning FE, such as mesh characteristics, stabilisation methods, numerical methods used to solve the linear system and parallel performance are discussed. For all these setups, the CPU time consumption in relation to the accuracy of the solution is analysed. Based on these findings, some general rules on optimising the computing time versus the accuracy of the results are deduced.

Citation: Gagliardini, O. and Zwinger, T.: The ISMIP-HOM benchmark experiments performed using the Finite-Element code Elmer, The Cryosphere, 2, 67-76,, 2008.
Publications Copernicus