[MITgcm-support] Specify W velocity with OBCS?
Matthew Mazloff
mmazloff at MIT.EDU
Mon Jan 14 14:43:55 EST 2008
Hi David,
It is just a nice feature of the c-grid and obcs implementation that
I really think the W and P in the center point of the obcs does not
effect the interior in the hydrostatic mode. If you are using non-
hydrostatic and have terms like Wxx this may become an issue...and as
the MITgcm currently does not have the capability of prescribing W
you could either
1) carefully choose tangential velocity so W get is what you want
2) implement the code so that W is overwritten after it is solved for
by divergence of tangential velocity
3) use sponge layers to limit effect of whatever is happening right
on the boundary
good luck
Matt
ps> Not sure what you mean by a "swirl", maybe someone else can help
you with that.
On Jan 14, 2008, at 2:17 PM, David Hebert wrote:
> Hi Matt,
>
> Thanks for your response. I can see that a point on a boundary
> wouldn't be solved for since it is specified, but wouldn't it
> affect the adjacent interior point somehow? I probably should note
> that I am using the model in non-hydrostatic mode.
>
> My previous question stems from wanting to specify a jet on the
> boundary. Implementing a "Top Hat" jet on the west boundary using
> OBWuFile seems to work. If I wanted to model a jet with swirl I was
> thinking I would specify all three velocity components on the
> boundary. A Google search for 'MITgcm OBWwFile' gives a few hits,
> so I thought it would be possible to specify vertical velocity on
> the boundary. Perhaps this was possible in the past? If not, it
> looks like something I will have to add to the code.
>
> As m. schaferkotter noted, the namelist in MITgcm/pkg/obcs/
> obcs_readparms.F does not have any W velocity boundary files.
>
> Thanks,
>
> David
>
> Matthew Mazloff wrote:
>> Hi David,
>>
>> I do not believe W (and P) on the open boundary comes into the
>> dynamics. If W on this point is solved for (and thus given as
>> output) it will be because there is some divergence in the
>> prescribed tangential velocity....but this doesn't effect the
>> interior.
>>
>> -Matt
>>
>>
>> On Jan 14, 2008, at 11:06 AM, David Hebert wrote:
>>
>>> All,
>>>
>>> Is there a way so specify W velocity using OBCS and EXF? When I
>>> specify 'OBWwFile = w.bin' in data.obcs, I get an error saying it
>>> is not found in namelist:
>>>
>>> PGFIO-F-239/namelist read/unit=11/entity name is not member of
>>> group.
>>> File name = /tmp/FTNwcP5jg7mUsXg8 formatted, sequential
>>> access record = 17
>>> In source file obcs_readparms.f, at line number 1989
>>>
>>> Thanks,
>>>
>>> David
>>>
>>>
>>> _______________________________________________
>>> 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