Fix J Integral adjacent to free surfaces ref #3898

If a crack intersects a free surface, we need to zero out the q function
on that surface for all crack front nodes except the one on that surface.
Add an option to list a set of free surfaces intersected by the crack
to do this.