PRO dd_load_ct, SHOW = show
; Load a nice color table to go from red to blue in a continous way

red = intarr(256)
green = intarr(256)
blue = intarr(256)

path = fltarr(3,256)

; ~Red to Yellow (3-96)
start = [1.0,0.3,0.]
stop = [1.,1.,0.]
steps = indgen(94)
; steps goes from 0-->1
steps = steps/FLOAT(n_elements(steps)-1)
steps = (steps)^0.75
for is=0,n_elements(steps)-1 do path(*,3+is) = (1.-steps(is))*start + $
	steps(is)*stop

; Yellow to Green (97-150)
start = [1.0,1.,0.]
stop = [0.5,1.,0.]
steps = indgen(54)
; steps goes from 0-->1
steps = steps/FLOAT(n_elements(steps)-1)
for is=0,n_elements(steps)-1 do path(*,97+is) = (1.-steps(is))*start + $
	steps(is)*stop

; Green to Cyan (151-234)
start = [0.0,1.0,0.5]
stop = [0.,1.,0.9]
steps = indgen(84)
; steps goes from 0-->1
steps = steps/FLOAT(n_elements(steps)-1)
for is=0,n_elements(steps)-1 do path(*,151+is) = (1.-steps(is))*start + $
	steps(is)*stop

; Cyan to White (235-254)
start = [0.3,1.0,0.9]
stop = [1.,1.,1.]
steps = indgen(20)
; steps goes from 0-->1
steps = steps/FLOAT(n_elements(steps)-1)
steps = steps^0.75
for is=0,n_elements(steps)-1 do path(*,235+is) = (1.-steps(is))*start + $
	steps(is)*stop

; add some others
; 0=black, 1 = white, 255 = white
path(*,0) = [0.,0.,0.]
path(*,1) = [1.,1.,1.]
path(*,2) = [0.,0.,0.5]
path(*,255) = [1.,1.,1.]

; change from path to RGB
red = FIX(path(0,*)*255)
green = FIX(path(1,*)*255)
blue = FIX(path(2,*)*255)

tvlct, red, green , blue

if KEYWORD_SET(SHOW) then begin
  y = indgen(256)
  plot, y, /NODATA, BACK=2, COLOR=1
  for ip=3,254 do oplot, [ip],[y(ip)],COLOR=ip,PSYM=4
end

RETURN
END