Edit on GitHub

FastIsostasy.jl

❄ Fast and friendly glacial isostatic adjustment on CPU and GPU.

FastIsostasy is a friendly and flexible software package for simulations of regional glacial isostatic adjustment (GIA) with laterally-variable mantle viscosity and lithospheric thickness. It is mainly adressed to ice-sheet modellers who seek for (1) a good representation of solid-Earth mechanics at virtually zero computational cost, (2) an approximation of the sea-level equation and (3) ready-to-use inversion tools to calibrate the model parameters to data. The simple interface of FastIsostasy allows to flexibly solve GIA problems within few lines of code. It is fully open-source under MIT license and was succesfully benchmarked against analytical, 1D GIA and 3D GIA model solutions.

FastIsostasy relies on a hybrid Fourier/finite-difference collocation of the problem introduced in [Cathles1975] and solved in [Lingle1985], [Bueler2007]. Thanks to a simplification of the full problem from 3D to 2D space and the use of optimized software packages, running kiloyears of regional GIA with $\Delta x = \Delta y \simeq 45 \, \mathrm{km}$ is a matter of seconds on a single CPU. For high resolution runs, the user can switch to GPU usage with minimal syntax change and enjoy the advantage of parallelization without requiring an HPC cluster. The central place of Fast-Fourier transforms in FastIsostasy's solving scheme inspired its name, along with a GitHub repository that eased the first steps of this package. For GIA "purists", this package is likely to miss interesting processes but we belive that its ridiculous run-time can help to fast-prototype a problem before transfering it to a more comprehensive model.

Star us on GitHub!

If you have found this library useful, please consider starring it on GitHub. This gives us a lower bound of the satisfied user count.

Getting started

FastIsostasy.jl is work under devlopment and is not a registered julia package yet. To install it, please run:

using Pkg
Pkg.add("https://github.com/JanJereczek/FastIsostasy.jl")

How to read the docs?

If you want a quick introduction to GIA, please go to A quick introduction to GIA. If you want to have a thorough but still accessbile introduction to GIA, we highly recommend reading [Whitehouse2019]. If you want to get started right away, feel free to directly go to the Examples. If you face any problem using the code or want to know more about the functionalities of the package, visit the API reference. If you face a problem you cannot solve or have feature requests, please open a GitHub issue with a minimal and reproduceable example.

Julia ecosystem

FastIsostasy.jl was written thanks to the shear amount of work that people invested in the vast and well-documented Julia ecosystem. Major help from packages deserves major appreciation, in particular for: