[MITgcm-support] Re: doubly-reentrant
baylor at MIT.EDU
Thu Mar 31 10:20:25 EST 2005
Yes Joe (and mitgcm-support for posterity),
That's true. As I understand the finite volume part of it, you
figure the flux through every face of a volume centered on the
(temperature, salt, divergence) point.
First, you also don't change tracer inside topography: in
gad_advection, when you figure the flux divergence you get lines like:
& *( af(i+1,j)-af(i,j)
The hFacC (from ini_masks_etc.F) here gives you the topography, i.e.,
in a wall you don't change temp. But, you might still have a flux
There is a mask hfacW that is 0 for blocked western cell faces, 1 for
open cell faces, and a fraction for lopped cell faces. Similarly hfacS
is for southern faces...
You then use this factor to calculate cell face area (which you do
every timestep in calc_common_factors, called from thermodynamics.F to
be sure that the free surface is treated properly):
C-- Calculate tracer cell face open areas
xA(i,j) = _dyG(i,j,bi,bj)
yA(i,j) = _dxG(i,j,bi,bj)
So, now xa has info about the topography built in.
For the advective fluxes, you also make sure topography is built into
uTrans and vTrans (also from calc_common_factors):
C-- Calculate velocity field "volume transports" through
C-- tracer cell faces.
uTrans(i,j) = uVel(i,j,k,bi,bj)*xA(i,j)
vTrans(i,j) = vVel(i,j,k,bi,bj)*yA(i,j)
So for diffusive fluxes, you find in gad_diff_x that xA comes in to set
dfx(i,j) = -diffKh*xA(i,j)
And for all of the various flux limiting advection schemes
(gad_*_adv_*.F) you find that the fluxes are always proportional in the
end to a uFlux which already has the b.c. built in...or that there is
an xA to make it work.
Similar detective work looking for hFacUp will tell you when vertical
fluxes are blocked by topography...
On Mar 31, 2005, at 8:51 AM, Joseph LaCasce wrote:
> Hey Baylor---
> What does that mean, exactly? If the boundary lies, say, in the
> middle of a channel, you still have to impose no flux on either
> side of it, don't you?
> J.H. LaCasce, Phone: +47 22 96 33 11; FAX: +47 22 69 63 55
> Norwegian Meteorological Institute, P.O. Box 43 Blindern, 0313 Oslo
> or (UPS) Niels Henrik Abelsvei 40, 0371 Oslo, Norway
More information about the MITgcm-support