Add forcing options for power of mutost.

Otherwise the function return an error for rejected particles with
rho.mc > 2*tau.u