From 89a77a9f327edad799e65e7aa5d296d77818f907 Mon Sep 17 00:00:00 2001 From: ignis Date: Mon, 28 Apr 2014 05:56:51 +0900 Subject: [PATCH] domain size given as parameter is used --- x_fftw.f90 | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/x_fftw.f90 b/x_fftw.f90 index a431471..7fc30b3 100644 --- a/x_fftw.f90 +++ b/x_fftw.f90 @@ -267,7 +267,7 @@ CONTAINS !------------------------------------------------------------------------------! ! in Fourier space it is (nx / 2 + 1) complex numbers along kx-axis do ix = 1, nx + 1, 2 - akx(ix) = real((ix - 1) / 2, 8) + akx(ix) = real((ix - 1) / 2, 8) * (two / lx) !ksj modification1 akx(ix + 1) = akx(ix) coskx2(ix) = dcos(half * akx(ix)) sinkx2(ix) = dsin(half * akx(ix)) @@ -278,21 +278,21 @@ CONTAINS end do ! in Fourier space ky-axis is distributed among the processors do iy = 1, nz - aky(iy) = real(myid * nz + iy - 1, 8) - if (aky(iy) > (0.5D0 * real(ny, 8))) aky(iy) = aky(iy) - real(ny, 8) + aky(iy) = real(myid * nz + iy - 1, 8) * (two / ly) !ksj modification1 + if (aky(iy) > (0.5D0 * real(ny, 8) * (two / ly))) aky(iy) = aky(iy) - real(ny, 8) * (two / ly) cosky2(iy) = dcos(half * aky(iy)) sinky2(iy) = dsin(half * aky(iy)) rezkay(iy) = 0 - if (dabs(aky(iy)) > (real(ny, 8)) / 3.0D0) rezkay(iy) = 1 + if (dabs(aky(iy)) > (real(ny, 8) * (two / ly)) / 3.0D0) rezkay(iy) = 1 end do ! in Fourier space the z wavenumbers are aligned along the second index do iz = 1, ny - akz(iz) = real(iz - 1, 8) - if (akz(iz) > (0.5D0 * real(nz_all, 8))) akz(iz) = akz(iz) - real(nz_all, 8) + akz(iz) = real(iz - 1, 8) * (two / lz) !ksj modification1 + if (akz(iz) > (0.5D0 * real(nz_all, 8) * (two / lz))) akz(iz) = akz(iz) - real(nz_all, 8) * (two / lz) coskz2(iz) = dcos(half * akz(iz)) sinkz2(iz) = dsin(half * akz(iz)) rezkaz(iz) = 0 - if (dabs(akz(iz)) > (real(nz_all, 8)) / 3.0D0) rezkaz(iz) = 1 + if (dabs(akz(iz)) > (real(nz_all, 8) * (two / lz)) / 3.0D0) rezkaz(iz) = 1 end do ! Definition of the array ialias.