* **Source Fortran File:** `MshGlacierSynthetic.f90`

This tool deforms vertically a 1m thick 2D or 3D mesh to conform with the bedrock and surface elevations given as two analytical expression in `fbed.f90`

and `fsurf.f90`

, respectively. These fortran functions are of the form:

FUNCTION fbed(x,y) USE types IMPLICIT NONE REAL(KIND=dp) :: x, y, fbed fbed = ... END FUNCTION fbed

For a 2D geometry, `fbed`

and `fsurf`

are only function of `x`

.

Informations (name and number of partitions) of mesh are given is the file `mesh_info.in`

. This tool works for already partitioned mesh.

Compilation is done as follow:

elmerf90-nosh MshGlacierSynthetic.f90 fbed.f90 fsurf.f90 -o MshGlacierSynthetic

An example using `MshGlacierSynthetic`

to construct ISMIP-HOM experiment A mesh can be found in `[ELMER_TRUNK]/elmerice/examples/Test_MshGlacierSynthetic`

.