[MITgcm-support] hfac question

[David Hebert]

Hi all,

I am attempting to balance barotropic inflow/outflow on the boundaries 
of a 2D simulation. I ran the model one timestep to obtain the hFac's in 
order to more accurately determine the depth the model uses. The result 
from Depth.data at one point is 100m.

In order to balance U, according to manual, I would use hFacW. I get a 
very different result if I do sum(drF * hFacW). Results from matlab (for 
2D hFac's)...

sum(drF .* hfacW(1,:)) = 90
sum(drF .* hfacC(1,:)) = 100

Seems I don't quite understand hfacW. Any ideas on what I did wrong here?

Thanks

David