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@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@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
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@lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________
Nek5000-users mailing list
Nek5000-users@lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________
Nek5000-users mailing list
Nek5000-users@lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________
Nek5000-users mailing list
Nek5000-users@lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users