;+
; NAME:
;       INDICES_ARRAY
; PURPOSE:
;       Does the inverse of ARRAY_INDICES
; CALLING SEQUENCE:
;       INDICES_ARRAY, x, y, array
; EXAMPLES:
;       array[indices_array(x, y, array)]
;           Gets the values at the location of the pair values
;           (x[i],y[i]) of array.
; KEYWORDS:
;       NONE
; MODIFICATION HISTORY:
;       Before 05/06/2007 Written by Eduard Westra
;-
FUNCTION indices_array, x, y, array
  IF N_ELEMENTS(x) GT 1 THEN xl = REFORM(x) ELSE xl = x
  IF N_ELEMENTS(y) GT 1 THEN yl = REFORM(y) ELSE yl = y
  IF SIZE(array, /N_DIMENSIONS) NE 2 THEN MESSAGE, 'Currently only capable of 2d arrays'
  sz = SIZE(array)
  RETURN, y##REPLICATE(sz[1], N_ELEMENTS(xl)) + x#REPLICATE(1, N_ELEMENTS(yl))
END