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

[Kevin Ha]

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.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190916/1fea5d62/attachment-0001.html>
-------------- next part --------------
# $Header: /u/gcmpack/MITgcm/verification/internal_wave/code/packages.conf,v 1.5 2007/02/19 23:24:40 jmc Exp $
# $Name:  $

gfd
#obcs
#timeave
#mnc
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gendata.m
Type: application/vnd.wolfram.mathematica.package
Size: 1117 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190916/1fea5d62/attachment-0001.m>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: eedata
Type: application/octet-stream
Size: 254 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190916/1fea5d62/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: data
Type: application/octet-stream
Size: 1482 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190916/1fea5d62/attachment-0003.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: external_forcing.F
Type: text/x-fortran
Size: 21013 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190916/1fea5d62/attachment-0006.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: FFIELDS.h
Type: text/x-chdr
Size: 10162 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190916/1fea5d62/attachment-0007.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: CPP_OPTIONS.h
Type: text/x-chdr
Size: 1804 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190916/1fea5d62/attachment-0008.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MOM_FLUXFORM_OPTIONS.h
Type: text/x-chdr
Size: 651 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190916/1fea5d62/attachment-0009.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ini_fields.F
Type: text/x-fortran
Size: 1850 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190916/1fea5d62/attachment-0010.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: SIZE.h
Type: text/x-chdr
Size: 2165 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190916/1fea5d62/attachment-0011.bin>