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

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

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

revision 1.17 by jmc, Mon May 14 14:29:49 2007 UTC revision 1.18 by mlosch, Tue May 15 14:32:56 2007 UTC
# Line 89  C     anything else we have to do it her Line 89  C     anything else we have to do it her
89             ENDDO             ENDDO
90            ENDDO            ENDDO
91    
92            DO j=1-Oly+1,sNy+Oly-1            DO j=0,sNy
93             DO i=1-Olx+1,sNx+Olx-1             DO i=0,sNx
94              eplus = e11(I,J,bi,bj) + e22(I,J,bi,bj)              eplus = e11(I,J,bi,bj) + e22(I,J,bi,bj)
95              eminus= e11(I,J,bi,bj) - e22(I,J,bi,bj)              eminus= e11(I,J,bi,bj) - e22(I,J,bi,bj)
96              sig11(I,J) = zeta(I,J,bi,bj)*eplus + eta(I,J,bi,bj)*eminus              sig11(I,J) = zeta(I,J,bi,bj)*eplus + eta(I,J,bi,bj)*eminus
97       &           - 0.5 _d 0 * PRESS(I,J,bi,bj)       &           - 0.5 _d 0 * PRESS(I,J,bi,bj)
98              sig22(I,J) = zeta(I,J,bi,bj)*eplus - eta(I,J,bi,bj)*eminus              sig22(I,J) = zeta(I,J,bi,bj)*eplus - eta(I,J,bi,bj)*eminus
99       &           - 0.5 _d 0 * PRESS(I,J,bi,bj)       &           - 0.5 _d 0 * PRESS(I,J,bi,bj)
100               ENDDO
101              ENDDO
102    
103              DO j=1,sNy+1
104               DO i=1,sNx+1
105              sig12(I,J) = 2. _d 0 * e12(I,J,bi,bj) *              sig12(I,J) = 2. _d 0 * e12(I,J,bi,bj) *
106       &           ( eta(I,J  ,bi,bj) + eta(I-1,J  ,bi,bj)       &           ( eta(I,J  ,bi,bj) + eta(I-1,J  ,bi,bj)
107       &           + eta(I,J-1,bi,bj) + eta(I-1,J-1,bi,bj) )       &           + eta(I,J-1,bi,bj) + eta(I-1,J-1,bi,bj) )
# Line 181  C     ice cover Line 186  C     ice cover
186         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
187          DO j=1,sNy          DO j=1,sNy
188           DO i=1,sNx           DO i=1,sNx
189            fuIceLoc=HALF*( DWATN(I,J,bi,bj)+DWATN(I,J+1,bi,bj) )*            fuIceLoc=HALF*( DWATN(I,J,bi,bj)+DWATN(I-1,J,bi,bj) )*
190       &         COSWAT *       &         COSWAT *
191       &         ( UICE(I,J,1,bi,bj)-GWATX(I,J,bi,bj) )       &         ( UICE(I,J,1,bi,bj)-GWATX(I,J,bi,bj) )
192       &         - SIGN(SINWAT, _fCori(I,J,bi,bj)) * 0.5 _d 0 *       &         - SIGN(SINWAT, _fCori(I,J,bi,bj)) * 0.5 _d 0 *
# Line 192  C     ice cover Line 197  C     ice cover
197       &         0.5 _d 0*(vIce(I-1,J  ,1,bi,bj)-GWATY(I-1,J  ,bi,bj)       &         0.5 _d 0*(vIce(I-1,J  ,1,bi,bj)-GWATY(I-1,J  ,bi,bj)
198       &                  +vIce(I-1,J+1,1,bi,bj)-GWATY(I-1,J+1,bi,bj))       &                  +vIce(I-1,J+1,1,bi,bj)-GWATY(I-1,J+1,bi,bj))
199       &         )       &         )
200            fvIceLoc=HALF*( DWATN(I,J,bi,bj)+DWATN(I+1,J,bi,bj) )*            fvIceLoc=HALF*( DWATN(I,J,bi,bj)+DWATN(I,J-1,bi,bj) )*
201       &         COSWAT *       &         COSWAT *
202       &         ( VICE(I,J,1,bi,bj)-GWATY(I,J,bi,bj) )       &         ( VICE(I,J,1,bi,bj)-GWATY(I,J,bi,bj) )
203       &         + SIGN(SINWAT,  _fCori(I,J,bi,bj)) * 0.5 _d 0 *       &         + SIGN(SINWAT,  _fCori(I,J,bi,bj)) * 0.5 _d 0 *
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22