99 |
hFacV = _maskS(i,j,k,bi,bj) - _maskS(i-1,j,k,bi,bj) |
hFacV = _maskS(i,j,k,bi,bj) - _maskS(i-1,j,k,bi,bj) |
100 |
|
|
101 |
e12(I,J,bi,bj) = e12(I,J,bi,bj) |
e12(I,J,bi,bj) = e12(I,J,bi,bj) |
102 |
& + recip_rAz(i,j,bi,bj) * |
& + recip_rAz(i,j,bi,bj) * 2. _d 0 * |
103 |
& ( hFacU * ( _dxC(i,j-1,bi,bj)*uFld(i,j ,1,bi,bj) |
& ( hFacU * ( _dxC(i,j-1,bi,bj)*uFld(i,j ,1,bi,bj) |
104 |
& + _dxC(i,j, bi,bj)*uFld(i,j-1,1,bi,bj) ) |
& + _dxC(i,j, bi,bj)*uFld(i,j-1,1,bi,bj) ) |
105 |
& + hFacV * ( _dyC(i-1,j,bi,bj)*vFld(i ,j,1,bi,bj) |
& + hFacV * ( _dyC(i-1,j,bi,bj)*vFld(i ,j,1,bi,bj) |