Copyright 2010--2024 Ed Bueler

This repository contains slides, notes, and computer programs about numerical glacier and ice sheet modeling. They have been used for the International Summer School in Glaciology, McCarthy, AK in years 2010, 2012, 2014, 2016, 2018, and 2022. The PDF notes (.pdf in notes/) plus their exercises plus the Matlab/Octave programs form a self-contained course. The PDF slides (.pdf in slides/) cover the same material in a more informal style.

This content will be completely revised for the June 2024 school. The current state is, at best, of historical interest.

The codes in subdirectory mfiles/ solve SIA and SSA problems. They should work in either Matlab or Octave; if not please report a bug, either by email or by using the issues for this repository. Download them either by cloning this repository or by getting a .zip or .tar.gz archive at the releases page. The PDFs mentioned above document these programs, but the programs also have help files (i.e. leading comments). You are encouraged to actually run and modify them!

The Python tools in stokes/ solve a free-surface 2D Glen-Stokes flow over a bedrock step. This is used to generate some images in the slides and for student projects. Note that the numerical technology used here is more advanced, and more experimental, than the mfiles/ content. The workflow uses the following tools: Firedrake finite element library, Gmsh mesh generator, PETSc solver library, and Paraview visualization. See stokes/README.md for more information.

Older versions (2009, 2010, 2012, 2014) lived in the repo https://github.com/bueler/karthaus