<div dir="ltr">Hi Morgan,<div><br></div><div>Sorry for the delay in replying. I've spent some time looking at the polar coordinate handling in yt-3.0 this morning and I believe I have come up with the current status.</div>
<div><br></div><div>As it stands, the actual creation of a nice plot like yt normally does -- through the PlotWindow object -- is not yet possible. There are a few more steps that need to happen before everything can be nicely integrated like that. However, I have refreshed my memory on the coordinate handling system, and I have been able to make some nice images relatively easily.</div>
<div><br></div><div>The code had to have a few changes to support this; I've put these changes into a pull request here:</div><div><br></div><div><a href="https://bitbucket.org/yt_analysis/yt-3.0/pull-request/82/fix-flash-polar-cylindrical-pixelization/diff">https://bitbucket.org/yt_analysis/yt-3.0/pull-request/82/fix-flash-polar-cylindrical-pixelization/diff</a><br>
</div><div><br></div><div>You can get them by doing:</div><div><br></div><div>cd $YT_DEST/src/yt-hg/</div><div><br></div><div>(or wherever your yt-3.0 installation is) and doing:</div><div><br></div><div>hg pull -r 12101b782ac8 <a href="https://bitbucket.org/MatthewTurk/yt-3.0">https://bitbucket.org/MatthewTurk/yt-3.0</a></div>
<div>hg up -r 12101b782ac8</div><div>python2.7 setup.py build_ext -i -f</div><div><br></div><div>(the -f will force every extension to rebuild; if you get funny errors with selector objects being the wrong size, you may need to do: "find . -name '*.pyx' | xargs touch" and then another build_ext run.)</div>
<div><br></div><div>This script then created the correct output for me: <a href="http://paste.yt-project.org/show/3788/">http://paste.yt-project.org/show/3788/</a></div><div><br></div><div>I've tested it on a few NIF items. The changes I made were specifically in how the data is expected to be ordered; I'd incorrectly had cylindrical as rtz and polar as rzt in FLASH. I've switched these in the PR.</div>
<div><br></div><div>Roman, Suoqing and Anthony, I've listed you as reviewers on the PR because I know you have all looked at polar and/or cylindrical data in the past. If you could give some feedback, and especially let me know if I have done something wrong.</div>
<div><br></div><div>In the future, I would like to have things working such that we can do SlicePlot( ... ) like normal. That will require a somewhat detailed discussion with PlotWindow experts about how to set up bounds for a window, but I think it should be doable in the not-too-distant future.</div>
<div><br></div><div>Let me know if that works for you, and if you need any more help!</div><div><br></div><div>Best,</div><div><br></div><div>Matt</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Aug 15, 2013 at 8:58 PM, Morgan MacLeod <span dir="ltr"><<a href="mailto:morganmacleod@gmail.com" target="_blank">morganmacleod@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">Dear yt gurus,<div><br></div><div>I've enjoyed using yt for a number of FLASH problems in cartesian geometry. Right now I'm just starting on a setup in FLASH 2D polar (eventually also in 3D cylindrical) geometries. </div>
<div><br></div><div>These are ordered </div><div>2D polar: r, theta</div><div>3D cylindrical: r, theta, z</div><div><br></div><div>I was intrigued when I noticed this image on the yt-blog (<a href="http://blog.yt-project.org/post/WhatsUpWith30.html" target="_blank">http://blog.yt-project.org/post/WhatsUpWith30.html</a>) -- which is almost exactly what I would like to do. </div>
<div><br></div><div><img alt="http://blog.yt-project.org/attachments/cylindrical_pixelizer.png" style="max-width:100%;min-height:200px;vertical-align:middle;border:0px;color:rgb(153,153,153);font-family:'Droid Sans',sans-serif;font-size:14px;line-height:20px;background-color:rgb(18,20,23);width:200px"><br>
</div><div><br></div><div><br></div><div>I have updated to yt-3.0, and tried the simplest "Sliceplot" of the data:</div><div><span style><br></span></div><div><font face="arial, helvetica, sans-serif" style><span>fn = "const_rho_g_atm_hdf5_plt_cnt_0013"</span><br>
</font></div><div><font color="#000000" face="arial, helvetica, sans-serif" style> pf = load(fn)</font></div><div><span style><span style="font-family:arial,helvetica,sans-serif"> </span><span style="font-size:13px;font-family:arial,helvetica,sans-serif">fig = SlicePlot(pf, 2,'dens').save("myfig.png")</span></span></div>
<div><br></div><div>This reads loads the 2D polar plot file correctly, but renders the result in a cartesian sense. <br></div><div><br></div><div>Does anyone have advice/recipes for how to proceed with Slices rendered "cylindrically"? I understand that this feature set is not yet fully developed but I was curious if there was a solution (or parts of a solution) floating around out there. </div>
<div><br></div><div>Thanks so much,</div><div><br></div><div>Morgan MacLeod</div><div>-------------------------------------------------</div><div>Graduate Student</div><div>Astronomy & Astrophysics</div><div>UC Santa Cruz</div>
</div>
<br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>