The computation of glacier movements leads to a system of nonlinear partial
differential equations. The existence and uniqueness of a weak solution is
established by using the calculus of variations. A discretization by the f
inite element method is done. The solution of the discrete problem is prove
d to be convergent to the exact solution. A first simple numerical algorith
m is proposed and its convergence numerically studied.