# Benchmark for Stokes solver in Truncated Anelastic Liquid approximation

[CitcomS]
solver = full


[CitcomS.solver]
stokes_flow_only = on
rayleigh = 1
dissipation_number = 0.5
gruneisen = 1


[CitcomS.solver.mesher]
levels = 5


[CitcomS.solver.vsolver]
Solver = multigrid
uzawa = bicg
#uzawa = cg


## This combination of ic and bc makes T=0 everywhere
## except one spherical harmonic load.
[CitcomS.solver.ic]
tic_method = 90
perturbl = 3
perturbm = 2


[CitcomS.solver.bc]
bottbcval = 0


[CitcomS.solver.output]
output_optional = surf, botm, geoid
self_gravitation = on
use_cbf_topo = off