The computation of glacier movements leads to a system of nonlinear partial differentialequations. The existence and uniqueness of a weak solution is established by using the calculus ofvariations. A discretization by the finite element method is done. Thesolution of the discrete problem is proved to be convergent to the exactsolution. A first simple numerical algorithm is proposed and its convergence numericallystudied.