# Benchmark for Stokes solver in Bousinessq approximation
[CitcomS]
solver = full
[CitcomS.solver]
stokes_flow_only = on
rayleigh = 1
[CitcomS.solver.mesher]
levels = 5
[CitcomS.solver.vsolver]
Solver = multigrid
## 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