Moving reprojections to PROBLEM_SOLVER_RESIDUAL_EVALUATE

Tracker item: 3675
Output field and projection results at each iteration for single region contact problems

Tracker item: 3675
Merge remote-tracking branch 'PrasadBabarendaGamage/single_region_contact' into linearisation

Conflicts:
	src/problem_routines.f90

Resolved the conflict temp, the next commit will add functions necessary to complete the iteration-dependent geometric term addition

Tracker item: 3675
Add iteration number when evaluating the contact metrics. For supporting only adding geometric term to contact stiffness matrix after some Newton iterations

Tracker item: 3675