Not sure yet. I have made by own Q2 code and implemented your experiments in it as well as the Heaviside advection pb. I am trying to compare my results with those of ASPECT but I am using a Crank-Nicolson algo while ASPECT relies on the BDF2 with implicit advection. More soon.