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

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

Hi Hank,

Thank you very much for the script. I will try this now. We have VisIt 2.1
on local machines but the cluster is powered with VisIt 1.12 only. Thanks
again !

Regards
Shriram Jagannathan


On 30 December 2010 12:30, <nek5000-users at lists.mcs.anl.gov> wrote:

> 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
>>
>>
>
> _______________________________________________
> 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/fe770fc7/attachment.html>