mccarthy

This repository contains slides, notes, and computer programs for numerical glacier and ice sheet modeling. These materials have been used for the International Summer School in Glaciology, McCarthy, AK in years 2010, 2012, 2014, 2016, 2018, 2022, and 2024.

Download all materials by one of these methods:

use the releases page for a .zip/.tar.gz archive
clone or shallow clone this repository:

git clone --depth=1 https://github.com/bueler/mccarthy.git

The optional depth=1 setting reduces the download size by not getting the history, which most users will not need.

slides and notes

The PDF slides slides/slides-2024.pdf are new this year. I believe they cover the essential material in a better way. They are supported by Python codes in py/.

The older material is in PDF notes notes/notes-2024.pdf, previous slides slides/slides-2022.pdf, and Matlab/Octave programs mfiles/. My plan is to update the notes soon!

Python programs

The codes in subdirectory py/ solve surface kinematic equation (SKE) and shallow ice approximation (SIA) problems. Python programs surface1d.py and shallowuw.py use only the numpy and matplotlib libraries. You are encouraged to actually run and modify them! Feel free to email me about how they work.

For now the slides (slides-2024.pdf) are the documentation, along with line comments in the .py source codes. Please report any bugs, either by email or by using the issues for this repository.

The surface2d.py program uses Firedrake. Additional Firedrake codes which solve the glaciological (Glen-law) Stokes model are in the stokes/ directory. See stokes/doc.pdf for some documentation.

Matlab/Octave programs (deprecated)

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. The current notes and the older slides document these programs, but the programs also have help files (leading comments).

Stokes solver

The Python tools in stokes/ are primarily for projects. By default they solve a 2D Glen-law Stokes flow over a bedrock step. The workflow uses the following tools: Firedrake (a finite element library), Gmsh (a mesh generator), PETSc (a solver library), and Paraview (for visualization). See stokes/doc.pdf for more information.

ancient versions

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

Mccarthy

Install / Use

README