fix a bug in the assembly so that now the bulk density is used on the rhs of the first equation
change pressure element to discontinous