<div dir="ltr"><div class="gmail_extra">Thank you Josh. This is very useful information. I will try it out on my problem.</div><div class="gmail_extra"><br></div><div class="gmail_extra">So this means visit does not make use of the high order definition of the solution but simply does linear interpolation ???</div>
<div class="gmail_extra"><br></div><div class="gmail_extra">Regards</div><div class="gmail_extra">praveen</div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Nov 8, 2013 at 11:59 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 dir="ltr">Praveen, <div><br></div><div>Another thing you could do is to output your solution onto a uniformly spaced grid within each element(as opposed to the standard GLL points).  I have found this to be particularly useful in cases where relatively large, high order elements are used, like I was seeing for my Euler solver.</div>

<div><br></div><div>To do this, in usrchk, include the following: (I'm doing this out of memory, so you may want to double check this)</div><div><br></div><div>subroutine userchk()</div><div><br></div><div>  include 'SIZE'</div>

<div>  include 'TOTAL'</div><div>  include 'RESTART'    ! Not 100% sure about this one</div><div><br></div><div>  integer mynxo</div><div><br></div><div>  mynxo = 10</div><div>  ifreguo = .true.    ! Tells prepost to interpolate to uniformly spaced grid in each element</div>

<div>  nxo = mynxo      ! nxo is the number of points (in each direction) to use within the elements</div><div><br></div><div>end</div><div><br></div><div>Theoretically, nxo can be anything from 2 (which would map your solution to linear finite elements) to lxo (which is set in your SIZE file).  Although, as I'm writing this, there may be a hard-limit on nxo set somewhere in the source code, so maybe the developers could verify this for me?</div>

<div><br></div><div>Using this will output your field files to  "_reg" files.  Keep in mind that these files will be for output purposes only--restarting from these would incur some potentially heavy interpolation errors.</div>

<div><br></div><div>Hope this helps!</div><div><br></div><div>Josh</div></div></blockquote></div><br><br></div></div>