[MITgcm-support] how to compute flux to get eta balance

[Martin Losch]

Hi,

in principle, you are right, but to be exact you should use hFacS for the southern and hFacW for the western (and eastern) boundary. In practice, I believe these are the same as hFacC on the boundary (but not 100% positive).

vin=vin+hfacs(i,2,k)*drf(1,1,k)*dxg(i,1)*masks(i,2,k)*vsouth(i,k);

uout=uout+hfacw(174,j,k)*drf(1,1,k)*dyg(174,j)*maskw(174,j,k)*ueast(j,k);