switch_cent.f90

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

        subroutine switch_cent(n1dfft,max2,m2,n2,lot,n1,lzt,zt,zw)
        implicit real*8 (a-h,o-z)
        dimension zw(2,lot,n2),zt(2,lzt,n1)

!       Here, zero and positive frequencies
        do 90,j=1,n1dfft
        do 90,i=1,max2+1
        zw(1,j,i)=zt(1,i,j)
        zw(2,j,i)=zt(2,i,j)
90      continue

!       Fill the center region with zeros
        do 100,i=max2+2,n2-m2+max2+1
        do 100,j=1,n1dfft
        zw(1,j,i)=0.d0
        zw(2,j,i)=0.d0
100     continue

!       Here, negative frequencies
        if(m2>=max2+2)then
         do 110,j=1,n1dfft
         do 110,i=max2+2,m2
         zw(1,j,i+n2-m2)=zt(1,i,j)
         zw(2,j,i+n2-m2)=zt(2,i,j)
110      continue
        endif

        return
        end