glaciome1D: python code for modeling ice mélange flow, 2024

glaciome1D is a quasi-one-dimensional continuum model for modeling the flow of dense packs of icebergs (ice mélange) found in some fjords. In many respects the model is similar to one-dimensional models of ice streams and ice shelves, except that it uses the nonlocal granular fluidity rheology of Henann and Kamrin (2013). The model was created with the intention of developing coupled glacier-ocean-melange models. This is reflected in the modeling framework, which mimics that used for ice streams and ice shelves (Schoof, 2007).
Using the model involves creating an instance of the glaciome class, which contains information on the glacier velocity, viscosity (granular fluidity), and geometry as well as model parameters and various external forcings. The glaciome class includes several basic and easy to use functions, such as: self.diagnostic(), self.prognostic(), self.steadystate(), self.save().
The model physics and numerics are described in detail in Amundson et al. (in press).
This data set includes a single python module that includes all of the functions for setting up and running the model, an example script that runs the model, and a conda environment list that contains python modules that the code has been tested on.