92 |
|
|
93 |
e12(I,J,bi,bj) = e12(I,J,bi,bj) |
e12(I,J,bi,bj) = e12(I,J,bi,bj) |
94 |
& + recip_rAz(i,j,bi,bj) * |
& + recip_rAz(i,j,bi,bj) * |
95 |
& ( hFacU * ( _dxC(i,j ,bi,bj)*uFld(i,j ,bi,bj) |
& ( hFacU * ( _dxC(i,j-1,bi,bj)*uFld(i,j ,bi,bj) |
96 |
& + _dxC(i,j-1,bi,bj)*uFld(i,j-1,bi,bj) ) |
& + _dxC(i,j, bi,bj)*uFld(i,j-1,bi,bj) ) |
97 |
& + hFacV * ( _dyC(i ,j,bi,bj)*vFld(i ,j,bi,bj) |
& + hFacV * ( _dyC(i-1,j,bi,bj)*vFld(i ,j,bi,bj) |
98 |
& + _dyC(i-1,j,bi,bj)*vFld(i-1,j,bi,bj) ) ) |
& + _dyC(i, j,bi,bj)*vFld(i-1,j,bi,bj) ) ) |
99 |
& - hFacU |
& - hFacU |
100 |
& * 0.25 _d 0 * (uFld(I,J,bi,bj)+uFld(I ,J-1,bi,bj)) |
& * 0.25 _d 0 * (uFld(I,J,bi,bj)+uFld(I ,J-1,bi,bj)) |
101 |
& * ( _tanPhiAtU(I,J,bi,bj) + _tanPhiAtU(I,J-1,bi,bj) ) |
& * ( _tanPhiAtU(I,J,bi,bj) + _tanPhiAtU(I,J-1,bi,bj) ) |