This is done through the `VariableBCDescriptor` type in C++ which reads
`variable_boundary_condition` in the scripting language.

Also add the following pgs-function:
  variable_bc: Vh * (boundary_condition) -> variable_boundary_condition