[MITgcm-support] open boundary in ocean modeling

[Matthew Mazloff]

Hello

The vertical velocity at the bottom boundary is 0, so it isn’t even written out. The W(k=NR) that you are looking at is the vertical velocity at the top of bottom cell. It is a combination of horizontal volume divergence and bottom form stress.

In short, U and V are on RC. But W is on RF. See 
http://mitgcm.org/sealion/online_documents/node48.html <http://mitgcm.org/sealion/online_documents/node48.html>

-Matt

 
> On Nov 7, 2017, at 9:20 AM, Yangxin He <y67he at uwaterloo.ca> wrote:
> 
> Hi,
> 
> I have a question about the open boundary package obcs in ocean modeling.
> 
> I am currently sending internal waves forced at the left/west boundary into the domain, which i coded up in obcs_calc.F and it seems working fine. 
> 
> I have added a bathymetry at the bottom/south boundary, to do this, I initialized a file through bathyFile in data and also I have set no_slip_bottom=.TRUE.. So the point is I want to have a closed boundary/wall at the bottom/south boundary.
> But the problem is that, with my current setting, the vertical velocity at the bottom boundary is not 0, which is not right.
> 
> How do you suggest I fix this problem?
> 
> Thanks a lot.
> 
> Yangxin
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org <mailto:MITgcm-support at mitgcm.org>
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support <http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20171107/6a43a8a9/attachment.html>