From ff4a338a98c5f2148b4e2ba94cd922d2469218f8 Mon Sep 17 00:00:00 2001
From: Antoine Cyril David Hoffmann <ahoffman@spcpc606.epfl.ch>
Date: Mon, 5 Oct 2020 14:52:21 +0200
Subject: [PATCH] Forgot Antialiasing...
---
src/compute_Sapj.F90 | 15 ++++++++++++---
1 file changed, 12 insertions(+), 3 deletions(-)
diff --git a/src/compute_Sapj.F90 b/src/compute_Sapj.F90
index a7a40a8c..e4f422f9 100644
--- a/src/compute_Sapj.F90
+++ b/src/compute_Sapj.F90
@@ -78,7 +78,12 @@ SUBROUTINE compute_Sapj
ENDIF
ENDDO nloope
- !CALL forward_FFT(Sepj(ip,ij,:,:)) !**then put the sum back to fourier space (spares n FFT)
+ !! Antialiasing
+ DO ikr = ikrs,ikre
+ DO ikz = ikzs, ikze
+ Sepj(ip,ij,ikr,ikz) = Sepj(ip,ij,ikr,ikz) * AA_r(ikr) * AA_z(ikz)
+ ENDDO
+ ENDDO
ENDDO jloope
ENDDO ploope
@@ -146,10 +151,14 @@ SUBROUTINE compute_Sapj
ENDDO nloopi
- !CALL forward_FFT(Sipj(ip,ij,:,:)) ! **then put the sum back to fourier space (spares n FFT)
+ !! Antialiasing
+ DO ikr = ikrs,ikre
+ DO ikz = ikzs, ikze
+ Sipj(ip,ij,ikr,ikz) = Sipj(ip,ij,ikr,ikz) * AA_r(ikr) * AA_z(ikz)
+ ENDDO
+ ENDDO
ENDDO jloopi
ENDDO ploopi
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-
END SUBROUTINE compute_Sapj
--
GitLab