#  Global parameters
set Dimension                              = 2
set Start time                             = 0
set End time                               = 100
set Use years in output instead of seconds = true
set Nonlinear solver scheme                = single Advection, iterated Stokes
set CFL number                             = 0.5
set Maximum time step                      = 2.5
set Output directory                       = /cluster/home/maaike/from_box_to_shell/case_shell_nr_2/case_shell_nr_2
set Timing output frequency                = 1
set Pressure normalization                 = no

# Model geometry
subsection Geometry model
  set Model name = spherical shell
  subsection Spherical shell
    set Inner radius = 3480.e3
    set Outer radius = 6371.e3
	set Opening angle = 360.
  end
end

## Mesh refinement specifications 
#subsection Mesh refinement
#  set Initial adaptive refinement        = 4
#  set Initial global refinement          = 1
#  set Time steps between mesh refinement = 1
#  set Strategy                           = minimum refinement function
#  subsection Minimum refinement function
#    set Coordinate system   = spherical
#	set Variable names      = r,phi
#	set Function expression = (phi<=0.0157 && phi>=-0.0157) ? 4 : 1
#  end
#end

# Mesh refinement specifications 
subsection Mesh refinement
  set Initial adaptive refinement        = 1
  set Initial global refinement          = 1
  set Time steps between mesh refinement = 1
  set Strategy                           = strain rate
end


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

# Formulation classification
subsection Formulation
  set Enable elasticity = true
end

subsection Mesh deformation
  set Mesh deformation boundary indicators = outer: free surface
  subsection Free surface
    set Surface velocity projection = normal
  end
end

# Velocity boundary conditions
subsection Boundary velocity model
  set Tangential velocity boundary indicators = inner
end

subsection Nullspace removal
  set Remove nullspace = angular momentum
end

# Prescribe a fixed vertical traction on the top boundary
subsection Boundary traction model
  set Prescribed traction boundary indicators = outer: function
  subsection Function
    set Coordinate system = cartesian
    set Variable names = x,y,t
    set Function constants  = alpha=0.0157, H0=1.e3, rhoi=931., g=9.8, t1=100., Re=6371.e3 # alpha is halfamplitude load of 100 km here
    set Function expression = x>=Re*cos(alpha) && y<=Re*sin(alpha) && y>=Re*sin(-alpha) ? -g*H0*rhoi*t/t1*cos((atan(y/x)*180/pi)) : 0; x>=Re*cos(alpha) && y<=Re*sin(alpha) && y>=Re*sin(-alpha) ? y>=0 ? -g*H0*rhoi*t/t1*cos((atan(y/x)*180/pi)) : g*H0*rhoi*t/t1*cos((atan(y/x)*180/pi)) : 0;
  end
end

# Number and names of compositional fields
subsection Compositional fields
  set Number of fields = 3
  set Names of fields = ve_stress_xx, ve_stress_yy, ve_stress_xy
end

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
subsection Boundary composition model
  set Allow fixed composition on outflow boundaries = true
  set List of model names = initial composition
end

# Temperature boundary conditions
subsection Boundary temperature model
  set Fixed temperature boundary indicators = inner, outer
  set List of model names = initial temperature
end

# Temperature initial conditions
subsection Initial temperature model
  set Model name = function
  subsection Function
    set Function expression = 293
  end
end

# Material model
subsection Material model
  set Model name = viscoelastic
  subsection Viscoelastic
    set Densities                   = 4450
    set Viscosities                 = 5.0e20
    set Elastic shear moduli        = 175.e9
    set Use fixed elastic time step = false
    set Fixed elastic time step     = 2.5
    set Use stress averaging        = false
    set Viscosity averaging scheme  = harmonic
  end

end

# Gravity model
subsection Gravity model
  set Model name = radial constant
  subsection Radial constant
    set Magnitude = 9.8
  end
end

# Post processing
subsection Postprocess
  set List of postprocessors = basic statistics, composition statistics, temperature statistics, topography, velocity statistics, visualization
  
  subsection Visualization
    set List of output variables = material properties, nonadiabatic pressure, strain rate, spherical velocity components
    subsection Material properties
      set List of material properties = density, viscosity
    end
    set Time between graphical output = 0
    set Interpolate output = false
    set Number of grouped files = 1
  end

  subsection Topography
     set Output to file = true
     set Time between text output = 0
  end
end

# Termination criteria
subsection Termination criteria
 set Termination criteria = end time
end