Submitted as: model description paper 12 Jan 2021

Submitted as: model description paper | 12 Jan 2021

Review status: this preprint is currently under review for the journal GMD.

icepack: a new glacier flow modeling package in Python, version 1.0

Daniel Shapero1, Jessica Badgeley2, Andrew Hoffmann2, and Ian Joughin1 Daniel Shapero et al.
  • 1Polar Science Center, Applied Physics Laboratory, University of Washington, Seattle, WA, USA
  • 2Department of Earth and Space Sciences, University of Washington, Seattle, WA, USA

Abstract. We introduce a new software package called icepack for modeling the flow of glaciers and ice sheets. Icepack is built on the finite element modeling library Firedrake, which uses the domain-specific language UFL for describing weak forms of partial differential equations. The diagnostic models in icepack are formulated through action principles that are specified in UFL. The components of each action functional can be substituted for different forms of the user’s choosing, which makes it easy to experiment with the model physics. The action functional itself can be used to define a solver convergence criterion that is independent of the mesh and requires little tuning on the part of the user. Icepack includes the 2D shallow ice and shallow stream models. We have also defined a 3D hybrid model based on spectral semi-discretization of the Blatter-Pattyn equations. Finally, icepack includes a Gauss-Newton solver for inverse problems that runs substantially faster than the BFGS method often used in the glaciological literature. The overall design philosophy of icepack is to be as usable as possible for a wide a swath of the glaciological community, including both experts and novices in computational science.

Daniel Shapero et al.

Status: final response (author comments only)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • RC1: 'Comment on gmd-2020-419', Douglas Brinkerhoff, 09 Feb 2021
    • AC1: 'Reply on RC1', Daniel Shapero, 07 Apr 2021
  • RC2: 'Comment on gmd-2020-419', Anonymous Referee #2, 19 Feb 2021
    • AC2: 'Reply on RC2', Daniel Shapero, 07 Apr 2021
    • AC4: 'Reply on RC2', Daniel Shapero, 07 Apr 2021
  • CC1: 'Community comment on gmd-2020-419', Ed Bueler, 19 Feb 2021
    • AC3: 'Reply on CC1', Daniel Shapero, 07 Apr 2021

Daniel Shapero et al.

Model code and software

icepack: glacier flow modeling with the finite element method in Python Daniel Shapero, Jessica Badgeley, and Andrew Hoffman

Daniel Shapero et al.


Total article views: 507 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
342 148 17 507 5 4
  • HTML: 342
  • PDF: 148
  • XML: 17
  • Total: 507
  • BibTeX: 5
  • EndNote: 4
Views and downloads (calculated since 12 Jan 2021)
Cumulative views and downloads (calculated since 12 Jan 2021)

Viewed (geographical distribution)

Total article views: 448 (including HTML, PDF, and XML) Thereof 447 with geography defined and 1 with unknown origin.
Country # Views %
  • 1
Latest update: 15 Apr 2021
Short summary
This paper describes a new software package called icepack for modeling the flow of ice sheets and glaciers. Glaciologists use tools like icepack to better understand how ice sheets flow, what role they have played in shaping earth's climate, and how much sea-level rise we can expect in the coming decades to centuries. Icepack includes several innovations to help researchers describe and solve interesting glaciological problems and to experiment with the underlying model physics.