[Nek5000-users] Post-Processing from avg_all files using Nek / Visit

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

Hi Shririam,

I would do:
OpenDatabase("file.nek5000") # use real filename
AddPlot("Pseudocolor", "pressure") # or whatever variable you want
DrawPlots()
Lineout((1,1,1),(2,2,2)) # takes lineout from point (1,1,1) to point (2,2,2)
SetActiveWindow(2)
sw = SaveWindowAttributes()
sw.format = sw.CURVE
sw.family = 0
sw.fileName = "myname_line1"
SetSaveWindowAttributes(sw)  # saves out myname_line1.curve
SaveWindow()
DeleteAllPlots()
SetActiveWindow(1)
Lineout((-1, 1, -1), (2, 4, 5))
SetActiveWindow(2)
sw.fileName = "myname_line2"
SetSaveWindowAttributes(sw)
SaveWindow()
DeleteAllPlots()
SetActiveWindow(1)
....

In VisIt 2.2, this will be much easier ... you won't have to make a
Pseudocolor plot and do a Lineout query ... you will be able to add the
Lineout directly via a curve plot.

Best,
Hank

On Thu, Dec 30, 2010 at 10:22 AM, <nek5000-users at lists.mcs.anl.gov> wrote:

> Hi Hank,
>
> Thanks. I think I can automate (with python) till I get Visit to do a line
> plot. Do you know a way I can just export the data from line plot to a text
> file ? I might have about 10 line plots (900 time steps) which I would like
> to automate.
>
> Regards
> Shriram Jagannathan
>
>
>
> On 30 December 2010 12:15, <nek5000-users at lists.mcs.anl.gov> wrote:
>
>>
>> Hi Shiriam,
>>
>> It is pretty easy to automate this in VisIt.  As we discussed at the
>> user's meeting, VisIt treats each spectral element as a 5x5x5, 7x7x7, etc
>> group of hexahedrons.  Within each hexahedron, there is linear
>> interpolation.
>>
>> So: if you want help automating this in VisIt, that should be easy.  As
>> far as accuracy, Nek will clearly do the best job, but hopefully you now
>> understand better what VisIt is doing...
>>
>>  Best,
>> Hank
>>
>> On Thu, Dec 30, 2010 at 10:12 AM, <nek5000-users at lists.mcs.anl.gov>wrote:
>>
>>> Hi,
>>>
>>> I would like to post-process quantities from avg_all files (rms, rm2,
>>> avg). Particularly I am looking at line plots (mainly) or contours on a
>>> plane. I had a question regarding the approach that I might be taking..
>>>
>>> If I load the files, say rmsblah.f0000* into nek using the load
>>> subroutine and fill the ux,uy,uz arrays with that of <u^2>, <v^2>, <w^2> and
>>> call hpts(), how would the interpolation work ? Would it interpolate <u^2>
>>> with the interpolation that is used for velocities ?. The other option would
>>> be to do a line plot from Visit and export it to a text file. Also, I am not
>>> sure how Visit interpolates. Which is more likely to be accurate ?
>>>
>>> Regards
>>> Shriram
>>>
>>> _______________________________________________
>>> 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
>>
>>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20101230/d3491634/attachment.html>