| 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) |