density
inverse of density regardless compressibility
function pointer to get 1/rho as a mask regardless incompressibility