PLEASE READ FOR INFORMATION ON DATA All data are in *.mat format, which contain Ice-Sheet and Sea-Level System Model 500-yr transient output for Thwaites and Pine Island glaciers as presented in the associated paper. Please refer to Schwans et al., 2023 for the particulars of each setting described below and finer details regarding workflow and output. Full list of models can be found in the publication's appendix. Naming convention for output: c/f = mesh resolution; either coarse (c) or fine (f) as described in section 2.2 m1/m8 = choice of basal rheology via power-law sliding exponent, m; either 1 (linear-viscous) or 8 (plastic) (see section 2.2.2) h/l/p = sub-shelf melt scheme used; either static ocean forcing with high (h) or low (l) inland-extrapolated values or dynamic forcing using PICOP module (p) from Pelle et al., 2019. (see section 2.3.1) subelon/subeloff = indicates whether subelement friction was enabled in model run if = ice-front prescription; either static (ifstat), dynamic (ifdyn), or dynamic with the basal friction coefficient, alpha, extrapolated seaward (ifdynxalpha) (see section 2.2.1/2.2.2) smf/smb = shear-margin perturbation, if applied; either reduction in basal friction (smf) or softening of the shear-margin ice fabric (smb) (see section 2.3.2) EX: f-m8-p-subelon-ifdyn.mat is a model that uses a fine mesh and simulates ice-flow over a plastic bed with the sliding-law coefficient set to 8, where ocean forcing is calculated by the PICOP module. This model resolves basal friction on a subelement basis, and uses a dynamic ice-front prescription, but without friction coefficient extrapolation. No shear-margin perturbation is applied.