<div dir="ltr">Thanks Martin. Yes, it makes sense when you put it that way. <div><br></div><div>Dhruv</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Mar 11, 2018 at 2:32 PM, Martin Losch <span dir="ltr"><<a href="mailto:Martin.Losch@awi.de" target="_blank">Martin.Losch@awi.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Dhruv,<br>
<br>
if you ignore other sources for your tracer, you’ll have from discretizing d Tracer/ dt  = 1/tauRelax * ( Tracer_r - Tracer ) in time (it’s more complicated, but just conceptually):<br>
Tracer(n+1) = Tracer(n) + dt/tauRelax * ( Tracer_r - Tracer(n) )<br>
for tauRelax = deltaT, this means Tracer(n+1) = Tracer_r, ie you replace your tracer with the restoring value. A shorter restoring time scale doesn’t make too much sense to me. And if you use the above logic for tauRelax << deltaT then (deltaT/tauRelax>>1) and I am not quite sure what you would be getting. Maybe your advection scheme limits the effect of the restoring after you get funny values from deltaT/tauRelax >> 1<br>
<br>
I suggest that you don’t use tauRelax<deltaT.<br>
<br>
Martin<br>
<div><div class="h5"><br>
> On 10. Mar 2018, at 21:27, Dhruv Balwada <<a href="mailto:db194@nyu.edu">db194@nyu.edu</a>> wrote:<br>
><br>
> Hello,<br>
><br>
> I am trying to do some surface tracer restoration experiments with relaxation time scales that are very fast (1s), extreme limit. This number is smaller than the time step in my experiment. I use the RBCS package for this. I am noticing that if the restoration time scale is smaller than the time step, then no tracer enters the domain. However, things work fine/as expected if the restoration time scale is larger than or equal to the time step. Is there some reason why this is the case? The restoration time scale seems to enter as a constant, not sure why it is tied to the numerics.<br>
><br>
> Thanks,<br>
> Dhruv<br>
</div></div>> ______________________________<wbr>_________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a><br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
<br>
______________________________<wbr>_________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
</blockquote></div><br></div>