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

Diff of /MITgcm/pkg/mom_common/mom_calc_tension.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 52  CEOP Line 52  CEOP
52  C Calculate tension of horizontal flow (ignoring lopping factors)  C Calculate tension of horizontal flow (ignoring lopping factors)
53  C *NOTE* that masking is implicit in the contents of the (u,v) fields.  C *NOTE* that masking is implicit in the contents of the (u,v) fields.
54          tension(i,j)=          tension(i,j)=
55       &    (dyf(i,j,bi,bj)*recip_dxf(i,j,bi,bj))       &    ( dyg( i , j ,bi,bj)*uFld( i , j )
56       &   *( uFld(i+1, j )*recip_dyg(i+1, j ,bi,bj)       &     -dyg(i-1, j ,bi,bj)*uFld(i-1, j )
57       &     -uFld( i , j )*recip_dyg( i , j ,bi,bj) )       &     -dxg( i , j ,bi,bj)*vFld( i , j )
58       &   -(dxf(i,j,bi,bj)*recip_dyf(i,j,bi,bj))       &     +dxg( i ,j-1,bi,bj)*vFld( i ,j-1) )*recip_rAz(i,j,bi,bj)
59       &   *( vFld( i ,j+1)*recip_dxg( i ,j+1,bi,bj)  c    &    (dyf(i,j,bi,bj)*recip_dxf(i,j,bi,bj))
60       &     -vFld( i , j )*recip_dxg( i , j ,bi,bj) )  c    &   *( uFld(i+1, j )*recip_dyg(i+1, j ,bi,bj)
61    c    &     -uFld( i , j )*recip_dyg( i , j ,bi,bj) )
62    c    &   -(dxf(i,j,bi,bj)*recip_dyf(i,j,bi,bj))
63    c    &   *( vFld( i ,j+1)*recip_dxg( i ,j+1,bi,bj)
64    c    &     -vFld( i , j )*recip_dxg( i , j ,bi,bj) )
65    
66         ENDDO         ENDDO
67        ENDDO        ENDDO
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22