We use a numerical flow model to simulate the behaviour of the Djankuat Glacier, a WGMS reference glacier situated in the North Caucasus (Republic of Kabardino-Balkaria, Russian Federation), in response to past, present and future climate conditions (1752–2100 CE). In particular, we adapt a more sophisticated and physically based debris model, which has not been previously applied in time-dependent numerical flow line models, to look at the impact of a debris cover on the glacier’s evolution.