Boundary Conditions

How do I create a Dirichlet boundary condition?

A Dirichlet boundary condition sets a constant value for the conserved quantities at the boundaries. A list containing the constant values for stage, xmomentum and ymomentum is constructed and used in the function call, e.g. Dirichlet_boundary([0.2,0.,0.]) .