Operator methods handling boundary grids

[Ivar Stefansson]

All methods returning an operator which might correspond to a boundary value need to handle boundary grids and return a

• TimeDependentDenseArray tied to the prescribed value if domain is boundary
• MixedDimVariable/other Operator if domain is internal