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

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

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
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20101230/47ede59a/attachment.html>