<div dir="ltr">This is great. Now I just need to figure out how to get the total velocity normal to the surface. </div><div class="gmail_extra"><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><div><div><font face="georgia, serif">415-860-0438</font><br><br></div></div></div></div></div>
<br><br><div class="gmail_quote">On Thu, Jul 3, 2014 at 7:00 PM, Nathan Goldbaum <span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@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 Melinda,<div><br></div><div>I believe you'll be able to construct that by extracting the x,y,z velocity fields and Density fields from the surface object:</div><div><br></div><div>vx, vy, vz = surf['x-velocity'], surf['y-velocity'], surf['z-velocity']</div>


<div>dens = surf['Density']</div><div><br></div><div>These will be four flat numpy arrays which correspond to the values of the velocity and density fields at the surface boundary.</div><div><br></div><div>Hope that helps,</div>


<div><br></div><div>Nathan</div></div><div class="gmail_extra"><br><br><div class="gmail_quote"><div><div class="h5">On Thu, Jul 3, 2014 at 2:19 PM, Melinda Soares-Furtado <span dir="ltr"><<a href="mailto:msoares.physics@gmail.com" target="_blank">msoares.physics@gmail.com</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 class="h5"><div dir="ltr"><span style="font-size:13.333333969116211px;font-family:arial,sans-serif">Hello, folks. I have a code that is computing the flux across a sphere of a given radius, as shown below</span><span style="font-size:13.333333969116211px;font-family:arial,sans-serif">. This code provides me with the summed flux across the surface, however what I am hoping to do is to bin this flux up into various velocities ranging from v_min to v_max. Thus, determine the mass flux across the surface for a chosen velocity bin.</span><div style="font-size:13.333333969116211px;font-family:arial,sans-serif">



<br></div><div style="font-size:13.333333969116211px;font-family:arial,sans-serif">Does anyone know of a simple way to read out flux for a given velocity (normal to the surface) rather than the summed total? </div><div style="font-size:13.333333969116211px;font-family:arial,sans-serif">



<br></div><div style="font-size:13.333333969116211px;font-family:arial,sans-serif">Code:</div><div style="font-size:13.333333969116211px;font-family:arial,sans-serif"><pre style="margin-top:0px;margin-bottom:0px;padding:5px 0px;font-family:'Bitstream Vera Sans Mono',monospace;font-size:13px;color:rgb(0,0,0)">

from yt.mods import *
pf=load("/trove/msoares/50Stars_LRefine5_Derefine0.05_velocity/50Stars_hdf5_plt_cnt_0200")
rad=0.1/pf['pc']
sp = pf.h.sphere([0,0,0],2.0*rad)
surf = pf.h.surface(sp,"Radius",rad)
flux = surf.calculate_flux("x-velocity","y-velocity","z-velocity","Density")
print flux</pre><pre style="margin-top:0px;margin-bottom:0px;padding:5px 0px;font-family:'Bitstream Vera Sans Mono',monospace;font-size:13px;color:rgb(0,0,0)"><br></pre></div><div><div dir="ltr"><div><font face="georgia, serif">Regards,</font></div>



<font face="georgia, serif">Melinda Soares-Furtado</font></div></div>
</div>
<br></div></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>
<br></blockquote></div><br></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>