<div dir="ltr"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Feb 13, 2014 at 4:10 PM,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id=":ls" style="overflow:hidden">Hi Praveen<br>
<br>
(Thu, Feb 13, 2014 at 01:02:38PM +0530) <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> :<br>
<div class="">>                         This is very inefficient of course, and I could<br>
> only run it in serial.<br>
<br>
</div>It turns out finally that I missed that in your first message, and now it<br>
works perfectly, thanks a lot!<br>
<br>
Best regards,<br>
<span class="HOEnZb"><font color="#888888">--<br>
Ismael</font></span></div></blockquote></div><br>Ok. But it is good to still vary delta and see how the vorticity changes. You can compute error in recomputed vorticity and try to find a good value of delta.</div><div class="gmail_extra">
<br></div><div class="gmail_extra">Also due to roundoff errors, it may be good to enforce continuity of computed velocity. Looking at the eddy_psi_omega example, I think this is done by</div><div class="gmail_extra"><br></div>
<div class="gmail_extra">call dsavg(vx)</div><div class="gmail_extra">call dsavg(vy)</div><div class="gmail_extra"><br></div><div class="gmail_extra">Best</div><div class="gmail_extra">praveen</div></div>