/[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.27 by gforget, Mon Oct 20 03:20:57 2014 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "SEAICE_OPTIONS.h"  #include "SEAICE_OPTIONS.h"
5    #ifdef ALLOW_AUTODIFF
6    # include "AUTODIFF_OPTIONS.h"
7    #endif
8    
9  C--  File seaice_jfnk.F: seaice jfnk dynamical solver S/R:  C--  File seaice_jfnk.F: seaice jfnk dynamical solver S/R:
10  C--   Contents  C--   Contents
# Line 150  C     cycle ice velocities Line 153  C     cycle ice velocities
153            vIceNm1(I,J,bi,bj) = vIce(I,J,bi,bj)            vIceNm1(I,J,bi,bj) = vIce(I,J,bi,bj)
154           ENDDO           ENDDO
155          ENDDO          ENDDO
156          IF ( .NOT.SEAICEuseIMEX ) THEN  C     As long as IMEX is not properly implemented leave this commented out
157    CML        IF ( .NOT.SEAICEuseIMEX ) THEN
158  C     Compute things that do no change during the Newton iteration:  C     Compute things that do no change during the Newton iteration:
159  C     sea-surface tilt and wind stress:  C     sea-surface tilt and wind stress:
160  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 166  C     FORCEX/Y0 - mass*(1.5*u/vIceNm1+0.
166       &         + seaiceMassV(I,J,bi,bj)*dvIcNm1(I,J,bi,bj)*recip_deltaT       &         + seaiceMassV(I,J,bi,bj)*dvIcNm1(I,J,bi,bj)*recip_deltaT
167           ENDDO           ENDDO
168          ENDDO          ENDDO
169          ENDIF  CML        ENDIF
170         ENDDO         ENDDO
171        ENDDO        ENDDO
172  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.27
  ViewVC Help
Powered by ViewVC 1.1.22