conservation laws can now take a function rho for a factor in the mass matrix term M = int_ rho Ni Nj simply add a function * massFactor in your conservation law to use it :-)