We present a numerically accurate, computationally efficient, (km-scale) high-resolution model for gravitationally consistent relative sea level that, unlike contemporary state-of-the-art models, operates efficiently on an unstructured mesh. The model is useful for earth system modeling and space geodesy. A straightforward and computationally less burdensome coupling to a dynamical ice-sheet model, for example, allows a refined and realistic simulation of fast-flowing outlet glaciers.