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
Keeping up with Juliane's branch.