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

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

Cool!  Keep the questions coming and I'll help out with whatever I can.
 It's very little effort for me to put together little scripts like that...

Best,
Hank

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

> Hi Hank,
>
> It works w/o any issues. Thanks.
>
> Regards
> Shriram Jagannathan
>
>
>
> On 30 December 2010 12:44, Shriram Jagannathan <shriramjegan at gmail.com>wrote:
>
>> 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
>>>
>>>
>>
>
> _______________________________________________
> 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/002cc91b/attachment.html>