[MITgcm-support] Specify W velocity with OBCS?

David Hebert hebert at nrlssc.navy.mil
Mon Jan 14 14:59:15 EST 2008


Thanks Matt, I appreciate the explanation.

By "swirl" I mean the motion of a jet induced by a rotating propeller. 
The jet would have some angular velocity in addition to an axial velocity.

David

Matthew Mazloff wrote:
> 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
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>




More information about the MITgcm-support mailing list