This type of results data is entered using routines DIEGH and DIEGR. For each element; routine DIEGH is called once only, and routine DIEGR is called once for each Gauss point at each surface. The number of Gauss points is the inner repeat, and the number of surfaces is the outer repeat when calling routine DIEGR. The co-ordinates for each Gauss point (in the global axes system) must be passed to routine DIEGR each time it is called.