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 |