Examples

The repository contains ready-to-run examples in the top-level examples directory.

For dedicated walkthrough pages, see:

• 2D Example

• 3D Example

• MTK_GMG Style

• MTK_GMG Example 1

• MTK_GMG Example 2 (Unzen)

• MTK_GMG 3D Examples (Unzen and Lanin)

Main Scripts

• Example2D.jl: 2D thermal evolution with repeated intrusions.

• Example3D.jl: 3D run with VTK output for ParaView.

• Example2D_ZASSy.jl: benchmark-style setup used in published comparisons.

• MTK_GMG_2D_example1.jl and MTK_GMG_2D_example2.jl: workflows integrating GeophysicalModelGenerator input.

Visual Outputs

2D example:

3D example:

MTK_GMG 3D examples:

Running an Example

From the repository root:

include("examples/Example2D.jl")

or

include("examples/Example3D.jl")

Code Snippets

The full README-equivalent code snippets are available in dedicated tabs:

• 2D Example

• 3D Example

Notes

• GPU execution typically requires loading CUDA before backend initialization.

• environment!(...) initializes package internals; if an example script uses ParallelStencil macros directly (@zeros, @parallel, etc.), it should also call @init_parallel_stencil(...) in script scope.

• 3D output is commonly explored with ParaView using VTK and PVD files.