[MITgcm-support] Implementing a sponge layer in the x-position

[Matthew Mazloff]

The obcs package should do this too.

Matt

On Mon, Sep 16, 2019 at 5:45 AM Edward W Doddridge <ewd at mit.edu> wrote:

> Dear Kevin,
>
> Is there as reason that you can’t use the RBCS package?
> https://mitgcm.readthedocs.io/en/latest/phys_pkgs/rbcs.html
>
> Cheers,
> Ed
>
>
> ------------------------------
> *Edward Doddridge*
> *Postdoctoral researcher*
> *Earth, Atmospheric and Planetary Sciences*
> *Massachusetts Institute of Technology*
>
> www.doddridge.me
>
> On 16 Sep 2019, at 04:48, Kevin Ha <kevin.ha.pro at gmail.com> wrote:
>
> Dear community,
>
> I am trying to implement a sponge layer denoted a(x) varying between 0 and
> 1 in the x-direction so that my velocity forcing relaxes at the targeted
> value (Uf,Vf)=(0.025,0) m/s at the two edges of my domain after a certain
> time (defined as some multiples of 3600 s) and (Uf,Vf)=(0,0) i.e. no
> forcing in the area of interest. My domain has an horizontal extent of Lx =
> sNx*NPx*dx = 50*40*30 m = 60 km long and I have added the mask function by
> modifying the right hand side of the Navier Stokes equation gU, gV and gT
> as followed (see lines implemented in external_forcing.F) to get an
> exponential decay after a certain time. It seems however that the code is
> not taking into account the sponge layer at all and keeps resolving Navier
> Stokes without any mask restrictions (which results in oscillations of U,V
> and theta in the longer times instead of the expected relaxation). Any idea
> about this issue ?
>
> Thank you in advance
>
> Kind regards
>
> Kevin.
>
>
> <packages.conf><gendata.m><eedata><data><external_forcing.F><FFIELDS.h>
> <CPP_OPTIONS.h><MOM_FLUXFORM_OPTIONS.h><ini_fields.F><SIZE.h>
> _______________________________________________
>
>
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
>
> _______________________________________________
> 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/20190916/9a8681ac/attachment-0001.html>