unswitch_cent.f90

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

        subroutine unswitch_cent(n1dfft,max2,m2,n2,lot,n1,lzt,zw,zt)
        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
        zt(1,i,j)=zw(1,j,i)
        zt(2,i,j)=zw(2,j,i)
90      continue

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

        return
        end