Hi Guys,
Is there a way to quantify values in yt by solid angle around a point?
For example I have a simulation of star clusters forming within a turbulent GMC. When a bubble is present I want to calculate the covering fraction of the shell. One way to do this would be using where the radiation pressure is absorbed and calculate the fraction of solid angle (around the oldest cluster) radiation pressure is absorbed above a minimum value out to a given radius. See attached for projections of this simulation at a time with a bubble shell! (bottom left shows radiation pressure absorption).
Has anyone written any code/script to calculate solid angle before?
Many thanks,
Libby
Hi Libby,
There are a couple ways to do this that come to mind right away. Unfortunately, calculating the exact solid angle that a cell subtends can be a bit tricky, but otherwise this is precisely the kind of analysis yt is designed to do.
1) Bin based on the solid angle, approximating this by, for instance, assuming that the cell is hit face on such that the area that it covers is dx^2. (This assumes cubical cells, though.) 2) Actually integrate out and accumulate the values you are interested in. This would require using the HEALpix camera, which is described here http://yt.enzotools.org/doc/visualizing/volume_rendering.html and on the blog. You could, for instance, take just a ProjectionTransferFunction and apply that to an integration of a field that has cutoffs based on the solid angle. The HEALpix camera performs an allsky integration calculated from a central point; this can be used for quantitative integration and analysis of field values.
There may be other ways, though  maybe somebody else has some other ideas?
Matt
On Mon, May 16, 2011 at 12:33 PM, Elizabeth HarperClark hclark@cita.utoronto.ca wrote:
Hi Guys,
Is there a way to quantify values in yt by solid angle around a point?
For example I have a simulation of star clusters forming within a turbulent GMC. When a bubble is present I want to calculate the covering fraction of the shell. One way to do this would be using where the radiation pressure is absorbed and calculate the fraction of solid angle (around the oldest cluster) radiation pressure is absorbed above a minimum value out to a given radius. See attached for projections of this simulation at a time with a bubble shell! (bottom left shows radiation pressure absorption).
Has anyone written any code/script to calculate solid angle before?
Many thanks,
Libby
 Elizabeth HarperClark MA MSci PhD Candidate, Astrophysics, UofT
www.astro.utoronto.ca/~hclark hclark@cita.utoronto.ca AIM: edphc1 MSN: edphc1@hotmail.com Skype: eharperclark Office phone: 4169785759
ytusers mailing list ytusers@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/ytusersspacepope.org
participants (2)

Elizabeth HarperClark

Matthew Turk