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/.本仓库包含了复现本文中MITgcm与ISSM试验所需的代码与输入文件，以及相关试验的输出结果。 1. STREAMICE STREAMICE文件夹包含以下子目录： code文件夹：存放未纳入MITgcm官方仓库的校准试验所用全部MITgcm源代码文件。 code_proj文件夹：与code文件夹功能类似，仅用于投影试验。 model_input文件夹：包含本文所有试验所用二进制输入文件与参数文件的子目录。 model_output文件夹：每个投影试验的输出结果均存储于.mat格式文件中，包含以下数组： X、Y：采用WGS 84 / 南极极射投影（EPSG:3031）的x、y坐标。 BED：冰床高程（单位：米）。 BETA：本文提及的β参数（单位：Pa^(1/2) (m/a)^{-1/6}）。 BGLEN：本文提及的$overline{B}$参数（单位：Pa^(1/2) a^{-1/6}）。 SURF：冰面高程（单位：米）。 THICK：冰厚（单位：米）。 THICK0：初始冰厚（单位：米）。 VX：x方向流速（单位：米/年）。 VY：y方向流速（单位：米/年）。 times：2004年起的时间（单位：年）。 本研究采用MITgcm源代码的checkpoint 68y版本（来源：mitgcm.org）。 本研究采用m1qn3优化工具（https://who.rocq.inria.fr/Jean-Charles.Gilbert/modulopt/optimization-routines/m1qn3/m1qn3.html）完成优化运算。 关于编译并运行正向与伴随STREAMICE代码的详细步骤，请参考https://mitgcm.readthedocs.io/en/latest/。 ------------------------------------------- 2. ISSM runme.m文件包含了运行所有ISSM试验的全部指令。 output/文件夹存放所有试验的mat格式文件，数组说明同上。 ISSM的源代码与文档可通过https://issm.jpl.nasa.gov/获取。