
Hello, I want to compute the magnetic flux over a disk which length is the minimum possible i.e. the resolution. I wanted to approximate this quantity using : sp=ds.disk(center,normal,(1,'pc'),(1*res,'pc')) flux=sp.calculate_flux("sqrt_magx","sqrt_magy","sqrt_magz","sqrt_magx") Where I defined the quantities sqrt_magx (root square of the magnetic field along the x direction) with units of Gauss^0.5, so the output of calculate_flux has Gauss*cm^2 units. However, it is not possible to use calculate_flux on YTDisks and only calculate_isocontour_flux is available for this data container. Do you know how I can obtain exactly the same result for the magnetic flux using the available calculate_isocontour_flux instead of the calculate_flux function? Cheers,

Hi Guido, Hmm, I think this is a limitation of how we have implemented the flux calculations, which are defined across triangles (and thus with built-in normal vectors). My suggestion would be to define a field that represents the edge of the area you want to use, and then use that as the surface you extract. You could do this with something that was along the disk_height field, for instance, and I believe if you extract from that data object the triangles will all be within it. On Tue, Sep 10, 2019 at 6:13 PM Guido Granda Muñoz <guidogranda@gmail.com> wrote:
Hello, I want to compute the magnetic flux over a disk which length is the minimum possible i.e. the resolution. I wanted to approximate this quantity using :
sp=ds.disk(center,normal,(1,'pc'),(1*res,'pc')) flux=sp.calculate_flux("sqrt_magx","sqrt_magy","sqrt_magz","sqrt_magx")
Where I defined the quantities sqrt_magx (root square of the magnetic field along the x direction) with units of Gauss^0.5, so the output of calculate_flux has Gauss*cm^2 units. However, it is not possible to use calculate_flux on YTDisks and only calculate_isocontour_flux is available for this data container.
Do you know how I can obtain exactly the same result for the magnetic flux using the available calculate_isocontour_flux instead of the calculate_flux function?
Cheers, _______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org
participants (2)
-
Guido Granda Muñoz
-
Matthew Turk