<div dir="ltr">Hi Nathan,<div><br></div><div>Yes!  I'll do that as soon as I can.</div><div><br></div><div>Britton</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, May 18, 2017 at 5:16 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>Hi Britton,</div><div><br></div><div>Any chance you can file an issue on github along with a runnable test script?</div><div><br></div><div>Thanks!</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Nathan</div></font></span><div class="HOEnZb"><div class="h5"><div><br><div class="gmail_quote"><div>On Thu, May 18, 2017 at 7:10 PM Britton Smith <<a href="mailto:brittonsmith@gmail.com" target="_blank">brittonsmith@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>Hi Kazem,<div><br></div><div>I have looked into this and I think there may be an bug here, but I am not 100% positive.  This may be obvious, but the reason the flux is coming out zero is because the surface object is empty.  Just doing surf[<some_field>] returns an empty array.</div><div><br></div><div>I have played around with this a bit using the IsolatedGalaxy dataset and have been unable to create a non-empty surface using the "radius" field.  My first suspicion was that this was a units issue as the marching_cubes function doesn't check that the contour value and the field values are in the same units.  This is true and probably should be fixed in its own right.  However, making sure my surface value was in the same units as what the radius field returns didn't help either.</div><div><br></div><div>Does anyone have any ideas?</div><div><br></div><div>Britton</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, May 16, 2017 at 11:15 PM, Kazem Ardaneh <span><<a href="mailto:kazem.arrdaneh@gmail.com" target="_blank">kazem.arrdaneh@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">Dear All, <div><br></div><div><div>I have a question about calculating the Mdot in yt.</div><div>I tried to set some concentric spheres shells with deferent radius and calculate the mdot through them.==>> Mdot vs R</div><div><br></div><div>I used an object in yt as:</div><div><br></div><div><div style="margin:0px;line-height:normal;font-family:menlo"><span style="font-variant-ligatures:no-common-ligatures">surf = ds.surface(sph,</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">"radius"</span><span style="font-variant-ligatures:no-common-ligatures">,r)</span></div></div><div><br></div><div>where</div><div><br></div><div><div style="margin:0px;line-height:normal;font-family:menlo"><span style="font-variant-ligatures:no-common-ligatures">sph = ds.sphere(cmax, (</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">2</span><span style="font-variant-ligatures:no-common-ligatures">, </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">"kpc"</span><span style="font-variant-ligatures:no-common-ligatures">))</span></div><div style="margin:0px;line-height:normal;font-family:menlo"><span style="font-variant-ligatures:no-common-ligatures">and r is radius of each shell.</span></div><div style="margin:0px;line-height:normal;font-family:menlo"><span style="font-variant-ligatures:no-common-ligatures"><br></span></div><div style="margin:0px;line-height:normal;font-family:menlo">After that I used:</div><div style="margin:0px;line-height:normal;font-family:menlo;color:rgb(195,55,32)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">flux = surf.calculate_flux(</span><span style="font-variant-ligatures:no-common-ligatures">"x-<wbr>velocity"</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">,</span><span style="font-variant-ligatures:no-common-ligatures">"y-velocity"</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">,</span><span style="font-variant-ligatures:no-common-ligatures">"z-<wbr>velocity"</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">,</span><span style="font-variant-ligatures:no-common-ligatures">"density"</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">)</span></div><div style="margin:0px;line-height:normal;font-family:menlo;color:rgb(195,55,32)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)"><br></span></div><div style="margin:0px;line-height:normal;font-family:menlo;color:rgb(195,55,32)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">which calculate the flux density, Mdot.</span></div><div style="margin:0px;line-height:normal;font-family:menlo;color:rgb(195,55,32)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)"><br></span></div><div style="margin:0px;line-height:normal;font-family:menlo;color:rgb(195,55,32)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">However it gives the mdot zero which is wrong. </span></div><div style="margin:0px;line-height:normal;font-family:menlo;color:rgb(195,55,32)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)"><br></span></div><div style="margin:0px;line-height:normal;font-family:menlo;color:rgb(195,55,32)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">If I set the surfaces as the density iso surface, then I got the mdot as function of density, which it work. However, as a function of radius it fives zero. </span></div><div style="margin:0px;line-height:normal;font-family:menlo;color:rgb(195,55,32)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)"><br></span></div><div style="margin:0px;line-height:normal;font-family:menlo"><font color="#000000">Does anyone has idea about this issue? how to fix it?</font></div></div><div><font color="#000000" face="menlo"><br></font></div><div><font color="#000000" face="menlo"><br></font></div><div>
<div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div>Best Regards,</div><div><br></div><div>kazem</div><div><br></div></div><br class="m_3680139490764529271m_-1715021517934428691m_5286831609696955821Apple-interchange-newline"><br class="m_3680139490764529271m_-1715021517934428691m_5286831609696955821Apple-interchange-newline">

</div>
<br></div></div><br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
<br></blockquote></div><br></div>
______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
</blockquote></div></div>
</div></div><br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
<br></blockquote></div><br></div>