Resistive stress / Nye crevasse formulation demo code for "Comprehensive Assessment of Stress Calculations for Crevasse Depths and Testing with Crevasse Penetration as Damage"

The jupyter notebook is set up to make plots of crevasse penetration with each resistive stress calculation from the paper for the Larsen B remnant and Pine Island Glacier ice shelves, as well as reproduce Figure 4e.

The script requires the numpy, xarray, and matplotlib python libaries. It has been confirmed to work in ghub's "earthaccess" conda environment, which can be used with a free account here: https://theghub.org/resources/4913.

Mirror of github repository: https://github.com/beynoldsj/crevCalcsDemo