Add an optional argument to O.forces.f() and O.forces.t() to get correct forces on clumps.