Enabled integration point data I/O in THM
The added functionality in this MR allows ogs to use the stress results at the integration points of a THM simulation as the initial stresses of other THM simulation. Besides, this MR added a shape function order checking in the process input data.
Feature description was added to the changelog
Tests covering your feature were added? Yes. A new benchmark with IP stored initial stress and BGRa creep model, which is based on a test of the ARESH project, is added.
Any new feature or behavior change was documented?