We present an open-source numerical tool to simulate the free-surface evolution of gravity-driven flows (e.g. glaciers) constrained by bed topography. No ad hoc post-processing is required to enforce positive ice thickness and mass conservation. We utilise finite elements, define benchmark tests, and showcase glaciological examples. In addition, we provide a thorough analysis of the applicability and robustness of different spatial stabilisation and time discretisation methods.