500-yr Projections of Thwaites Glacier, Antarctica, with MALI, including glacial isostatic adjustment

This archive contains model code, results, and analysis scripts for
reproducing the material presented in the manuscript "Stabilizing effect of
bedrock uplift on retreat of Thwaites Glacier, Antarctica, at centennial
timescales" by Cameron Book, et al.  Questions should be directed to Matt
Hoffman (mhoffman@lanl.gov).

This archive contains the following directories:

|-- MALI_code
|-- PIGL_control
|-- analysis
|-- control
|-- run_setup
|-- N1
|-- N2
|-- N3
|-- N4
`-- PIGL_N3

'MALI_code' is a snapshot of the MALI repository used for these simulations,
commit 454e0fc8bf384bee1c1560d6c3eaa5fae43bfdda.
This commit is present in an older MALI repository that is no longer
maintained, at https://github.com/MPAS-Dev/MPAS-Model
MALI is currently maintained on Github at https://github.com/MALI-Dev/E3SM
Building MALI requires the Albany multiphysics library, which is available
at https://github.com/sandialabs/Albany.  The simulations presented use Albany
master from March 5, 2021.

'analysis' contains the scripts used to process the model output and produce
the figures and results presented in the manuscript.  Filepaths will have to
be adjusted to your local layout.

'run_setup' is a directory of files and scripts necessary to reproduce the
model simulations presented.  The GIA model giapy is in the file
'giascript.py'.  giapy can also be found on Github at https://github.com/skachuck/giapy

'control' is the control run with the GIA model disabled.  It corresponds to
the run labeled CTRL in the manuscript.

'PIGL_control' is the control run using the high melt forcing.  It corresponds
to the run labeled HM-CTRL in the manuscript.

The five run directories included here (N1-N4, PIGL_N3) are the standard
ensemble described in the manuscript.  (There is a separate archive for the
runs briefly mentioned that exclude the elastic response of the lithosphere.)
The individual runs have the following correspondence to the manuscript:
N1=TYP
N2=BEST2
N3=VLV-THIN
N4=VLV
PIGL_N3=HM-VLV-THIN
Within each run directory, are the following model output files:
globalStats.nc: MALI global, scalar time-series
iceload_all.nc: Thwaites Glacier ice load on the GIA grid
output_*.nc: MALI spatial output fields, separated by century
uplift_GIA_all.nc: GIA output on GIA grid