[Nek5000-users] Visit reader for double precision f files

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

Hi Philipp, Hank,

I _think_ there may still be an issue for the Pn-Pn-2 case if we are
dumping out for the full-restart capability, in which case pressure
does not occupy (N+1)^3 memory sites, but only (N-1)^3, with the remainder
being zero.

This would not influence analysis of velocity fields, however.

Paul

PS - it's also possible that this deficit of pressure data exists
      only for the restart files, rather than for standard files that
      are written in double precision, but I'd have to check on that.


On Fri, 14 Sep 2012, nek5000-users at lists.mcs.anl.gov wrote:

> Dear Hank,
> thanks a lot for your answer. Indeed, now that I checked myself it turns out 
> that double precision files work perfectly in version 2.4.0 and onwards of 
> visit, so this was false alarm, sorry! I will inquire why we thought in the 
> first place that 64 bit would not work.
>
> Thanks again,
> Philipp
>
>
>
>
> On Fri, 14 Sep 2012 00:46:11 +0200, <nek5000-users at lists.mcs.anl.gov> wrote:
>
>> Hi Philip,
>> 
>> With regards to VisIt, I am happy to investigate.  IIRC, this issue came up
>> before and we decided that the reader did work.  (But possibly it slipped
>> off my queue and it doesn't work!!)
>> 
>> If you can make a file available to me (hchilds at lbl.gov), I can pursue.
>> 
>> With respect to the file format issue, someone else more knowledgable will
>> need to answer.
>> 
>> Best,
>> Hank
>> 
>> On Thu, Sep 13, 2012 at 3:41 PM, <nek5000-users at lists.mcs.anl.gov> wrote:
>> 
>>> Dear all,
>>> We were wondering whether it would be possible to include support for 64
>>> bit
>>> .f files in the Visit reader for Nek5000. So far, only 32 bit files are
>>> supported, which we are hardly using any longer except for those files to
>>> visualise with Visit; however we have quite a number of f files 
>>> (snapshots)
>>> from various simulations that we use for different purposes such as the
>>> computation of spectra, POD modes etc. for which we naturally chose 64 bit
>>> due to accuracy consideration. To be able to directly visualise such files
>>> in Visit without converting them to 32 bit would be very helpful indeed.
>>> 
>>> Another aspect that came up recently was the question on how to indicate 
>>> in
>>> a f file whether or not the pressure is mapped to the Pn mesh if the
>>> Pn/Pn-2
>>> method is used. Normally, restart files would not perform this mapping,
>>> snapshot files however would do. We were wondering whether a simple
>>> solution
>>> could be to choose different letters in the header (say "P" for pressure 
>>> on
>>> the Pn mesh and "Q" for pressure on the Pn-2 mesh)?
>>> 
>>> We could give a try for both things here if you think it would not be too
>>> difficult...
>>> 
>>> Best regards,
>>> Philipp
>>> 
>>> 
>>> _______________________________________________
>>> Nek5000-users mailing list
>>> Nek5000-users at lists.mcs.anl.gov
>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users