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
builtin 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
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, _______________________________________________ ytusers mailing list  ytusers@python.org To unsubscribe send an email to ytusersleave@python.org
participants (2)

Guido Granda Muñoz

Matthew Turk