[yt-users] calculate the Mdot in yt.

[Kazem Ardaneh]

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




-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20170517/9d317ed4/attachment.htm>