# 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