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 |
& ( dyg( i , j ,bi,bj)*uFld( i , j ) |
& ( dyg(i+1, j ,bi,bj)*uFld(i+1, j ) |
56 |
& -dyg(i-1, j ,bi,bj)*uFld(i-1, j ) |
& -dyg( i , j ,bi,bj)*uFld( i , j ) |
57 |
& -dxg( i , j ,bi,bj)*vFld( i , j ) |
& -dxg( i ,j+1,bi,bj)*vFld( i ,j+1) |
58 |
& +dxg( i ,j-1,bi,bj)*vFld( i ,j-1) )*recip_rAz(i,j,bi,bj) |
& +dxg( i , j ,bi,bj)*vFld( i , j ) )*recip_rA(i,j,bi,bj) |
59 |
c & (dyf(i,j,bi,bj)*recip_dxf(i,j,bi,bj)) |
c & (dyf(i,j,bi,bj)*recip_dxf(i,j,bi,bj)) |
60 |
c & *( uFld(i+1, j )*recip_dyg(i+1, 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) ) |
c & -uFld( i , j )*recip_dyg( i , j ,bi,bj) ) |