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

• Main Authors: Book, Cameron, Hoffman, Matthew J., Kachuck, Samuel B., Hillebrand, Trevor R., Price, Stephen, Perego, Mauro, Bassis, Jeremy
• Format: Other/Unknown Material
• Language: English
• Published: Zenodo 2022
• Subjects: ice sheet modeling; glacier isostatic adjustment; Antarctica; sea-level rise; Thwaites Glacier; Antarc*; Ice Sheet
• Online Access: https://doi.org/10.5281/zenodo.6607706

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 ([e-mail]). 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 ...