[Nek5000-users] visit visualization question

[nek5000-users at lists.mcs.anl.gov]

Praveen,

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.

To do this, in usrchk, include the following: (I'm doing this out of
memory, so you may want to double check this)

subroutine userchk()

  include 'SIZE'
  include 'TOTAL'
  include 'RESTART'    ! Not 100% sure about this one

  integer mynxo

  mynxo = 10
  ifreguo = .true.    ! Tells prepost to interpolate to uniformly spaced
grid in each element
  nxo = mynxo      ! nxo is the number of points (in each direction) to use
within the elements

end

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?

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.

Hope this helps!

Josh


On Fri, Nov 8, 2013 at 1:53 AM, <nek5000-users at lists.mcs.anl.gov> wrote:

> The latest visit version 2.7 seems to be unable to read my results, I get
>
> VisIt could not read from the file "/private/tmp/re10000/channel.nek5000".
>
>
> The generated error message was:
>
>
> There was an error opening /private/tmp/re10000/channel.nek5000. It may be
> an invalid file.
>
> VisIt tried using the following file format readers to open the file:
> Nek5000, Silo
>
>
> If you know the specific format reader VisIt should use to
>
> read this data, you can use Open As... (GUI) or
>
> '-o <file>,<plugin> (CL arg.) and identify that specific reader
>
> for VisIt to try. This will possibly give more information on
>
> the exact error.
>
>
> Has anybody had success with 2.7 ?
>
>
> My file looks like this
>
>
> filetemplate:         channel.fld%02d
>
> firsttimestep: 1
>
> numtimesteps: 161
>
>
> Thanks
>
> praveen
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>


-- 
Josh Camp

"All that is necessary for the triumph of evil is that good men do nothing"
-- Edmund Burke
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20131108/528de7e6/attachment.html>