Supplementary Data and Code from: Sea-level contribution from Thwaites Glacier doubles in an observation-constrained ice-sheet model.

README The content of this repository is from: Sea-level contribution from Thwaites Glacier doubles in an observation-constrained ice-sheet model. 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. This can be done with perform(org,’Projection.GZ-Melt-Sensitivity’), after introducing this global variable 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-ISSM ‘model’ object. Usage notice: Before running the model, make sure that ISSM is properly installed on your machine. To do so, follow the instruction on the ISSM's Github page. After installation, you will have all routines needed to read and run the code. 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.