For a Spherical cap:
https://en.wikipedia.org/wiki/Spherical_cap

Let d = (r-h)
d = sqrt(r*r - a*a)
r = (a*a + h*h) / 2*h
theta = arcsin(a/r)

- Brian
Hasty trig math, check for errors:-)