# Test the viscosity refinement criterion
set Adiabatic surface temperature          = 1623               # default: 0
set Dimension                              = 2

set End time                               = 0

subsection Boundary temperature model
  set List of model names = initial temperature

  subsection Initial temperature
    set Maximal temperature = 3773
    set Minimal temperature = 273
  end

end


subsection Geometry model
  set Model name = box # default: 

  subsection Box
    set X extent  = 500000
    set Y extent  = 500000
  end

end


subsection Gravity model
  set Model name = vertical # default: 

  subsection Vertical
    set Magnitude = 10.0
  end

end


subsection Initial temperature model
  set Model name = function # default: 
  subsection Function
    set Function constants  = l=250000
    set Function expression = if(x < l, 1873, 1623)
    set Variable names      = x,y
  end
end


subsection Material model

  set Model name = simple # default: 
  subsection Simple model
    set Thermal viscosity exponent = 5.0
    set Reference temperature = 1750
  end
end


subsection Mesh refinement
  set Coarsening fraction                      = 0.05
  set Refinement fraction                      = 0.3

  set Initial adaptive refinement              = 3                    # default: 2
  set Initial global refinement                = 2                    # default: 2

  set Strategy                                 = viscosity
end


# The parameters below this comment were created by the update script
# as replacement for the old 'Model settings' subsection. They can be
# safely merged with any existing subsections with the same name.

subsection Boundary temperature model
  set Fixed temperature boundary indicators   = 0,3        # default: 
end

subsection Boundary velocity model
  set Tangential velocity boundary indicators = 0,1,2,3
end

subsection Boundary velocity model
  set Zero velocity boundary indicators       =           # default: 
end


subsection Postprocess
  set List of postprocessors = 
end