This repository contains research code for "Simulation-Based Inference of Surface Accumulation and Basal Melt Rates of an Antarctic Ice shelf from Isochronal Layers". The repository used here is used to preprocess ice shelf thickness and velocity data from the real world, and to create synthetic ice shelves, in order to extract realistic transects. These transects are used in the main repository for the paper. The picked IRH transect data required to run calculate_flowtube.mis found here.

This repository requires an installation of the finite element solver firedrake, and the ice flow solver icepack to be installed in this environment. We direct users to follow the installation instructions of firedrake first.

Once a working firedrake environment is available on your machine, activate the virtual and environment and install this repository with:

  git clone https://github.com/mackelab/preprocessing-ice-data.git
  cd preprocessing_ice_data
  pip install -e setup.py

Create and spin up synthetic ice shelves using python preproc/synthetic_flowlines.py. Configurations for the synthetic ice shelf can be controlled with configs/.

Loading and processing of data from Ekström Ice Shelf can be done with python preproc/loading_flowlines.py. Note that you need access to a version of a Bedmachine dataset and an its_live velocity map.

Creating the flow tube coordinates for Ekström ice shelves is done with preproc/calculate_flowtube.m.