Inverse methods test cases

• Modification date : 04/07/2024

This test cases for the adjoint inverse methods have been updated in April 2020 from the material presented for the Elmer/Ice course held in Oslo in 2016.

The description of the test cases can be found in the associated presentation available here. Note that the solvers have been updated so that the implementation differs and use the new solvers for the inverse methods (see documentation here).

A new Regularisation solver that can be used to prescribed prior error statistics under the form of a covariance matrix has been introduced in 07/2024. The documentation is available here. A test case for the Mass conservation Method can be found here.

• Content of this directory:

  • DATA: data sets and processing tools required to run the experiments.

  • src: user functions needed by the test cases.

  • SCRIPTS: python scripts used for post-processing.

TEST CASES USING THE SSA SOLVER:

• MacAyeal_SSA: Optimisation of the basal friction coefficient based on a synthetic test case.

• RonneFilchner_SSA: Optimisation of the viscosity on the Ronne-Filchner ice-shelf. You will need data-sets to run this experiment, see in the DATA directory.

• RonneFilchner2_SSA: Optimisation of both the viscosity and fiction on the Ronne-Filchner ice-shelf (including ice rises). External data sets are required to run this experiment, see in the DATA directory.

TEST CASES USING THE STOKES SOLVER:

• MacAyeal_Stokes: Optimisation of the basal friction coefficient based on a synthetic test case.

• StokesWeertman: non-linear weertman friction law with the vectorised incompressible Stokes solver

TEST CASES FOR THE THICKNESS SOLVER:

• MassConservation: Validation of the adjoint of the steady state thickness solver; i.e. the mass conservation method.