fmt  = "%9.5f"   ; use to 'pretty print' values

    x    = (/ 1002, 1017, 1018, 1020, 1018, 1027, \ 
              1028, 1030, 1012, 1012,  982, 1012, \
              1001,  996,  995, 1011, 1027, 1025, \ 
              1030, 1016,  996, 1006, 1002,  982  /)
    N    = dimsizes(x)

    cf   = ezfftf (x)
    ptint("")
    print("++++++:   EZFFTF")
    print("")
    printVarSummary(cf)

    cfn  = cf/(N/2)       ; normalize

    p    = cf(0,1:)^2  + cf(1,1:)^2     ; do not square the mean
    pn   = cfn(0,1:)^2 + cfn(1,1:)^2
                         
    print("")
    print("======: real/imag")
    print( sprintf(fmt, cf(0,:)) +" "+sprintf(fmt, cf(1,:))+" norm=>"+ \
           sprintf(fmt,cfn(0,:)) +" "+sprintf(fmt,cfn(1,:)) )       
    print("")
    print("======: periodogram: ezfftf")
    print( sprintf(fmt, p) +" "+sprintf(fmt, pn) )

    print("")
    print("++++++:   CFFTF")
    print("")
    CF   = cfftf(x, 0.0, 0)
    printVarSummary(CF)
    
    CFn  = CF/N            ; normalize
    print("")
    print("======")
    print( sprintf(fmt, CF(0,:)) +" "+sprintf(fmt, CF(1,:))+" "+ \
           sprintf(fmt,CFn(0,:)) +" "+sprintf(fmt,CFn(1,:))      )

    P    = CF(0,1:)^2  + CF(1,1:)^2     ; do not square the mean
    Pn   = CFn(0,1:)^2 + CFn(1,1:)^2
    print("")
    print("======: periodogram: cfftf")
    print( sprintf(fmt, P) +" "+sprintf(fmt, Pn) )