153 |
IF ( .NOT.SEAICEuseIMEX ) THEN |
IF ( .NOT.SEAICEuseIMEX ) THEN |
154 |
C Compute things that do no change during the Newton iteration: |
C Compute things that do no change during the Newton iteration: |
155 |
C sea-surface tilt and wind stress: |
C sea-surface tilt and wind stress: |
156 |
C FORCEX/Y0 - mass*(abA*u/vIceNm1+abB*(u/vIceNm1-u/vIceNm2))/deltaT |
C FORCEX/Y0 - mass*(1.5*u/vIceNm1+0.5*(u/vIceNm1-u/vIceNm2))/deltaT |
157 |
DO J=1-OLy,sNy+OLy |
DO J=1-OLy,sNy+OLy |
158 |
DO I=1-OLx,sNx+OLx |
DO I=1-OLx,sNx+OLx |
159 |
FORCEX(I,J,bi,bj) = FORCEX0(I,J,bi,bj) |
FORCEX(I,J,bi,bj) = FORCEX0(I,J,bi,bj) |