function CubeFrame, nL, color=Color, corners=Corners
;+
; NAME:
;	CubeFrame
; PURPOSE:
;	Draws outline of 3D data cube
; CATEGORY:
;	Graphics
; CALLING SEQUENCE:
;	Cube = CubeFrame(nDim, color=Color)
; INPUTS:
;	nDim		scalar		cube size in data units
; OPTIONAL INPUT PARAMETERS:
;	color=Color	scalar		color index used for drawing
;							(if omitted then !p.color is used)
; OUTPUTS:
;	(None)
; CALLS:
; COMMON BLOCKS:
; SIDE EFFECTS:
; RESTRICTIONS:
; PROCEDURE:
;	Make sure the proper 3D transformation is set up.
;	The cube drawn extends from 0 to nDim-1 in all three dimensions.
; MODIFICATION HISTORY:
;	Feb-1997, Paul Hick (UCSD)
;-
on_error, 2

if keyword_set(Corners) then begin
	Cube = intarr(3,16)

	Cube(0,*) = [0,0,0,0,0,1,1,0,1,0,0,1,1,1,1,1]
	Cube(1,*) = [0,1,1,0,0,0,1,1,0,0,1,1,0,0,1,1]
	Cube(2,*) = [0,0,1,1,0,0,0,0,1,1,1,1,1,0,0,1]

	Cube = nL*Cube

	if n_elements(Color) ne 0 then begin
		plots, Cube(*,0: 7), /t3d, color=Color
		plots, Cube(*,8:15), /t3d, color=Color
	endif

endif else begin

	s = (1.+findgen(nL))/(nL)	; 1./nL to 1 along cube side
	t1 = replicate(1.,nL)
	t0 = replicate(0.,nL)

	Cube = fltarr(3,1+7*nL,2)

	j = 0
	k = 0

	Cube(0,j,k) = 0.
	Cube(1,j,k) = 0.
	Cube(2,j,k) = 0.

	i = j+1  &  j = j+nL
	Cube(0,i:j,k) = t0
	Cube(1,i:j,k) = s
	Cube(2,i:j,k) = t0

	i = j+1  &  j = j+nL
	Cube(0,i:j,k) = t0
	Cube(1,i:j,k) = t1
	Cube(2,i:j,k) = s

	i = j+1  &  j = j+nL
	Cube(0,i:j,k) = t0
	Cube(1,i:j,k) = 1-s
	Cube(2,i:j,k) = t1

	i = j+1  &  j = j+nL
	Cube(0,i:j,k) = t0
	Cube(1,i:j,k) = t0
	Cube(2,i:j,k) = 1-s

	i = j+1  &  j = j+nL
	Cube(0,i:j,k) = s
	Cube(1,i:j,k) = t0
	Cube(2,i:j,k) = t0

	i = j+1  &  j = j+nL
	Cube(0,i:j,k) = t1
	Cube(1,i:j,k) = s
	Cube(2,i:j,k) = t0

	i = j+1  &  j = j+nL
	Cube(0,i:j,k) = 1-s
	Cube(1,i:j,k) = t1
	Cube(2,i:j,k) = t0

	j = 0
	k = 1

	Cube(0,j,k) = 1.
	Cube(1,j,k) = 0.
	Cube(2,j,k) = 1.

	i = j+1  &  j = j+nL
	Cube(0,i:j,k) = 1-s
	Cube(1,i:j,k) = t0
	Cube(2,i:j,k) = t1

	i = j+1  &  j = j+nL
	Cube(0,i:j,k) = t0
	Cube(1,i:j,k) = s
	Cube(2,i:j,k) = t1

	i = j+1  &  j = j+nL
	Cube(0,i:j,k) = s
	Cube(1,i:j,k) = t1
	Cube(2,i:j,k) = t1

	i = j+1  &  j = j+nL
	Cube(0,i:j,k) = t1
	Cube(1,i:j,k) = 1-s
	Cube(2,i:j,k) = t1

	i = j+1  &  j = j+nL
	Cube(0,i:j,k) = t1
	Cube(1,i:j,k) = t0
	Cube(2,i:j,k) = 1-s

	i = j+1  &  j = j+nL
	Cube(0,i:j,k) = t1
	Cube(1,i:j,k) = s
	Cube(2,i:j,k) = t0

	i = j+1  &  j = j+nL
	Cube(0,i:j,k) = t1
	Cube(1,i:j,k) = t1
	Cube(2,i:j,k) = s

	Cube = nL*Cube
endelse

return, Cube  &  end