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,


Elizabeth Harper-Clark MA MSci
PhD Candidate, Astrophysics, UofT

AIM: edphc1
MSN: edphc1@hotmail.com
Skype: eharperclark
Office phone: 416-978-5759