picardResidual = 0. _d 0
& picardResidual,myThid)
picardResidual = SQRT(picardResidual)
picardResidualKm1 = picardResidual
picardConverged = picardResidual.LT.picardTol
& .OR.picardResidual.EQ.0. _d 0
FGMRESeps = krylovLinTol * picardResidual
& krylovLinTol,picardResidual
picardTol=SEAICEnonLinTol*picardResidual
& JFNKres_t = picardResidual * JFNKres_tFac
_RL picardResidual