[MITgcm-support] boundary conditions, in general
Martin Losch
Martin.Losch at awi.de
Fri Dec 22 04:34:12 EST 2006
Hi Jeff,
Helen is right: the periodicity is handled in the "exchange"
routines. Whenever you see a _EXCH_XYZ_RL (uVel,myThid)
(e.g. in model/src/do_fields_blocking_exchanges.F) the tile edges are
updated including periodic (wrap around) boundary conditions. You
supress those by inserting a grid cell row and/or column that is
land, eg. in matlab code:
h=-5000*ones(nx,ny);
h(:,end) = 0;
h(end,:) = 0;
is a closed box with a depth of 5000 meters.
periodicity in the vertical is not supported, there is always
"bottom" at one end and free surface at the other end and the
boundary conditions are tailored to this situation. In particular
there are no "overlaps" in the vertical direction. After all, it's an
ocean model and the ocean is hardly ever periodic in z (o:
Martin
On 21 Dec 2006, at 18:15, helen hill wrote:
> Hi Jeff
>
> I think the model is doubly periodic by default so you may not find
> anything to indicate that the model is being made to be doubly
> periodic.
>
> As to your second question - I'm not sure - maybe someone else knows
>
> Helen
>
> On 12/21/06, Haferman, Jeffrey (Jeff) (CIV) <jlhaferm at nps.edu> wrote:
>
> Hi -
> I'm looking at exp5 supplied with MITgcm (imhomogeneously forced ocean
> convection in a doubly-periodic box).
>
> I've spent quite a bit of time looking through the documentation
> and the
> input/data file, and it's not clear to me what settings are making the
> lateral (x, y) boundary conditions to be periodic. Can someone
> clarify?
>
> Also, is it possible to make the boundary conditions periodic in the
> vertical (z, or in the input code "r") dimension?
>
> Thanks,
> Jeff
>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list