unfill_cent.f90

基于linux操作下的fortran快速付立变换的程序

        subroutine unfill_cent(md1,md3,lot,n1dfft,max3,m3,n3,zw,zf)
        implicit real*8 (a-h,o-z)
        dimension zw(2,lot,n3),zf(2,md1,md3)

!       Here, zero and positive frequencies
        do 90,i3=1,max3+1
        do i1=1,n1dfft
        zf(1,i1,i3)=zw(1,i1,i3)
        zf(2,i1,i3)=zw(2,i1,i3)
        enddo
90      continue

!       Here, negative frequencies
        do 110,i3=max3+2,m3
        do i1=1,n1dfft
        zf(1,i1,i3)=zw(1,i1,i3+n3-m3)
        zf(2,i1,i3)=zw(2,i1,i3+n3-m3)
        enddo
110     continue

        return
        end