From 5650694bf692b06462124d2425c24b9f760e83c1 Mon Sep 17 00:00:00 2001 From: Antoine Hoffmann <antoine.hoffmann@epfl.ch> Date: Tue, 26 Sep 2023 16:17:29 +0200 Subject: [PATCH] cleaning --- src/fourier_mod.F90 | 9 +-------- 1 file changed, 1 insertion(+), 8 deletions(-) diff --git a/src/fourier_mod.F90 b/src/fourier_mod.F90 index 50b44d5a..091525d5 100644 --- a/src/fourier_mod.F90 +++ b/src/fourier_mod.F90 @@ -344,13 +344,6 @@ END SUBROUTINE fft1D_plans CALL apply_ExB_NL_factor(ikyG,ExB_NL_factor) CALL apply_ExB_NL_factor(ikyF,ExB_NL_factor) CALL apply_ExB_NL_factor(ikxG,ExB_NL_factor) - ! ! TEST apply inverse to test the identity - ! invfactor = 1._xp/ExB_NL_factor - ! CALL apply_ExB_NL_factor(ikxF,invfactor) - ! CALL apply_ExB_NL_factor(ikyG,invfactor) - ! CALL apply_ExB_NL_factor(ikyF,invfactor) - ! CALL apply_ExB_NL_factor(ikxG,invfactor) - ! ! END TEST ENDIF ! Anti Aliasing DO iky = 1,local_nky @@ -508,7 +501,7 @@ END SUBROUTINE fft1D_plans ! Treat the result with the ExB NL factor DO iky = 1,NY_/2+1 DO ix = 1,local_nx_ - tmp_kyx(iky,ix) = tmp_kyx(iky,ix)!*inv_ExB_NL_factor(iky,ix) + tmp_kyx(iky,ix) = tmp_kyx(iky,ix)*inv_ExB_NL_factor(iky,ix) ENDDO ENDDO ! Back to Fourier space in the third buffer -- GitLab