<div><div dir="auto">The obcs package should do this too.</div></div><div dir="auto"><br></div><div dir="auto">Matt</div><div dir="auto"><br></div><div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Sep 16, 2019 at 5:45 AM Edward W Doddridge <<a href="mailto:ewd@mit.edu">ewd@mit.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">



<div style="word-wrap:break-word;line-break:after-white-space">
Dear Kevin,
<div><br>
</div>
<div>Is there as reason that you can’t use the RBCS package? <a href="https://mitgcm.readthedocs.io/en/latest/phys_pkgs/rbcs.html" target="_blank">https://mitgcm.readthedocs.io/en/latest/phys_pkgs/rbcs.html</a></div>
<div><br>
</div>
<div>Cheers,</div>
<div>Ed</div>
<div><br>
<div>
<div style="color:rgb(0,0,0);font-family:Palatino;font-size:16px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">
<span class="m_-5386335672115366828Apple-style-span" style="font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;border-collapse:separate;border-spacing:0px"><br class="m_-5386335672115366828Apple-interchange-newline">
<hr style="font-family:Tahoma;font-size:13px">
<span class="m_-5386335672115366828Apple-style-span" style="font-family:Tahoma;font-size:13px"><font face="Tahoma" size="2"><span dir="ltr" style="font-size:10pt"><font size="1"><span style="font-size:13px"><font size="2"><b>Edward
 Doddridge</b><br>
<i>Postdoctoral researcher</i></font></span></font></span></font></span></span></div>
<div style="color:rgb(0,0,0);font-family:Palatino;font-size:16px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">
<span class="m_-5386335672115366828Apple-style-span" style="font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;border-collapse:separate;border-spacing:0px"><span class="m_-5386335672115366828Apple-style-span" style="font-family:Tahoma;font-size:13px"><font face="Tahoma" size="2"><span dir="ltr" style="font-size:10pt"><font size="1"><span style="font-size:13px"><font size="2"><i>Earth,
 Atmospheric and Planetary Sciences</i></font></span></font></span></font></span></span></div>
<div style="color:rgb(0,0,0);font-family:Palatino;font-size:16px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">
<span class="m_-5386335672115366828Apple-style-span" style="font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;border-collapse:separate;border-spacing:0px"><span class="m_-5386335672115366828Apple-style-span" style="font-family:Tahoma;font-size:13px"><font face="Tahoma" size="2"><span dir="ltr" style="font-size:10pt"><font size="1"><span style="font-size:13px"><font size="2"><i>Massachusetts
 Institute of Technology</i></font></span></font></span></font></span></span></div>
<div style="color:rgb(0,0,0);font-family:Palatino;font-size:16px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">
<div style="font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal">
<span class="m_-5386335672115366828Apple-style-span" style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-align:-webkit-auto;border-spacing:0px"><span class="m_-5386335672115366828Apple-style-span" style="font-family:Tahoma;font-size:13px"><font face="Tahoma" size="2"><span dir="ltr" style="font-size:10pt"><font size="1"><span style="font-size:13px"><font size="2"><i><br>
</i></font></span></font></span></font></span></span></div>
<div style="font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal">
<span class="m_-5386335672115366828Apple-style-span" style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-align:-webkit-auto;border-spacing:0px"><span class="m_-5386335672115366828Apple-style-span" style="font-family:Tahoma;font-size:13px"><font face="Tahoma" size="2"><a href="http://www.doddridge.me" target="_blank">www.doddridge.me</a></font></span></span></div>
</div>
</div>
<div><br>
<blockquote type="cite">
<div>On 16 Sep 2019, at 04:48, Kevin Ha <<a href="mailto:kevin.ha.pro@gmail.com" target="_blank">kevin.ha.pro@gmail.com</a>> wrote:</div>
<br class="m_-5386335672115366828Apple-interchange-newline">
<div>
<div dir="ltr">
<div>Dear community,</div>
<div><br>
</div>
<div>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 ?</div>
<div><br>
</div>
<div>Thank you in advance</div>
<div><br>
</div>
<div>Kind regards</div>
<div><br>
</div>
<div>Kevin.</div>
<div><br>
</div>
<div><br>
</div>
</div>
<span id="m_-5386335672115366828cid:f_k0m5zkyd8"><packages.conf></span><span id="m_-5386335672115366828cid:f_k0m5zck72"><gendata.m></span><span id="m_-5386335672115366828cid:f_k0m5zcjv1"><eedata></span><span id="m_-5386335672115366828cid:f_k0m5zcj80"><data></span><span id="m_-5386335672115366828cid:f_k0m5zkx64"><external_forcing.F></span><span id="m_-5386335672115366828cid:f_k0m5zkxi5"><FFIELDS.h></span><span id="m_-5386335672115366828cid:f_k0m5zkwu3"><CPP_OPTIONS.h></span><span id="m_-5386335672115366828cid:f_k0m5zky27"><MOM_FLUXFORM_OPTIONS.h></span><span id="m_-5386335672115366828cid:f_k0m5zkxt6"><ini_fields.F></span><span id="m_-5386335672115366828cid:f_k0m5zkyk9"><SIZE.h></span>_______________________________________________</div></blockquote></div></div></div><div style="word-wrap:break-word;line-break:after-white-space"><div><div><blockquote type="cite"><div><br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</div>
</blockquote>
</div>
<br>
</div>
</div>

_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></div></div>