[MITgcm-support] OBCs over topography!

Van Thinh Nguyen vtnguyen at moisie.math.uwaterloo.ca
Thu Oct 20 12:13:24 EDT 2005


It worked, I did a test already.
I wonder it may ignore cells at North & South where are defined as side
walls.

Vth
-----
On Thu, 20 Oct 2005, Matthew Mazloff wrote:

> Actually I would not put OB_Ieast=0, 78*-1, 0,  as I don't think this
> will work properly...i'm not sure what this will do?....did you find it
> does work?
>
> I would use OB_Ieast=80*-1
>
> -Matt
>
>
> Matthew Mazloff wrote:
>
> > Hi Van Thinh,
> >
> > I actually think either would be fine.  I believe the masks will
> > over-rule the obcs-prescribe so that if you prescribe a T,S,U,V where
> > there is land the model will not be effected.
> > Please let me know if you find otherwise.
> >
> > -Matt
> >
> >
> >
> > Van Thinh Nguyen wrote:
> >
> >>Hi Matt,
> >>
> >> Instead of:
> >>     OB_Ieast=80*-1
> >>     OB_Iwest=80*1
> >>
> >> I should write:  OB_Ieast=0, 78*-1, 0
> >>                  OB_Iwest=0, 78*1, 0
> >>
> >> Is that what you meaned?
> >>
> >> Thanks
> >>
> >> Van Thinh
> >>
> >>On Thu, 20 Oct 2005, Matthew Mazloff wrote:
> >>
> >>
> >>
> >>>Hi Van Thinh!
> >>>
> >>>You can give all cells and the masks will take care of it.  (e.g.
> >>>OB_Ieast=80*-1,OB_Iwest=80*1 is fine)
> >>>
> >>>Other things I have learned:
> >>>1) You should make the topogrophy flat right at the obcs.  (e.g. make
> >>>the height at (i,j) = (i+1,j) for the western boundary)
> >>>2) Do not use the CD scheme with obcs
> >>>
> >>>-Matt
> >>>
> >>>
> >>>Van Thinh Nguyen wrote:
> >>>
> >>>
> >>>
> >>>>Hi
> >>>>
> >>>>I'm confused how to set up the boundary condition at a OB where a
> >>>>topography intersects with the grid block (please see my figure enclosed).
> >>>>I would like to simulate internal waves due to a tidal flow over a
> >>>>depression Gaussian bump in a channel, which has OBCs at West and East,
> >>>>North & South are closed sides. It looks like the geometry in exp4, the
> >>>>difference is only the bump (my case is depression bump!).
> >>>>Model is 3D, but the Fig. enclosed is a (X,Z) section at middle of bump.
> >>>>
> >>>>At West/East OBs, topography intersects grid block at cell 55 (Z direction
> >>>>with 80 cells). It means:
> >>>>
> >>>>- from free surface to cell 55 is OB
> >>>>- from cell 55 to the bottom is closed
> >>>>
> >>>>In this case, should I define cells (1-55) are OBCs in data.obcs or I can
> >>>>define all cells at West/East are OBCs (with OB_Ieast=80*-1,
> >>>>OB_Iwest=80*1)
> >>>>
> >>>>
> >>>>Could someone please give a hint for this?
> >>>>
> >>>>Thanks so much!
> >>>>
> >>>>Van Thinh
> >>>>
> >>>>------------------------------------------------------------------------
> >>>>
> >>>>_______________________________________________
> >>>>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
> >>
> >>
> >------------------------------------------------------------------------
> >
> >_______________________________________________
> >MITgcm-support mailing list
> >MITgcm-support at mitgcm.org
> >http://mitgcm.org/mailman/listinfo/mitgcm-support
> >
> >
>



More information about the MITgcm-support mailing list