stress_build-up_maxwell_2d.prm

#  Global parameters
set Dimension                              = 2
set Start time                             = 0
set End time                               = 1e6
set Use years in output instead of seconds = true
set Use direct solver for Stokes system    = false
set Linear solver tolerance                = 1e-7
set Nonlinear solver scheme                = IMPES
set Nonlinear solver tolerance             = 1e-4
set Max nonlinear iterations               = 10
set Number of cheap Stokes solver steps    = 500
set CFL number                             = 0.5
set Maximum time step                      = 1000
set Output directory                       = output
set Timing output frequency                = 1
set Pressure normalization                 = surface
set Surface pressure                       = 0.

# Model geometry (100x100 km, 1 km spacing)
subsection Geometry model
  set Model name = box

  subsection Box
    set X repetitions = 100
    set Y repetitions = 100
    set X extent      = 100e3
    set Y extent      = 100e3
  end
end

# Mesh refinement specifications 
subsection Mesh refinement
  set Initial adaptive refinement        = 0
  set Initial global refinement          = 0
  set Time steps between mesh refinement = 0
end

# Boundary classifications
subsection Model settings
  set Fixed composition boundary indicators   = bottom, top, left, right
  set Fixed temperature boundary indicators   = bottom, top, left, right
  set Prescribed velocity boundary indicators = bottom y:function, right x:function 
  set Free surface boundary indicators        = 
  set Tangential velocity boundary indicators = top, left 
  set Zero velocity boundary indicators       = 
  set Enable elasticity                       = true
end

subsection Boundary velocity model
  subsection Function
    set Variable names      = x,y
    set Function constants  = cm=0.01, year=1
    set Function expression = if (x<50e3 , -3.154*cm/year, 0.); if (y<50e3 , -3.154*cm/year, 0.); 
  end
end

# Number and name of compositional fields
subsection Compositional fields
  set Number of fields = 3
  set Names of fields  = s11, s22, s12
end

# Spatial domain of different compositional fields
subsection Initial composition model
  set Model name = function
  subsection Function
    set Variable names      = x,y
    set Function constants  = 
    set Function expression = 0; 0; 0;
  end
end

# Composition boundary conditions
# (Note: should change boundary values to be defined using analytical solution)
subsection Boundary composition model
  set Model name = initial composition
end

# Temperature boundary and initial conditions
subsection Boundary temperature model
  set Model name = box
  subsection Box
    set Bottom temperature = 273
    set Left temperature   = 273
    set Right temperature  = 273
    set Top temperature    = 273
  end
end
subsection Initial temperature model
  set Model name = function
  subsection Function
    set Function expression = 273
  end
end

# Material model
subsection Material model

  set Model name = viscoelastic

  subsection Viscoelastic
    set Densities            =  2800
    set Viscosities          = 1.e22
    set Elastic shear moduli = 1.e10
    set Elastic time step = 1e3
    set Viscosity averaging scheme = harmonic
  end

end

# Gravity model
subsection Gravity model
  set Model name = vertical
  subsection Vertical
    set Magnitude = 0.
  end
end

# Element types
subsection Discretization
  set Composition polynomial degree     = 2
  set Stokes velocity polynomial degree = 2
  set Temperature polynomial degree     = 1
end

# Post processing
subsection Postprocess
  set List of postprocessors = basic statistics, composition statistics, temperature statistics, velocity statistics, visualization
  subsection Visualization
    set List of output variables = material properties, strain rate 

    subsection Material properties
      set List of material properties = density, viscosity
    end

    set Time between graphical output = 1e3
    set Interpolate output = true
  end

end

# Termination criteria
subsection Termination criteria
 set End step = 1000
 set Termination criteria = end time
end


subsection Checkpointing
  set Time between checkpoint = 100
end