Cleanups.

src/omp/poisson_fft.f90

@@ -30,7 +30,6 @@ function init(xdirps, ydirps, zdirps) result(poisson_fft)
       class(dirps_t), intent(in) :: xdirps, ydirps, zdirps
 
       type(omp_poisson_fft_t) :: poisson_fft
-      integer :: nx, ny, nz
 
       call poisson_fft%base_init(xdirps, ydirps, zdirps, SZ)
 

src/poisson_fft.f90

@@ -55,8 +55,6 @@ subroutine base_init(self, xdirps, ydirps, zdirps, sz)
       class(dirps_t), intent(in) :: xdirps, ydirps, zdirps
       integer, intent(in) :: sz
 
-      integer :: nx, ny, nz
-
       self%nx = xdirps%n; self%ny = ydirps%n; self%nz = zdirps%n
 
       allocate (self%ax(self%nx), self%bx(self%nx))
@@ -84,7 +82,7 @@ subroutine waves_set(self, xdirps, ydirps, zdirps, sz)
       real(dp) :: w, wp, rlexs, rleys, rlezs, xtt, ytt, ztt, xt1, yt1, zt1
       complex(dp) :: xt2, yt2, zt2, xyzk
 
-      integer :: i, j, k, ka, kb, ix, iy, iz
+      integer :: i, j, ka, kb, ix, iy, iz
 
 
       nx = xdirps%n; ny = ydirps%n; nz = zdirps%n
@@ -156,12 +154,7 @@ subroutine waves_set(self, xdirps, ydirps, zdirps, sz)
          zk2(i) = cmplx(1._dp, 1._dp, kind=dp)*(nz*wp/zdirps%L)**2
       end do
 
-      !do k = !sp%xst(3), sp%xen(3)
-         !do j = sp%xst(2), sp%xen(2)
-      ! kxyz array is in z pencil orientation, with nz=nz/2+1 !!nooo
-      ! kxyz array is in x pencil orientation.
-      !do k = 1, nx*ny/SZ
-      print*, 'set kxyz array'
+      print*, 'set waves array'
       ! TODO: do loop ranges below are valid only for single rank runs
       do ka = 1, nz/2 + 1
          do kb = 1, ny/sz