Abstract.In this paper we present the GRISLI (Grenoble Ice Sheet and Land Ice) model in its newest revision (version 2.0). Whilst GRISLI is applicable to any given geometry, we focus here on the Antarctic ice sheet because it highlights the importance of grounding line dynamics. Important improvements have been implemented since its original version (Ritz et al., 2001) including notably an explicit flux computation at the grounding line based on the analytical formulations of Schoof (2007) 5 and Tsai et al. (2015) and a basal hydrology model. A calibration of the mechanical parameters of the model based on an ensemble of 150 members sampled with a Latin Hypercube method is used. The ensemble members performance is assessed relative to the deviation from present-day observed Antarctic ice thickness. The model being designed for multi-millenial longterm integrations, we also present glacial-interglacial ice sheet changes throughout the last 400 kyr using the best ensemble members. To achieve this goal, we construct a simple climatic perturbation of present-day climate forcing fields based on 10 two climate proxies, both atmospheric and oceanic. The model is able to reproduce expected grounding line advances during glacials and subsequent retreats during terminations with reasonable glacial-interglacial ice volume changes.