<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><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Nov 8, 2013 at 1:53 AM, <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">The latest visit version 2.7 seems to be unable to read my results, I get<div><br></div><div>
<p style="margin:0px">VisIt could not read from the file "/private/tmp/re10000/channel.nek5000".</p>
<p style="margin:0px"><br></p>
<p style="margin:0px">The generated error message was:</p>
<p style="margin:0px"><br></p>
<p style="margin:0px">There was an error opening /private/tmp/re10000/channel.nek5000. It may be an invalid file.</p>
<p style="margin:0px">VisIt tried using the following file format readers to open the file: Nek5000, Silo</p>
<p style="margin:0px"><br></p>
<p style="margin:0px">If you know the specific format reader VisIt should use to</p>
<p style="margin:0px">read this data, you can use Open As... (GUI) or</p>
<p style="margin:0px">'-o <file>,<plugin> (CL arg.) and identify that specific reader</p>
<p style="margin:0px">for VisIt to try. This will possibly give more information on</p>
<p style="margin:0px">the exact error.</p><p style="margin:0px"><br></p><p style="margin:0px">Has anybody had success with 2.7 ?</p><p style="margin:0px"><br></p><p style="margin:0px">My file looks like this</p><p style="margin:0px">
<br></p><p style="margin:0px">filetemplate: channel.fld%02d</p><p style="margin:0px">firsttimestep: 1</p><p style="margin:0px">numtimesteps: 161</p><p style="margin:0px"><br></p><p style="margin:0px">Thanks</p><p style="margin:0px">
praveen</p></div></div>
<br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>Josh Camp<br><br>"All that is necessary for the triumph of evil is that good men do nothing" -- Edmund Burke
</div>