From c3d043a4cddeee72baeebb70af0f0ed40e6a810b Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Thu, 28 Sep 2023 10:15:03 +0200
Subject: [PATCH] cleaning
---
src/fourier_mod.F90 | 15 +++++++--------
1 file changed, 7 insertions(+), 8 deletions(-)
diff --git a/src/fourier_mod.F90 b/src/fourier_mod.F90
index 091525d5..3eb22833 100644
--- a/src/fourier_mod.F90
+++ b/src/fourier_mod.F90
@@ -306,7 +306,7 @@ END SUBROUTINE fft1D_plans
! module variable (convolution theorem)
SUBROUTINE poisson_bracket_and_sum( ky_array, kx_array, inv_Ny, inv_Nx, AA_y, AA_x,&
local_nky, total_nkx, F_, G_,&
- ExB_NL_CORRECTION, ExB_NL_factor,sky_ExB,sum_real_)
+ ExB_NL_CORRECTION,ExB_NL_factor,sum_real_)
IMPLICIT NONE
INTEGER, INTENT(IN) :: local_nky,total_nkx
REAL(xp), INTENT(IN) :: inv_Nx, inv_Ny
@@ -318,12 +318,10 @@ END SUBROUTINE fft1D_plans
COMPLEX(xp), DIMENSION(total_nkx,local_nky), &
INTENT(IN) :: ExB_NL_factor
LOGICAL, INTENT(IN) :: ExB_NL_CORRECTION
- REAL(xp),DIMENSION(local_nky), INTENT(IN) :: sky_ExB
real(c_xp_r), pointer, INTENT(INOUT) :: sum_real_(:,:)
! local variables
INTEGER :: ikx,iky
COMPLEX(xp), DIMENSION(total_nkx,local_nky) :: ikxF, ikyG, ikyF, ikxG
- COMPLEX(xp), DIMENSION(total_nkx,local_nky) :: invfactor !TEST
REAL(xp):: kxs, ky
! Build the fields to convolve
@@ -347,10 +345,12 @@ END SUBROUTINE fft1D_plans
ENDIF
! Anti Aliasing
DO iky = 1,local_nky
- ikxF(:,iky) = ikxF(:,iky)*AA_y(iky)*AA_x(:)
- ikyG(:,iky) = ikyG(:,iky)*AA_y(iky)*AA_x(:)
- ikyF(:,iky) = ikyF(:,iky)*AA_y(iky)*AA_x(:)
- ikxG(:,iky) = ikxG(:,iky)*AA_y(iky)*AA_x(:)
+ DO ikx = 1,total_nkx
+ ikxF(ikx,iky) = ikxF(ikx,iky)*AA_y(iky)*AA_x(ikx)
+ ikyG(ikx,iky) = ikyG(ikx,iky)*AA_y(iky)*AA_x(ikx)
+ ikyF(ikx,iky) = ikyF(ikx,iky)*AA_y(iky)*AA_x(ikx)
+ ikxG(ikx,iky) = ikxG(ikx,iky)*AA_y(iky)*AA_x(ikx)
+ ENDDO
ENDDO
!-------------------- First term df/dx x dg/dy --------------------
#ifdef SINGLE_PRECISION
@@ -372,7 +372,6 @@ END SUBROUTINE fft1D_plans
call fftw_mpi_execute_dft_c2r(planb, ikxG, real_data_g)
#endif
sum_real_ = sum_real_ - real_data_f*real_data_g*inv_Ny*inv_Nx
-
END SUBROUTINE poisson_bracket_and_sum
!******************************************************************************!
!******************************************************************************!
--
GitLab