Implement output of fault tractions based on Lagrange multiplier faults

[baagaard-usgs]

This requires projection from the cohesive cells onto the fault submesh. This requires changes to OutputSubfield. We need to pass in the label and label value of the physics object.

TODO

• [ ] Output of change in fault tractions using Lagrange multiplier
• [ ] Output of normal and tangential directions