Fortran结果中的FFTW只包含零

我一直在尝试编写一个简单的程序,使用fftw3在一维输入阵列上执行fft。这里我用地震记录作为输入。然而,输出数组只包含零。

program fft

    use functions
    implicit none
    include 'fftw3.f90'
    integer nl,row,col
    double precision, allocatable :: data(:,:),time(:),amplitude(:)
    double complex, allocatable :: out(:)
    integer*8 plan


    open(1,file='test-seismogram.xy')
    nl=nlines(1,'test-seismogram.xy')
    allocate(data(nl,2))
    allocate(time(nl))
    allocate(amplitude(nl))
    allocate(out(nl/2+1))
    do row = 1,nl
        read(1,*,end=101) data(row,1),data(row,2)
        amplitude(row)=data(row,2)
    end do
    101 close(1)


    call dfftw_plan_dft_r2c_1d(plan,nl,amplitude,out,FFTW_R2HC,FFTW_PATIENT)
    call dfftw_execute_dft_r2c(plan, amplitude, out)
    call dfftw_destroy_plan(plan)


    do row=1,(nl/2+1)
        print *,out(row)
    end do


    deallocate(data)
    deallocate(amplitude)
    deallocate(time)
    deallocate(out)
end program fft

这个 nlines() 函数是一个用来计算文件行数的函数,它可以正常工作。它在名为函数的模块中定义。

http://www.fftw.org/fftw3_doc/Fortran-Examples.html

我可能只是犯了一个非常简单的逻辑错误,但我实在搞不清楚这里到底出了什么问题。任何提示都会很有帮助。

这几乎就是整个输出的样子比如:-

           .
           .
           .
           (0.0000000000000000,0.0000000000000000)
           (0.0000000000000000,0.0000000000000000)
           (0.0000000000000000,0.0000000000000000)
           (0.0000000000000000,0.0000000000000000)
           (0.0000000000000000,0.0000000000000000)
           .
           .
           .

我的疑问是关于fftw的,因为在SO上有fftw的标签,所以我希望这个问题不是离题