[instat. advection] use numpymatrix with sparse identity for mass