c    na2 = na/2, 2**lp = na
       double complex ffttbl(namh*lpm,2), wa(nam), wb(nam)
       double complex wc(0:nam,0:nam)
       integer fftlst(namh*lpm)
       common / fftwrk / ffttbl, fftlst