[1] A thermomechanical ice flow model is formulated using the finite volume method. Separate submodels solve the full, two-dimensional momentum equations, the advective-diffusive heat equation, and evolution of the free surface. A unique aspect of the method is the use of a boundary-fitted, orthogonal, curvilinear coordinate system, which simplifies the implementation of boundary conditions, le...