/[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.25 by mlosch, Fri Feb 7 14:27:21 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*(1.5*u/vIceNm1+0.5*(u/vIceNm1-u/vIceNm2))/deltaT  C     FORCEX/Y0 - mass*(1.5*u/vIceNm1+0.5*(u/vIceNm1-u/vIceNm2))/deltaT
# Line 162  C     FORCEX/Y0 - mass*(1.5*u/vIceNm1+0. Line 163  C     FORCEX/Y0 - mass*(1.5*u/vIceNm1+0.
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.25  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.22