| Summary: | Code accompanying manuscript submission to The Cryosphere journal. surfacespeed_ratios22.m computes surface speeds predicted by internal deformation and Weertman sliding mechanisms, comparing predictions from the two mechanisms to observed speeds. surfacespeed_ratios23.m repeats the computation fromsurfacespeed_ratios22.m, but with more smoothing on several variables input to the internal deformation and Weertman sliding models. This code requires that the user already has a workspace open with the latitude and longitude data for a flowline. surfacespeed_ratios24.m computesthe misfit between predictions of each mechanism and observed data. This code requires that the user load the results from surfacespeed_ratios22.m or surfacepseed_ratios23.m into the workspace, choosing a particular flowline and a region within it where the user has already concluded (qualitatively) that one mechanism provides closer agreement with observed data than the other. This code is to verify qualitative conclusions made on the figure comparisons produced by surfacespeed_ratios23.m and surfacespeed_ratios22.m. The ERA-Interim reanalysis data used to estimate the ice temperature at depth are publicly available, and were downloaded from https://www.ecmwf.int/en/forecasts/datasets/archive-datasets/reanalysis-datasets/era-interim . Ice thickness, ice elevation, bedrock topography data and relevant Matlab functions were sourced from the Bedmap2 data set, available from Matlab File Exchange: https://www.mathworks.com/matlabcentral/fileexchange/42353-bedmap2-toolbox-for-matlab . Ice velocity and grounding line data were sourced from MEaSUREs data set, available from Matlab File Exchange: https://www.mathworks.com/matlabcentral/fileexchange/47329-measures?s_tid=ta_fx_results .
|
|---|