<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<meta name="Generator" content="Microsoft Exchange Server">
<!-- converted from text --><style><!-- .EmailQuote { margin-left: 1pt; padding-left: 4pt; border-left: #800000 2px solid; } --></style>
</head>
<body>
<meta content="text/html; charset=UTF-8">
<style type="text/css" style="">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
<p>Hi Martin,</p>
<p><br>
</p>
<p>So when I set M_rbc=1 and tauRelaxT to be O(deltaT), the simulation quickly generates NaN. Is this forcing considered too strong<font size="3"> or I probably did something wrong?</font></p>
<p>I tried other values of tauRelaxT, if the values are much larger, say deltaT=23 and tauRelaxT=230, I have waves generated but with a much smaller amplitude, which is probably expected since the forcing is pretty weak.</p>
<p>From what I have seen so far, it seems like if I want to generate waves at xx=500pt, I should probably start forcing the waves at xx=400pt and slowly ramp up the forcing?</p>
<p><br>
</p>
<p>I would appreciate some insights in this.</p>
<p><br>
</p>
<p>Thanks</p>
<p><br>
</p>
<p>yangxin</p>
<p><br>
</p>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> MITgcm-support <mitgcm-support-bounces@mitgcm.org> on behalf of Martin Losch <Martin.Losch@awi.de><br>
<b>Sent:</b> Wednesday, April 4, 2018 5:24:27 AM<br>
<b>To:</b> MITgcm Support<br>
<b>Subject:</b> Re: [MITgcm-support] setting up with RBCS package</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">Yangxin,<br>
<br>
OBCS and RBCS are fundamentally different. OBCS is an explicit package for open boundary conditions. It has some restoring options for a sponge layer. The R in RBCS is for “restoring” or “relaxation". RBCS can be used to do 3D restoring in the entire domain,
or according to a mask that you have to specify. I think S/R rbcs_add_tendency is the file to look at to understand how it works. For your problem you may want to use a restoring time scale O(deltaT) and a mask that is nonzero only at xx=500. More answers
below,<br>
<br>
Martin<br>
<br>
> On 3. Apr 2018, at 20:21, Yangxin He <y67he@uwaterloo.ca> wrote:<br>
> <br>
> Hello there,<br>
> <br>
> Here is the model I want to set up:<br>
> Assume I have nx=3000 pt, ny=50 pt, nz=50 pt. From xx=1 to 409 pt, I want to set up a layer of stretch grid to damp out any reflections coming from the right. From xx=2501 to 3000 pt, I will also set up a layer of stretch grid. I have rigid lid at the top
and no slip at the bottom. The domain where all the physics comes from is in xx=500 to 2500 pt. These can be done by initializing the grid at the start and it is not a problem.<br>
> Now I want to force a mode-one internal wave at xx=500 pt so that the waves propagate to the right into the domain. Since the forcing is not at the boundary, OBCS does not seem to be the right tool so I took a look at the RBCS package, which got me quite
confused.<br>
> Here are my questions:<br>
> 1) I checked the exp4 in the MITgcm verification folder and also a sample code from Jody. It seems that as long as RBCS is used, OBCS is also used. Do they have to go together or is this an coincidence? Why?<br>
They can be used independently.<br>
> 2) In my head, I thought setting my routine using RBCS would be similar to use OBCS prescribe except telling the model which location I want the force to turn on. But after reading the MITgcm documentation, the RBCS package seems much more complicated. From
my understanding, <br>
> 2.1) M_rbc / \tau_T acts like a ramp up function? Any specific rules to pick these numbers or just something slow and smooth will do?<br>
Not sure what you are refering to. I only see a ramp down flag (rbcsVanishingTime). This probably applies to specific cases, but is not generally required.<br>
> 2.2) T/U/V_rbc is the T/U/V I want to force at the location? So T_rbc is a function of (x,y,z,t), which means the relaxTFile is a matrix of (:,:,:,:)?<br>
Yes<br>
> 3) If anyone has any simple codes/files set up similar to my case and are willing to share with me, I would really appreciate it!<br>
> <br>
> Thanks<br>
> <br>
> Yangxin<br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> MITgcm-support@mitgcm.org<br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
<br>
_______________________________________________<br>
MITgcm-support mailing list<br>
MITgcm-support@mitgcm.org<br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</div>
</span></font>
</body>
</html>