Added function to compute flux in direction of normal
Store only diagonal mass matrix.
Make top wall to be slip wall