Add mesh visualization functionality

[simulkade]

Not really important, but may be helpful for preparing presentations and documents.

[simulkade]

If boundary condition is passed to the mesh visualization, label each boundary with the boundary condition type, i.e., 'D' for Dirichlet, 'N' for Neumann, 'R' for Robin, and 'P' for periodic.