FIND_LIMB also finds the location of the center of the optical sun and it's
radius. A sample calling sequence is:
IDL> find_limb, img, x0, y0, r0
where img contains the input image and x0,y0,r0
receive the output center
location and diameter of the disk. FIND_LIMB only works on whole images. Also see FIT_LIMB.