To reduce uncertainties associated with sea level rise projections, an accurate representation of ice flow is paramount. Most ice sheet models rely on simplified versions of the underlying ice flow equations. Due to the high computational costs, ice sheet models based on the complete ice flow equations have been restricted to < 1000 years. Here, we present a new model setup that extends the applicability of such models by an order of magnitude, permitting simulations of 40 000 years.