/[MITgcm]/MITgcm/pkg/seaice/seaice_jfnk.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_jfnk.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.24 by mlosch, Tue Feb 4 18:30:31 2014 UTC revision 1.26 by mlosch, Thu Mar 20 09:24:49 2014 UTC
# Line 150  C     cycle ice velocities Line 150  C     cycle ice velocities
150            vIceNm1(I,J,bi,bj) = vIce(I,J,bi,bj)            vIceNm1(I,J,bi,bj) = vIce(I,J,bi,bj)
151           ENDDO           ENDDO
152          ENDDO          ENDDO
153          IF ( .NOT.SEAICEuseIMEX ) THEN  C     As long as IMEX is not properly implemented leave this commented out
154    CML        IF ( .NOT.SEAICEuseIMEX ) THEN
155  C     Compute things that do no change during the Newton iteration:  C     Compute things that do no change during the Newton iteration:
156  C     sea-surface tilt and wind stress:  C     sea-surface tilt and wind stress:
157  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
158          DO J=1-OLy,sNy+OLy          DO J=1-OLy,sNy+OLy
159           DO I=1-OLx,sNx+OLx           DO I=1-OLx,sNx+OLx
160            FORCEX(I,J,bi,bj) = FORCEX0(I,J,bi,bj)            FORCEX(I,J,bi,bj) = FORCEX0(I,J,bi,bj)
# Line 162  C     FORCEX/Y0 - mass*(abA*u/vIceNm1+ab Line 163  C     FORCEX/Y0 - mass*(abA*u/vIceNm1+ab
163       &         + seaiceMassV(I,J,bi,bj)*dvIcNm1(I,J,bi,bj)*recip_deltaT       &         + seaiceMassV(I,J,bi,bj)*dvIcNm1(I,J,bi,bj)*recip_deltaT
164           ENDDO           ENDDO
165          ENDDO          ENDDO
166          ENDIF  CML        ENDIF
167         ENDDO         ENDDO
168        ENDDO        ENDDO
169  C     Start nonlinear Newton iteration: outer loop iteration  C     Start nonlinear Newton iteration: outer loop iteration
Legend:
Removed from v.1.24  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.22