ThicknessSmoothing
When this parameter is defined to a
non-zero value E, the
smoothing length for the evaluation over the mesh of the potential
arising from each planet is E*H, where H is the (linearly interpolated)
value of the disk thickness at the location of the planet. This
same disk thickness fraction applies to all planets simultaneously. The
use of this parameter is mutually exclusive with the use of RocheSmoothing