Added Field3D::applyBoundary(region, condition)

e.g. applyBoundary("target", "neumann")
Fixes for DND equilibrium communications

Problems duplicating MPI_Comm handles. Work-around to
copy rather than duplicate.