Crank-Nicolson is a popular method for solving parabolic equations because it is unconditionally stable and second-order accurate. One drawback of CN is that it responds to jump discontinuities in the initial conditions with oscillations which are weakly damped and therefore may persist for a long time. We compare a selection of methods to reduce the amplitude of these oscillations.