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