U ndiagmx,
mdiag(mId,listId) = ndiagmx + 1
ndiagmx = ndiagmx + kdiag(mate)*averageCycle(listId)
IF ( ndiagmx.GT.numDiags ) THEN
INTEGER ndiagmx
idiag(mId,listId) = ndiagmx + 1
ndiagmx = ndiagmx + kdiag(ndId)*averageCycle(listId)
IF (ndiagmx+kdiag(mate).GT.diagSt_size) THEN
mSdiag(mId,listId) = ndiagmx + 1
ndiagmx = ndiagmx + kdiag(mate)
IF (ndiagmx+kdiag(ndId).GT.diagSt_size) THEN
iSdiag(mId,listId) = ndiagmx + 1
ndiagmx = ndiagmx + kdiag(ndId)