[MITgcm-support] oscillating mean flow
Klymak Jody
jklymak at uvic.ca
Thu Jun 10 11:53:14 EDT 2010
Hi Mark,
On 10-Jun-10, at 8:33 AM, Mark Inall wrote:
> Jody, I tried adding a constant mean flow at the boundaries. This
> worked fine with f=0: the non-zero U BCs propagated into the domain
> at the barotropic phase speed. But with rotation, the imposed U-
> momentum wants to rotate, and then there's a mismatch between the N
> and S boundaries (at least that's what I think happens) ... that's
> why I wondered about imposing a cross slice pressure gradient (or
> moving to a channel with walls) - but I'm not sure how to do that ....
Oh yeah, you are right! Shouldn't post before coffee...
If your proscribed boundary flow is U_o, then you'd get high in the
west and low in the east and a southward BT flow. An eastward bottom
Ekman layer would allow a steady state and the model is behaving as it
should. However, you wouldn't have a mean eastward flow in the
interior.
I haven't played with it too much, but you should be able to set up a
body force to mimic a N/S pressure gradient in the interior that would
set up a eastward mean flow. Might be hard to match with your E/W
boundaries though. I think it can be done in external_forcing.F where
you would modify gV in your case. There maybe a better way to do it,
and to be honest I'm not sure of the units, but I think the docs
explain it.
Anyways, thats what I'd do. A steady gV to the north will give you an
eastwards geostrophically balanced flow. Cheaper than making a
channel, but thats your other alternative.
Cheers, Jody
>
> Mark.
>
>>>> Samar Khatiwala <spk at ldeo.columbia.edu> 10 June 2010 16:00 >>>
> Hi Mark,
>
> Yes, this has been done many times, for e.g., see here: http://www.ldeo.columbia.edu/~spk/Papers/KhatiwalaInternalTidesDSR03.pdf
>
> This is a bit dated, but I know others have implemented it in more
> recent versions of the code.
>
> Samar
>
> On Jun 10, 2010, at 10:16 AM, Mark Inall wrote:
>
>> Hi,
>>
>> Does anyone have experience with setting up an oscillating mean
>> flow - either in 2D or in a narrow channel - with rotation?
>>
>> Presently I'm running an oscillating flow over a bump based on the
>> internal wave verification setup, and I want to add a mean flow -
>> is the best way to impose a barotropic pressure gradient across the
>> section and wait for a geostrophic balance??
>>
>> Any thoughts/help much appreciated.
>>
>> Thanks,
>>
>> Mark.
>>
>>
>>
>>
>> Dr Mark Inall
>> Reader in Marine Physics
>> Scottish Association for Marine Science
>> Scottish Marine Institute
>> Oban PA37 1QA
>>
>> E: mei at sams.ac.uk
>> T:+44 (0) 1631 559267
>> M:+44 (0) 7795 573475
>> F:+44 (0) 1631 559001
>>
>> +++++++++++++++++++++++++++++++
>> The Scottish Association for Marine Science (SAMS) is registered in
>> Scotland as a Company Limited by Guarantee (SC009292) and is a
>> registered charity (9206). There are 2 wholly owned subsidiary
>> companies: SAMS Research Services Ltd a Limited Company (SC224404)
>> and the European Centre for Marine Biotechnology, a Company Limited
>> By Guarantee (SC205318). All Companies in the group are registered
>> in Scotland and share a registered office at Scottish Marine
>> Institute, Oban Argyll PA37 1QA.
>>
>> The content of this message may contain personal views which are
>> not the views of SAMS unless specifically stated.
>>
>> Please note that all email traffic is monitored for purposes of
>> security and spam filtering. As such individual emails may be
>> examined in more detail.
>> +++++++++++++++++++++++++++++++
>>
>>
>>
>> _______________________________________________
>> 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