[MITgcm-support] obcs balance
elisabei at yahoo.com.cn
Thu Dec 18 21:23:45 EST 2008
I had set exactConserv=.true. in data (PARM01),but it seems that the result had not a major improvement.
Also ,I had the tests as you showed.The sea level hardly raised if all boundaries had been close or did not prescript the velocities normal to the boundaries.Meanwhile,I had another test.If only the west boundary velocities been prescribed,the sea level hardly raised.But, If only the south boundary velocities been prescribed,the sea level evidently raised.I think the larger flux of volume caused the sea level raise,beacause the south boundary flux of volume is larger 15.6% than the west boundary flux of volume.As the exp4 the open boundary flux of volume is small,so the sea level hardly raised after ten years integrated.
Also ,I have a question: which routine apply the parameter "obcs balance" or which routine balance the boundary flux of volume ?
I've had a quick look through your files. There are a few things you can
try to narrow down the problem:
1. set exactConserv=.true. in data (PARM01); this should improve the
conservation of volume.
2. close all of you boundaries (e.g. by creating a topography with a wall at
i=1 and j=1) and see what happens
3. with the open boundaries turned on, do not prescript the velocities normal
to the boundaries (the default should be zero, so that there should not be any
flux of volume through the boundaries).
PS. can you try to reproduce your problem with a smaller domain for faster
integrations (say a 10x20 or so box). that would make debugging simpler for you
and then we could try this on one of our machine and by that help you much
On 16 Dec 2008, at 10:53, 名 徐 wrote:
> Thank you for your answer,I have sent the all files you need on enclosure.
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
MITgcm-support mailing list
MITgcm-support at mitgcm.org
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the MITgcm-support