[yt-users] calculate the Mdot in yt.

[Nathan Goldbaum]

Hi Britton,

Any chance you can file an issue on github along with a runnable test
script?

Thanks!

Nathan

On Thu, May 18, 2017 at 7:10 PM Britton Smith <brittonsmith at gmail.com>
wrote:

> Hi Kazem,
>
> 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.
>
> 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.
>
> Does anyone have any ideas?
>
> Britton
>
> On Tue, May 16, 2017 at 11:15 PM, Kazem Ardaneh <kazem.arrdaneh at gmail.com>
> wrote:
>
>> Dear All,
>>
>> I have a question about calculating the Mdot in yt.
>> I tried to set some concentric spheres shells with deferent radius and
>> calculate the mdot through them.==>> Mdot vs R
>>
>> I used an object in yt as:
>>
>> surf = ds.surface(sph,"radius",r)
>>
>> where
>>
>> sph = ds.sphere(cmax, (2, "kpc"))
>> and r is radius of each shell.
>>
>> After that I used:
>> flux = surf.calculate_flux("x-velocity","y-velocity","z-velocity",
>> "density")
>>
>> which calculate the flux density, Mdot.
>>
>> However it gives the mdot zero which is wrong.
>>
>> 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.
>>
>> Does anyone has idea about this issue? how to fix it?
>>
>>
>> Best Regards,
>>
>> kazem
>>
>>
>>
>>
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20170519/596e0c8f/attachment.html>