The ISMIP-HOM benchmark experiments performed using the Finite-Element code Elmer

• Main Authors: Gagliardini, O., Zwinger, T.
• Other Authors: Laboratoire de glaciologie et géophysique de l'environnement (LGGE), Observatoire des Sciences de l'Univers de Grenoble (OSUG), Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP )-Institut national des sciences de l'Univers (INSU - CNRS)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Savoie Mont Blanc (USMB Université de Savoie Université de Chambéry )-Centre National de la Recherche Scientifique (CNRS)-Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP )-Institut national des sciences de l'Univers (INSU - CNRS)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Savoie Mont Blanc (USMB Université de Savoie Université de Chambéry )-Centre National de la Recherche Scientifique (CNRS)-Centre National de la Recherche Scientifique (CNRS), Scientific Computing Ltd (CSC)
• Format: Article in Journal/Newspaper
• Language: English
• Published: HAL CCSD 2008
• Subjects: [SDU.ENVI]Sciences of the Universe [physics]/Continental interfaces; environment; [SDU.OCEAN]Sciences of the Universe [physics]/Ocean; Atmosphere; [SDU.STU]Sciences of the Universe [physics]/Earth Sciences; Ice Sheet; The Cryosphere
• Online Access: https://hal.science/hal-00331157; https://hal.science/hal-00331157/document; https://hal.science/hal-00331157/file/tc-2-67-2008.pdf

International audience 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 (http://www.csc.fi/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.