From 4a8a67d241bd2e139f3063caa9dd4209c98036c8 Mon Sep 17 00:00:00 2001
From: ignis <ignis@postech.ac.kr>
Date: Sun, 24 Dec 2017 20:19:22 +0900
Subject: [PATCH] fixed bug in ialias array construction

---
 x_fftw.f90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/x_fftw.f90 b/x_fftw.f90
index 7fc30b3..4efc0b2 100644
--- a/x_fftw.f90
+++ b/x_fftw.f90
@@ -298,7 +298,7 @@ CONTAINS
     ! Definition of the array ialias.
     ! The array ialias is just the number of wavenumbers at (i,j,k) that have
     ! their magnitude higher than nx/3.  This is needed in dealiasing procedures.
-    rnx3 = real(nx/3, 8)
+    rnx3 = real(2*nx, 8) / real(3*lx, 8)
     do k = 1,nz
        if (abs(aky(k)) .gt. rnx3) ialias(:,:,k) = 1
        do j = 1,ny