The 3-dimensional finite difference thermodynamic coupled model on Antarctic ice sheet,GLIMMER model,is described. An ide-alized ice sheet numerical test was conducted under the EISMINT-1 benchmark,and the characteristic curves of ice sheets under steady state were obtained. Based on this,this model was simplified from a 3-dimensional one to 2-dimensional one. Improvement of the dif-ference method and coordinate system was proposed. Evolution of the 2-dimensional ice flow was simulated under coupled temperature field conditions. The results showed that the characteristic curves deriving from the conservation of the mass,momentum and energy agree with the results of ice sheet profile simulated with GLIMMER model and with the theoretical results. The application prospect of the simplified 2-dimensional ice flow model to simulate the relation of age-depth-accumulation in Dome A region was discussed. The 3-dimensional finite difference thermodynamic coupled model on Antarctic ice sheet, GLIMMER model, is described. An ide-alized ice sheet numerical test was conducted under the EISMINT-1 benchmark, and the characteristic curves of ice sheets under steady state were obtained. Based on this, this model was simplified from a 3-dimensional one to 2-dimensional one. Improvement of the dif-ference method and coordinate system was proposed. Evolution of the 2-dimensional ice flow was simulated under coupled temperature field conditions. results showed that the characteristic curves deriving from the conservation of the mass, momentum and energy agree with the results of ice sheet profile simulated with GLIMMER model and with the theoretical results. The application prospect of the simplified 2-dimensional ice flow model to simulate the relation of age-depth-accumulation in Dome A region was discussed.