[yt-users] calculate the Mdot in yt.

Britton Smith brittonsmith at gmail.com
Thu May 18 17:33:50 PDT 2017


Hi Nathan,

Yes!  I'll do that as soon as I can.

Britton

On Thu, May 18, 2017 at 5:16 PM, Nathan Goldbaum <nathan12343 at gmail.com>
wrote:

> 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
>>
>
> _______________________________________________
> 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/20170518/d247240e/attachment.html>


More information about the yt-users mailing list