assignment of boundary conditions in case of Euler flow

[highlando]

By now, the call of condense_sysmats_bybcs assumes that we gonna include the stiffness matrix.

For the numerical tests it didn't pose a problem, since v=0 at the boundary.