The Generalized Interpolation Material Point Method (GIMPM) for the Shallow Shelf Approximation (SSA) of ice flow with Damage
Developer: Alex Huth ([email protected] ; [email protected])
This repository contains the GIMPM-SSA-Damage model detailed in:
Huth, A., Duddu, R., & Smith, B. (2021a). A generalized interpolation material point method for shallow ice shelves. 1: Shallow shelf approximation and ice thickness evolution. Journal of Advances in Modeling Earth Systems, 13, e2020MS002277. https://doi.org/10.1029/2020MS002277
Huth, A., Duddu, R., & Smith, B. (2021b). A generalized interpolation material point method for shallow ice shelves. 2: Anisotropic nonlocal damage mechanics and rift propagation. Journal of Advances in Modeling Earth Systems, 13, e2020MS002292. https://doi.org/10.1029/2020MS002292
In addition to the GIMPM, this code also includes the standard Material Point Method (sMPM)
Damage models included (see Huth et al., 2021b for details):
- SSA creep damage (Huth et al., 2021b)
- SSA "zero-stress" damage (Sun et al., 2017)
- The SSA "zero-stress" damage model + a modification to include necking/mass balance effects (Bassis & Ma, 2015)
See README in PROG, which contains the main source code.
The examples from Huth et al., 2021a are found in test1d, test2d, and mismip/steady
The examples from Huth et al., 2021b are found in mismip/damage
Each directory contains a README with instructions for running the examples.
- An installation and knowledge of Elmer FEM and Elmer/Ice is required
- Tested with Elmer v 9.0