point_value_02.prm

# A variation of the cookbook/sinker-with-averaging/sinker-with-averaging.prm
# file for evaluating the solution at chosen points
#
# This variation tests that we correctly convert output to years as well.

set Dimension                              = 2
set Start time                             = 0
set End time                               = 0.1
set Use years in output instead of seconds = true

set Linear solver tolerance                = 1e-7
set Pressure normalization                 = volume
set Use conduction timestep		   = false

subsection Geometry model
  set Model name = box
  subsection Box
    set X extent  = 1.0000
    set Y extent  = 1.0000
  end
end

subsection Model settings
  set Tangential velocity boundary indicators = 
  set Zero velocity boundary indicators       = left, right, bottom, top
end

subsection Material model
  set Model name = simple

  subsection Simple model
    set Reference density             = 1
    set Viscosity                     = 1
    set Thermal expansion coefficient = 0.0
    set Composition viscosity prefactor = 1e6
    set Density differential for compositional field 1 = 10
  end

  set Material averaging = harmonic average
end

subsection Gravity model
  set Model name = vertical
  subsection Vertical
    set Magnitude = 1
  end
end


############### Parameters describing the temperature field
# Note: The temperature plays no role in this model

subsection Boundary temperature model
  set List of model names = box
end

subsection Initial temperature model
  set Model name = function
  subsection Function
    set Function expression = 0
  end
end


############### Parameters describing the compositional field
# Note: The compositional field is what drives the flow
# in this example

subsection Compositional fields
  set Number of fields = 1
end

subsection Initial composition model
  set Model name = function
  subsection Function
    set Variable names      = x,y
    set Function expression = if( (sqrt((x-0.5)^2+(y-0.5)^2)>0.22) , 0 , 1 )
  end
end


############### Parameters describing the discretization

subsection Mesh refinement
  set Initial global refinement          = 4
  set Initial adaptive refinement        = 0
end



############### Parameters describing what to do with the solution

subsection Postprocess
  set List of postprocessors = visualization, point values

  subsection Point values
    set Evaluation points = 0.25, 0.5  ; \
                            0.5 , 0.25 ; \
                            0.75, 0.75
  end

  subsection Visualization
    set Output format                 = vtu
    set Time between graphical output = 0
    set List of output variables = density, viscosity
  end
end