Parsed fem dirichlet

This seems to be working as intended now.  In particular we can locally switch between Dirichlet and Neumann temperature boundaries based on whether or not the velocity indicates inflow or outflow.  In the long run that might be a reasonable condition to hard-code, but in the short run I suspect the CPU cost of recomputing constraints, etc. will outweigh the cost of calling fparser.