_RL diagCG_resTarget
& diagCG_resTarget, diagCG_pcOffDFac
diagNormaliseRHS = diagCG_resTarget.GT.0.
residCriter = diagCG_resTarget
residCriter = a2dNorm * ABS(diagCG_resTarget)
diagCG_resTarget = cg2dTargetResidual
CALL WRITE_0D_RL( diagCG_resTarget, INDEX_NONE,
& diagCG_maxIters, diagCG_resTarget,