;  Given a number of points and a spacing, return an
;  array of the wavelet fourier series and the
;  array of frequencies (2 pi/wavelength) for the
;  j, k wavelet as defined by Sato.

function sato_jk_wL, n, d, j, k, wL

  i = complex(0,1)

  psihatL = sato_j0_wL(n, d, j, wL)  
  psihatL = psihatL * exp(- i * k * 2.^(-float(j)) * wL )
  return, psihatL

end

pro sato_jk_wL_test

  n = 256
  d = 20./n
  j=1
  k=1
  psihatL = sato_mother_wL(n, d, wL)
  psijkhatL = sato_jk_wL(n, d, j,k, wL)
  window, 0, xs=1000,ys=500
  !p.multi=[0,2,0]
  plot, wL, abs(psihatL), xr=[-20,20], yr=[0,1.2], /ys
  oplot, wL, abs(psijkhatL), /li, ps=-4

end