fill_cent.f90

来自「基于linux操作下的fortran快速付立变换的程序」· F90 代码 · 共 31 行

        subroutine fill_cent(md1,md3,lot,n1dfft,max3,m3,n3,zf,zw)
        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
        zw(1,i1,i3)=zf(1,i1,i3)
        zw(2,i1,i3)=zf(2,i1,i3)
        enddo
90      continue

!       Fill the center region with zeros
        do 100,i3=max3+2,n3-m3+max3+1
        do i1=1,n1dfft
        zw(1,i1,i3)=0.d0
        zw(2,i1,i3)=0.d0
        enddo
100     continue

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

        return
        end