Boundary function specification

[bendudson]

would like to be able to specify more complex boundary conditions by specifying an algebraic expression. This is currently done for the initial conditions using the FieldFactory class, and would be good to extend this to boundary regions only.

Even better would be applying time-dependent boundary conditions this way, so in addition to 'x', 'y' and 'z' add 't' for time to the parser.