<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div>Hi Melinda and all,</div><div><br></div>Just for clarification, Jill's code works in yt 2.x. The 3.0 beta version (to be released soon) of this code would be:<div><br></div><div>import yt</div><div>pf = yt.load(fname)</div><div>sphere = pf.sphere([0.,0.,0.], (pf.domain_right_edge[0]/3.1e18, "pc"))<br>rad_profile = yt.create_profile(sphere, "radius", ["density","temperature","velocity_x"],</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>extrema={"radius":(0.0,pf.domain_right_edge[0]/3.1e18*yt.units*pc)},</div><div> n_bins=ninterp, units={"radius":"pc"}, logs={"radius":False})</div><div><br></div><div>radius = rad_profile.x</div><div><br></div><div>Then you can access the fields like so:</div><div><br>rad_profile["density"]<br>rad_profile["temperature"]<br>rad_profile["velocity_x"]</div><div><br></div><div>Also, yt 3.0 is aware of spherical geometries, so you could also just use the unbinned quantities:</div><div><br></div><div>sp = pf.sphere([0.0,0.0,0.0], (pf.domain_right_edge[0]/3.1e18, "pc"))</div><div><div><div>radius = sp["radius"].in_units("pc")</div><div>density = sp["density"]</div><div><br></div><div>etc...</div><div><br></div><div>Best,</div><div><br></div><div>John ZuHone</div><div><br></div><div>On Jul 9, 2014, at 5:25 PM, Jill Naiman <<a href="mailto:jnaiman@gmail.com">jnaiman@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div dir="ltr"><div>The following has worked for me:<br><br>pf = load(fname)<br>sphere = pf.h.sphere([0.,0.,0.], (pf.domain_right_edge[0]/3.1e18, "pc"))<br>rad_profile = BinnedProfile1D(sphere, ninterp, "Radiuspc", 0.0, \<br>
pf.domain_right_edge[0]/3.1e18, log_space=False)<br>radius = rad_profile['Radiuspc']<br> <br>rad_profile.add_fields("Density")<br>rad_profile.add_fields("Temperature")<br>
rad_profile.add_fields("x-velocity")<br><br><br></div>And then I just plot using "fig"<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jul 9, 2014 at 2:47 PM, Melinda Soares-Furtado <span dir="ltr"><<a href="mailto:msoares.physics@gmail.com" target="_blank">msoares.physics@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">Is there a simple way to create density/temp/velocity profiles using 1D FLASH data? I've done this in the past for 3-dimensional data, which is very user friendly.<div>
<br></div><div>Is the best method of attack is <span style="font-family:arial,sans-serif;font-size:13px">dumping the data to np arrays and then plotting from these arrays?</span><br>
<div><br clear="all"><div><div dir="ltr"><div><font face="georgia, serif">Regards,</font></div><font face="georgia, serif">Melinda Soares-Furtado</font><div><font face="georgia, serif"><a href="mailto:msoares@princeton.edu" target="_blank">msoares@princeton.edu</a><br>
</font><div></div></div><div><font face="georgia, serif"><a href="http://cargo.ucsc.edu/msoares" target="_blank">http://cargo.ucsc.edu/msoares</a></font></div><div><font face="georgia, serif">415-860-0438</font><br>
<br></div></div></div>
</div></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>
_______________________________________________<br>yt-users mailing list<br><a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org<br></blockquote></div><br></div></body></html>