[MITgcm-support] Zonally integrated momentum budget

[Noriyuki Yamamoto]

Dear all.

I'm having a problem when I try to check the zonal momentum budget
that is zonally and vertically integrated.

The model has reached nearly steady state and
zonally and vertically integrated (at the c-grid V points) meridonal
volume transport is almost zero.
However, zonnally and vertically integrated Um_Cori or V (interpolated
on the U points) have nonzero values.
The way I integrate them is like following

do i = 1, imax
do k = 1, kmax
integrated_value(j) = integrated_value(j) + Um_Cori(i,j,k) * dx(i,j,k) *
dz(k) * hFacW(i,j,k)
enddo
enddo

This occurs probably because of eliminations of half grid next to
topography (here hFacW = 0).
Other momentum diagnostic terms such as Um_Advec, Um_Diss or
ADVy_Um, VISCy_Um also seem to have zero values where hFacW = 0.

Is there any way to include the contribution of half grids next to
topography?

I will be grateful for any help.
Noriyuki

-- 
Noriyuki Yamamoto
PhD Student - Physical Oceanography Group 
Division of Earth and Planetary Sciences,
Graduate School of Science, Kyoto University.
Mail:nymmto at kugi.kyoto-u.ac.jp
Tel:+81-75-753-3924