/[MITgcm]/MITgcm/pkg/mom_common/mom_calc_strain.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_common/mom_calc_strain.F

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

revision 1.1 by adcroft, Fri May 14 17:43:11 2004 UTC revision 1.2 by adcroft, Mon May 24 18:44:27 2004 UTC
# Line 55  CEOP Line 55  CEOP
55    
56  C       Strain of horizontal flow field (ignoring lopping factors)  C       Strain of horizontal flow field (ignoring lopping factors)
57          strain(I,J)=          strain(I,J)=
58       &    dyu(I,J,bi,bj)*recip_dxv(I,J,bi,bj)*(       &    ( dyc( i , j ,bi,bj)*vFld( i , j )
59       &              vFld( I , J )*recip_dyc( I , J ,bi,bj)       &     -dyc(i-1, j ,bi,bj)*vFld(i-1, j )
60       &             -vFld(I-1, J )*recip_dyc(I-1, J ,bi,bj) )       &     +dxc( i , j ,bi,bj)*uFld( i , j )
61       &   +dxv(I,J,bi,bj)*recip_dyu(I,J,bi,bj)*(       &     -dxc( i ,j-1,bi,bj)*uFld( i ,j-1) )*recip_rAz(i,j,bi,bj)
62       &             +uFld( I , J )*recip_dxc( I , J ,bi,bj)  c    &    dyu(I,J,bi,bj)*recip_dxv(I,J,bi,bj)*(
63       &             -uFld( I ,J-1)*recip_dxc( I ,J-1,bi,bj) )  c    &              vFld( I , J )*recip_dyc( I , J ,bi,bj)
64    c    &             -vFld(I-1, J )*recip_dyc(I-1, J ,bi,bj) )
65    c    &   +dxv(I,J,bi,bj)*recip_dyu(I,J,bi,bj)*(
66    c    &             +uFld( I , J )*recip_dxc( I , J ,bi,bj)
67    c    &             -uFld( I ,J-1)*recip_dxc( I ,J-1,bi,bj) )
68    
69  C       Set strain to zero on boundaries (free-slip)  C       Set strain to zero on boundaries (free-slip)
70          IF (hFacZ(i,j).EQ.0.) THEN          IF (hFacZ(i,j).EQ.0.) THEN
# Line 72  C       Set strain to zero on boundaries Line 76  C       Set strain to zero on boundaries
76                    
77  C     Special stuff for Cubed Sphere  C     Special stuff for Cubed Sphere
78        IF (useCubedSphereExchange) THEN        IF (useCubedSphereExchange) THEN
79          I=1         STOP 'S/R MOM_CALC_STRAIN: We should not use strain on the cube!'
         J=1  
         strain(I,J)=0.  
         I=sNx+1  
         J=1  
         strain(I,J)=0.  
         I=1  
         J=sNy+1  
         strain(I,J)=0.  
         I=sNx+1  
         J=sNy+1  
         strain(I,J)=0.  
80        ENDIF        ENDIF
81    
82        RETURN        RETURN
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22