We consider the paper by Brugnot and Pochat (1981), which describes a one-dimensional model applied to a snow avalanche. The main advance made here is the introduction of the second dimension in the runout zone. Indeed, in the channelled course, we still use the one-dimensional model, but, when the avalanche spreads before stopping, we apply a (x, y) grid on the ground and six equations have to be solved: (1) for the avalanche body, one equation for continuity and two equations for momentum conservation, and (2) at the front, one equation for continuity and two equations for momentum conservation. We suppose the front to be a mobile jump, with longitudinal velocity varying more rapidly than transverse velocity.

We solve these equations by a finite difference method. This involves many topological problems, due to the actual position of the front, which is defined by its intersection with the reference grid (SI, YJ). In the near future our two directions of research will be testing the code on actual avalanches and improving it by trying to make it cheaper without impairing its accuracy.