[MITgcm-support] Open boundary conditions for Internal Waves.

[Jeremy Miller]

Hi Jody,
Thanks for coming back to me.
To begin with, I want to understand (in detail) what are the default OBCs
being imposed when I run the internal wave verification with default
settings.
In particular, the docs in section 8.3.1.4.1 say that a detailed
calculation is used to set OBCs for the IW verification in "obcs_calc.f". I
want to understand what those same OBCs are.

Taking your example of  a wide sponge, I'm not clear on the notation

L ~ mode -1*wavelength

But nonetheless, how would I set that using a matlab or python snippet?
Best
Jeremy




On Tue, 25 May 2021 at 17:16, Jody Klymak <jklymak at uvic.ca> wrote:

> Hi Jeremy,
>
> It would be helpful to know what you are trying to simulate.  Overall, I
> have found internal waves are best set using a wide sponge (L\sim mode-1
> wavelength) either with obcs or with rbcs packages.
>
> Cheers,   Jody
>
> On 25 May 2021, at 00:14, Jeremy Miller <jeremysharonmiller at gmail.com>
> wrote:
>
> Dear MITgcm forum,
>
> I'm sorry for sending multiple emails. The above error I did resolve. I
> would still be happy to learn what the default OBCs are, and how to custom
> create OBCs.
> Best
> Jeremy
>
>
> On Tue, 25 May 2021 at 08:50, Jeremy Miller <jeremysharonmiller at gmail.com>
> wrote:
>
>> Dear MITgcm forum,
>>
>> Hope everyone's well. Following yesterday's email and reading Section
>> 8.3.1.4.1 of the docs I attempted to change the file data.obcs myself,
>> because I wanted to move from the default grid to Nx ✕ Ny = 120 ✕ 60 .
>> Accordingly I changed data.obcs to
>>
>> # Open-boundaries
>>  &OBCS_PARM01
>>  OB_Ieast = 60*-1,
>>  OB_Iwest = 60*1,
>>  OB_Jnorth = 120*-1,
>>  OB_Jsouth = 120*1,
>>  &
>>
>> After running I get this error:
>>
>> At line 2509 of file obcs_readparms.f (unit = 11, file =
>> 'scratch1.000000000')
>> Fortran runtime error: Cannot match namelist object name c
>>
>> Please, could you advise.
>>
>> Best
>> Jeremy
>>
>> On Mon, 24 May 2021 at 13:27, Jeremy Miller <jeremysharonmiller at gmail.com>
>> wrote:
>>
>>> Dear MITgcm forum,
>>>
>>> Hope all are well. I want to input open boundary conditions (OBCs) for
>>> the internal wave verification. I have a few questions.
>>>
>>> 1. What are the default OBCs ?
>>> 2. How does one change the OBCs? I was reading through obcs_calc.F and I
>>> found it somewhat difficult to understand (not coming with a Fortran
>>> background).
>>> 3. Is it possible to output the OBCs from a short python snippet?
>>>
>>> Best
>>> Jeremy Miller
>>>
>>>
>>>
>>> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
>
> --
> Jody Klymak
> http://ocean-physics.seos.uvic.ca/~jklymak/
>
>
>
>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20210526/0be435f4/attachment.html>