research/dns-hit3d-fdm
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

fixed bug in ialias array construction

Browse source
This commit is contained in:
ignis 2017-12-24 20:19:22 +09:00
parent db7732025a
commit 4a8a67d241
1 changed files with 1 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
x_fftw.f90
View file

@ -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