Data and code for Recent observations of Thwaites Glacier, West Antarctica are consistent with high rates of loss in next 50 years

This repository contains the code and input files needed for running the MITgcm and ISSM experiments in this manuscript, and the relevant outputs from the experiments.

1. STREAMICE

the STREAMICE folder contains the following subfolders:

code: contains all MITgcm source files for calibration experiments that are not included in the MITgcm repository

code_proj: similar to code/ but for projection experiments

model_input: contains subdirectories containing binary inputs and parameter files for all experiments in the manuscript

model_output: output from each projection experiment is contained in a .mat file with the following arrays:

X, Y: x-and y-coordinates in WGS 84 / Antarctic Polar Stereographic (EPSG:3031)

BED: bed elevation (m)

BETA: \beta described in the manuscript (units: Pa^(1/2) (m/a)^{-1/6})

BGLEN: \overline{B} described in the manuscript (units: Pa^(1/2) a^{-1/6})

SURF: ice surface elevation (m)

THICK: ice thickness (m)

THICK0: initial ice thickness (m)

VX: x-velocity (m/a)

VY: y-velocity (m/a)

times:  time in years after 2004

checkpoint 68y of the MITgcm source code (mitgcm.org) was used for this study.

m1qn3 (https://who.rocq.inria.fr/Jean-Charles.Gilbert/modulopt/optimization-routines/m1qn3/m1qn3.html) was used for optimisation.

details of how to compile and run the forward and adjoint STREAMICE code can be found at https://mitgcm.readthedocs.io/en/latest/

-------------------------------------------

2. ISSM

runme.m contains all ISSM instructions to run the experiments.

output/ contains mat files for all experiments, with array descriptions as above.

ISSM source and documentation can be found at https://issm.jpl.nasa.gov/.