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

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

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

revision 1.4 by mlosch, Tue Apr 24 02:59:31 2007 UTC revision 1.5 by mlosch, Tue May 15 14:32:55 2007 UTC
# Line 52  C  hFacU, hFacV - determine the no-slip Line 52  C  hFacU, hFacV - determine the no-slip
52  C  C
53        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
54         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
55          DO j=1-Oly+1,sNy+Oly-1          DO j=1-Oly,sNy+Oly-1
56           DO i=1-Olx+1,sNx+Olx-1           DO i=1-Olx,sNx+Olx-1
57  C NOW EVALUATE STRAIN RATES  C     evaluate strain rates
58            e11(I,J,bi,bj)= _recip_dxF(I,J,bi,bj) *            e11(I,J,bi,bj)= _recip_dxF(I,J,bi,bj) *
59       &         (uFld(I+1,J,bi,bj)-uFld(I,J,bi,bj))       &         (uFld(I+1,J,bi,bj)-uFld(I,J,bi,bj))
60       &         -HALF*       &         -HALF*
# Line 63  C NOW EVALUATE STRAIN RATES Line 63  C NOW EVALUATE STRAIN RATES
63            e22(I,J,bi,bj)= _recip_dyF(I,J,bi,bj) *            e22(I,J,bi,bj)= _recip_dyF(I,J,bi,bj) *
64       &         (vFld(I,J+1,bi,bj)-vFld(I,J,bi,bj))       &         (vFld(I,J+1,bi,bj)-vFld(I,J,bi,bj))
65  C     one metric term is missing  C     one metric term is missing
66             ENDDO
67            ENDDO
68            DO j=1-Oly+1,sNy+Oly
69             DO i=1-Olx+1,sNx+Olx
70            e12(I,J,bi,bj)=HALF*(            e12(I,J,bi,bj)=HALF*(
71       &         (uFld(I  ,J  ,bi,bj) * _dxC(I  ,J  ,bi,bj)       &         (uFld(I  ,J  ,bi,bj) * _dxC(I  ,J  ,bi,bj)
72       &         -uFld(I  ,J-1,bi,bj) * _dxC(I  ,J-1,bi,bj)       &         -uFld(I  ,J-1,bi,bj) * _dxC(I  ,J-1,bi,bj)
# Line 81  C     one metric term is missing Line 85  C     one metric term is missing
85          ENDDO          ENDDO
86          IF ( SEAICE_no_slip ) THEN          IF ( SEAICE_no_slip ) THEN
87  C     no slip boundary conditions apply only to e12  C     no slip boundary conditions apply only to e12
88           DO j=1-Oly+1,sNy+Oly-1           DO j=1-Oly+1,sNy+Oly
89            DO i=1-Olx+1,sNx+Olx-1            DO i=1-Olx+1,sNx+Olx
90             hFacU = _maskW(i,j,k,bi,bj) - _maskW(i,j-1,k,bi,bj)             hFacU = _maskW(i,j,k,bi,bj) - _maskW(i,j-1,k,bi,bj)
91             hFacV = _maskS(i,j,k,bi,bj) - _maskS(i-1,j,k,bi,bj)             hFacV = _maskS(i,j,k,bi,bj) - _maskS(i-1,j,k,bi,bj)
92    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22