Amundsen Sea Embayment ice sheet model simulations 2021-2300 forced with Paris2C and RCP8.5 climate scenarios

This repository contains model simulation results for the Amundsen Sea Embayment region from 2021 until the year 2300 obtained using the Ua ice-flow model. Each file represents a single simulation and is named according to the climate forcing (either RCP8.5 or Paris 2C) along with a unique 4 digit number identifying the simulation number. Note: Since the finite element mesh is constantly adapting, every timestep will have a different number of nodes, different nodal coordinates and a different finite element connectivity. All fields in the netcdf files are padded with NaNs in each timestep to account for this changing field size.

The contents of each file are as follows:

time: Simulation time in years (outputs saved every 2 model years)

/modelFields: Contains all main evolving model fields of interest

/modelFields/s: Ice sheet surface elevation (m)

/modelFields/h: Ice thickness (m)

/modelFields/b: Ice sheet base (m)

/modelFields/ub and vb: Ice sheet velocity in the u and v directions (m/yr)

/modelFields/as and ab: Mass balance at the ice sheet surface and base, respectively (m/yr)

/nodalCoordinates/x and y: The x and y coordinates of all nodes in the finite element mesh

/meshConnectivity: Contains three connectivty matrices to form the finite element triangulation. To assemble the triangulation for one point in time run the following in matlab: conn = cat(1,connectivity1(t,:),connectivity2(t,:),connectivity3(t,:))';

/parameters: Contains scalar values for all values of model parameters for this particular experiment