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...<div><br></div><div>Best,</div><div>Hank<br><br><div class="gmail_quote">
On Thu, Dec 30, 2010 at 11:03 AM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div dir="ltr">Hi Hank,<div><br></div><div>It works w/o any issues. Thanks.<br><div><br clear="all">Regards<br>Shriram Jagannathan<div><div></div><div class="h5"><br>
<br><br><div class="gmail_quote">On 30 December 2010 12:44, Shriram Jagannathan <span dir="ltr"><<a href="mailto:shriramjegan@gmail.com" target="_blank">shriramjegan@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Hi Hank,<div><br></div><div>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 !</div><div><br clear="all">
Regards<br><font color="#888888">Shriram Jagannathan</font><div><div></div><div><br>
<br><br><div class="gmail_quote">On 30 December 2010 12:30, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Shririam,<div><br></div><div>I would do:</div><div>OpenDatabase("file.nek5000") # use real filename</div><div>AddPlot("Pseudocolor", "pressure") # or whatever variable you want</div><div>
DrawPlots()</div>
<div>Lineout((1,1,1),(2,2,2)) # takes lineout from point (1,1,1) to point (2,2,2)</div><div>SetActiveWindow(2)</div><div>sw = SaveWindowAttributes()</div><div>sw.format = sw.CURVE</div><div>sw.family = 0</div><div>sw.fileName = "myname_line1" </div>
<div>SetSaveWindowAttributes(sw) # saves out myname_line1.curve</div><div>SaveWindow()</div><div>DeleteAllPlots()</div><div>SetActiveWindow(1)</div><div>Lineout((-1, 1, -1), (2, 4, 5))</div><div>SetActiveWindow(2)</div>
<div>
sw.fileName = "myname_line2"</div><div>SetSaveWindowAttributes(sw)</div><div>SaveWindow()</div><div>DeleteAllPlots()</div><div>SetActiveWindow(1)</div><div>....</div><div><br></div><div>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.</div>
<div><br></div><div>Best,</div><div>Hank<div><div></div><div><br><br><div class="gmail_quote">On Thu, Dec 30, 2010 at 10:22 AM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Hank,<div><br></div><div>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. </div>
<div><br clear="all">Regards<br>Shriram Jagannathan<div><div></div><div><br>
<br><br><div class="gmail_quote">On 30 December 2010 12:15, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>Hi Shiriam,<div><br></div><div>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.</div>
<div><br></div><div>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...</div><div><br>
</div>
<div>Best,</div><div>Hank</div><div><br><div class="gmail_quote"><div><div></div><div>On Thu, Dec 30, 2010 at 10:12 AM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div></div><div><div dir="ltr">Hi,<div><br></div><div>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..</div>
<div><br></div><div>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 ?</div>
<div><br clear="all">Regards<br>Shriram<br>
</div></div>
<br></div></div>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div>
<br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div></div></div></div>
<br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div></div></div>
<br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div></div></div></div>
</blockquote></div><br></div></div></div></div></div>
<br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div>