pro oc_search_hd_test

; Program makes calls to oc_search_hd* routines to test them.

; Author: HBT
;
; History:
;   10-Oct-2005: Original version HBT

  ra_test	= 0.5d
  dec_test	= 0.5d
  radius_test	= 0.1d

  result=oc_search_hd_single(ra_test, dec_test, radius_test, nstars)

  print, 'Found ' + st(nstars) + ' stars at RA = ' + st(ra_test) + $
    ', dec = ' + st(dec_test) + ', radius = ' + st(radius_test) + ' radians'

  print, 'Found: ' + st(nstars) + ' stars, RA_sum = ' + st(total(result[*].ra)) + ', Dec_sum = ' + st(total(result[*].dec))

  plot, result[*].ra, result[*].dec, psym=3

  print

  print, 'Results by HBT 10-Oct-2005: '
  print, 'Input: 0.5d, 0.5d, 0.1d radians'
  print, 'Output: 234 stars, RA_sum = 117.03723, Dec_sum = 116.43132'

  print
  print, 'RA_sum, Dec_sum are the sums of the RA and Dec of all individual stars, and can be used as a checksum'

end