| 80 | 
 C     u/vIceRes :: residual of sea-ice momentum equations | 
 C     u/vIceRes :: residual of sea-ice momentum equations | 
| 81 | 
       _RL uIceRes(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
       _RL uIceRes(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
| 82 | 
       _RL vIceRes(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
       _RL vIceRes(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
| 83 | 
  | 
 C     vector version of the residuals | 
| 84 | 
  | 
       _RL resTmp (nVec,1,nSx,nSy) | 
| 85 | 
 C     du/vIce   :: ice velocity increment to be added to u/vIce | 
 C     du/vIce   :: ice velocity increment to be added to u/vIce | 
| 86 | 
       _RL duIce  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
       _RL duIce  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
| 87 | 
       _RL dvIce  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
       _RL dvIce  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
| 92 | 
       _RL etaZPre (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
       _RL etaZPre (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
| 93 | 
       _RL dwatPre (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
       _RL dwatPre (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | 
| 94 | 
 CEOP | 
 CEOP | 
 | 
       INTEGER n | 
  | 
 | 
       _RL resTmp (nVec,1,nSx,nSy) | 
  | 
| 95 | 
  | 
  | 
| 96 | 
 C     Initialise | 
 C     Initialise | 
| 97 | 
       newtonIter          = 0 | 
       newtonIter          = 0 |