Data from: Seawater intrusions control the retreat of Thwaites Glacier, West Antarctica

README

The content of this repository is from: Seawater intrusions control the retreat of Thwaites Glacier, West Antarctica (under review).

Description of the data and file structure

File description:

runme.m: MATLAB script to reproduce results presented in the manuscript. Notably, the sensitivity simulations discussed in the paper are initialized and run in step 16.

Thwaites.mat: All data associated with the ice sheet model initial state is included in Thwaites.mat, after spin-up (step 11 in readme.m). The file contains a MATLAB ‘model’ object. 

Usage notice:

In MATLAB, the model can be loaded and displayed by running load(‘Thwaites.mat’), which will load in the model variable ‘md’. Of particular interest will be the following data contained in md: md.mesh (mesh information), md.geometry (initial ice sheet geometry, ice shelf geometry, and bed topography), and md.mask (ice mask and grounded ice mask). Note that all fields are defined on the mesh nodes, and one can plot a given field in MATLAB using the ISSM tool ‘plotmodel’.

Ice sheet modeling results and initial states are compatible with the open source, NASA funded Ice-sheet and Sea-level System Model (ISSM, Larour et al. 2012), which is freely available for download on Github. In addition, the data files provided in the publication are available as *.mat files, which are compatible with MATLAB but can be accessed using most scripting languages.