continental_extension.prm

####
# This test is based on the Continental Extension Cookbook
# There, the initial temperature is computed through a
# somewhat cumbersome function specified in the input file.
# With the continental geotherm plugin, one only has to
# specify the boundary temperature constraints and layer
# thicknesses of the upper crust, lower crust and mantle
# lithosphere. Other parameters such as layer densities
# are read from other sections in the input file, ensuring
# consistency between the parameters set for the initial
# temperature and those in the material model and compositional
# heating plugins. Significantly, the continental geotherm
# initial temperature plugin enforces that specific names are
# used for crustal and mantle layers, which are specified
# below in the Compositional fields subsection.

include $ASPECT_SOURCE_DIR/cookbooks/continental_extension/continental_extension.prm

set Timing output frequency                = 100
set End time                               = 1e4
set Max nonlinear iterations               = 1

# Perform only one global refinement and no adaptive refinements
subsection Mesh refinement
  set Initial adaptive refinement        = 0
  set Initial global refinement          = 1
  set Time steps between mesh refinement = 0
end

# Update the names of the compositional fields associated with the upper crust, lower crust,
# and lithospheric mantle to be consistent with the naming scheme enforced by the continental
# geothermal initial temperature plugin.
subsection Compositional fields
  set Names of fields = noninitial_plastic_strain, plastic_strain, upper_crust, lower_crust, lithospheric_mantle
end

# Initial temperature field
subsection Initial temperature model
  set Model name = continental geotherm

  subsection Continental geotherm
    set Lithosphere-Asthenosphere boundary isotherm = 1613.
    set Surface temperature = 273.
    set Layer thicknesses = 20000., 20000., 60000.
  end
end

# Post processing
subsection Postprocess
  set List of postprocessors = temperature statistics
end