# Listing of Parameters
# ---------------------

# At the top, state that the program should be run for a two-dimensional 
# simulation
set Dimension = 2 

# CFL number (related to time step)
set CFL number                             = 1.0 

# The end time of the simulation. Units: years if the 'Use years in output
# instead of seconds' parameter is set; seconds otherwise.
set End time                               = 1e10 

# Output directory (may be an absolute or a relative path).
set Output directory                       = output_b3_2ds_phi360

# Resume computation parameter (true or false)
set Resume computation                     = false

# The start time of the simulation. Units: years if the 'Use years in output
# instead of seconds' parameter is set; seconds otherwise.
set Start time                             = 0 

# Surface adiabatic parameters (note that adiabatic surface temperature does
# not have to match surface temperature boundary condition)
set Adiabatic surface temperature          = 0 
set Surface pressure                       = 0 

# Time units in years instead of seconds (true or false)
set Use years in output instead of seconds = true  # default: true

# Select one of the following solver schemes:
set Nonlinear solver scheme                = IMPES

# Temperature boundary model
subsection Boundary temperature model
  set Model name = spherical constant
  subsection Spherical constant
    # Temperature at the inner boundary (core mantle boundary). Units: K.
    set Inner temperature = 2850
    # Temperature at the outer boundary (lithosphere water/air). Units: K.
    set Outer temperature = 300
  end
end


subsection Discretization
  # The polynomial degree to use for the velocity variables in the Stokes system.
  set Stokes velocity polynomial degree       = 2
  # The polynomial degree to use for the temperature variable
  set Temperature polynomial degree           = 2
  set Use locally conservative discretization = false

  subsection Stabilization parameters
    # The exponent $\alpha$ in the entropy viscosity stabilization.
    set alpha = 2
    # The $\beta$ factor in the artificial viscosity stabilization. An
    # appropriate value for 2d is 0.052 and 0.078 for 3d. Units: None.
    set beta  = 0.078
    # The $c_R$ factor in the entropy viscosity stabilization. Units: None.
    set cR    = 0.11
  end

end

# Geometry section
subsection Geometry model
  set Model name = spherical shell
  subsection Spherical shell
    # Inner radius of the spherical shell. Units: m.
    set Inner radius  = 3480000
    # Opening angle in degrees of the section of the shell that we want to
    # build. Units: degrees.
    set Opening angle = 360
    # Outer radius of the spherical shell. Units: m.
    set Outer radius  = 6370000
  end

end

# Gravity section
subsection Gravity model
  set Model name = radial earth-like
  subsection Radial constant
    set Magnitude = 10.0
  end

end

# Initial conditions
subsection Initial conditions
  set Model name = spherical hexagonal perturbation
  subsection Spherical gaussian perturbation
    # The amplitude of the perturbation.
    set Amplitude             = 0.01
    # The angle where the center of the perturbation is placed.
    set Angle                 = 0e0
    # The non-dimensional radial distance where the center of the perturbation
    # is placed.
    set Non-dimensional depth = 0.7
    # The standard deviation of the Gaussian perturbation.
    set Sigma                 = 0.2
    # The sign of the perturbation.
    set Sign                  = 1
  end
end


subsection Material model
  set Model name = simple
  subsection Simple model
    # Reference density $\rho_0$. Units: $kg/m^3$.
    set Reference density             = 4500
    # The value of the specific heat $cp$. Units: $JG/kgK$.
    set Reference specific heat       = 1000    
    # The reference temperature $T_0$. Units: $K$.
    set Reference temperature         = 300
    # The value of the thermal conductivity $k$. Units: $W/m/K$.
    set Thermal conductivity          = 4.0
    # The value of the thermal expansion coefficient $\beta$. Units: $1/K$.
    set Thermal expansion coefficient = 2.5e-5
    # The value of the constant viscosity. Units: $kg/m/s$.
    set Viscosity                     = 1.00e23
  end
end


subsection Mesh refinement
  set Additional refinement times        = 
  # The number of adaptive refinement steps performed after initial global
  # refinement but while still within the first time step.
  set Initial adaptive refinement        = 3

  # The number of global refinement steps performed on the initial coarse
  # mesh, before the problem is first solved there.
  set Initial global refinement          = 4

  # The fraction of cells with the largest error that should be flagged for
  # refinement.
  set Refinement fraction                = 0.3

  # The fraction of cells with the smallest error that should be flagged for
  # coarsening.
  set Coarsening fraction                = 0.05

  # The method used to determine which cells to refine and which to coarsen.
  set Strategy                           = Temperature

  # The number of time steps after which the mesh is to be adapted again based
  # on computed error indicators.
  set Time steps between mesh refinement = 5
end


subsection Model settings
  set Fixed temperature boundary indicators   = 0,1
  set Include shear heating                   = false
  set Prescribed velocity boundary indicators = 
  set Radiogenic heating rate                 = 0e0
  set Tangential velocity boundary indicators = 0, 1
  set Zero velocity boundary indicators       = 
end


subsection Postprocess
  set List of postprocessors = visualization,velocity statistics,temperature statistics,heat flux statistics, depth average
  subsection Visualization
    set Number of grouped files       = 0
    # The file format to be used for graphical output.
    set Output format                 = vtu
    # step. Units: years if the 'Use years in output instead of seconds'
    # parameter is set; seconds otherwise.
    set Time between graphical output = 1e8
  end
  subsection Depth average
    set Time between graphical output = 1e8
  end
end

subsection Checkpointing
  set Time between checkpoint = 0
  set Steps between checkpoint = 50
end