<div dir="ltr">My bad, yeah, this is yt-2.x stuff. Also, ninterp is just some number to bin things in. In my code I have the note:<br><br>ninterp = 1400 # for example: nblockx =1500 => ninterp ~ 500-1400<br><br></div>
<div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jul 9, 2014 at 3:58 PM, John ZuHone <span dir="ltr"><<a href="mailto:jzuhone@gmail.com" target="_blank">jzuhone@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 style="word-wrap:break-word"><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 style="white-space:pre-wrap"> </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><div class="h5"><div><br></div><div>On Jul 9, 2014, at 5:25 PM, Jill Naiman <<a href="mailto:jnaiman@gmail.com" target="_blank">jnaiman@gmail.com</a>> wrote:</div><br><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" target="_blank">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" target="_blank">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>
</blockquote></div></div></div><br></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>