'yc = 0.0

xc = 5.0

nsel, s, loc, x, 0 , xc

esln,s,

esel, r, cent, y, -1e-3, 1e-3

cm, cenelem, elem

nelem = 1000

iel = 0

Phi = 0.0

Psi = 0.0

*do, i, 1, nelem, 1

iel = elnext(iel)

*if, iel, ne, 0, then

*do, j, 1, 4, 1

nd = nelem(iel,j)

Phi = ny(nd) - yc

Psi = nx(nd) - xc

xfdata, ENRICH1, LSM, iel, nd, Phi, Psi

*enddo

*endif

*enddo'

This code is extracted from APDL, XEFM based stationary crack analysis.

I didn't understand this looping section of the code and would mean a lot if someone from the related field could clear my confusion, as I am doing my final year project related to crack growth in pelton bucket.