a2dNorm = 1. _d 0
IF ( a2dMax .GT. 0. _d 0 ) a2dNorm = 1. _d 0/a2dMax
b2dNorm = a2dNorm
x2dNorm = a2dNorm*rhsMax
residCriter = a2dNorm * ABS(diagCG_resTarget)
aW2d(i,j,bi,bj) = aW2d(i,j,bi,bj)*a2dNorm
aS2d(i,j,bi,bj) = aS2d(i,j,bi,bj)*a2dNorm
& ' diag_cg2d (it=', myIter,') a2dNorm,x2dNorm=', a2dNorm,
_RL a2dMax, a2dNorm