############### Global parameters
# This is a test for the radiogenic heating and the internal
# heating statistics postprocessor. Since density and radiogenic
# heating are constant there is an analytical solution, which
# is reproduced in this testcase. 

set Dimension = 2

set Start time                             = 0
set End time                               = 1000
set Use years in output instead of seconds = false

set Output directory = radiogenic_heating

subsection Heating model
  set List of model names = constant heating

  subsection Constant heating
    set Radiogenic heating rate = 1
  end
end

subsection Geometry model
  set Model name = box

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


subsection Gravity model
  set Model name = vertical
  subsection Vertical
    set Magnitude = 10.0
  end
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
  # We only fix the temperature at the upper boundary, the other boundaries
  # are isolating. To guarantuee a steady downward flow, we fix the velocity
  # at the top and bottom, and set it to free slip on the sides. 
  set Fixed temperature boundary indicators   = 3
end

subsection Boundary velocity model
  set Prescribed velocity boundary indicators = 2:function, 3:function
end

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


############### Boundary conditions
# We set the top temperature to T1=1000K. 
subsection Boundary temperature model
  set List of model names = box
  subsection Box
    set Top temperature = 1000
  end
end

# We prescribe a constant downward flow.
subsection Boundary velocity model
  subsection Function
    set Function expression = 0;-1
    set Variable names      = x,y
  end
end

subsection Initial temperature model
  set Model name = function
  subsection Function
    set Function expression = 1000.0
    set Variable names      = x,y
  end
end


subsection Material model
  set Model name = simple

  subsection Simple model

    set Reference density                              = 3400
    set Reference specific heat                        = 1000
    set Reference temperature                          = 1000
    set Thermal conductivity                           = 0.0

    set Thermal expansion coefficient                  = 0.0

    set Thermal viscosity exponent                     = 0.0
    set Viscosity                                      = 1e21
  end
end


subsection Mesh refinement
  set Initial adaptive refinement        = 0 
  set Initial global refinement          = 5
  set Time steps between mesh refinement = 0

end


subsection Postprocess

  set List of postprocessors = temperature statistics, heating statistics
  subsection Visualization
    set Interpolate output = false
    set List of output variables = heating
    set Time between graphical output = 0
  end

end