;  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 phase_jk_wL, H, d, xh, j, k, wL

  psirefrhatL = refr_jk_wL(H, d, xh, j, k, wL)
  ihm = complex(0,1) * H * wL
  psihatL = (((H*wL)^2.+1.)^(-0.25))*((2+ihm)/(2-ihm))^(0.5)*psirefrhatL
  return, psihatL

end

pro phase_jk_wL_test

  H = 32.
  d = 1.
  j=1
  k=1
  psihatL = phase_jk_wL(H, d, 0., 0, 0, wL)
  psijkhatL = phase_jk_wL(H, d, 0., j, k, wL)
  window, 0, xs=1000,ys=500
  !p.multi=[0,2,0]
  plot, wL, abs(psihatL), xr=[-20,20]/(4.*H)
  oplot, wL, abs(psijkhatL), /li, ps=-4

end