diff --git a/src/fourier_mod.F90 b/src/fourier_mod.F90
index 9eec0d8c7c3c5c685d768ebc05a91b3b2f4ce58b..53dfde72f2c86ef4807c3320e743f8101ae4a6d4 100644
--- a/src/fourier_mod.F90
+++ b/src/fourier_mod.F90
@@ -92,7 +92,8 @@ MODULE fourier
     ENDDO
     call fftw_mpi_execute_dft_c2r(planb, cmpx_data_f, real_data_f)
     call fftw_mpi_execute_dft_c2r(planb, cmpx_data_g, real_data_g)
-    bracket_sum_r = bracket_sum_r + real_data_f*inv_Ny*inv_Nx  * real_data_g*inv_Ny*inv_Nx
+    ! bracket_sum_r = bracket_sum_r + real_data_f*inv_Ny*inv_Nx  * real_data_g*inv_Ny*inv_Nx
+    bracket_sum_r = bracket_sum_r + real_data_f  * real_data_g*inv_Ny*inv_Nx
     ! Second term -df/dy x dg/dx
     DO ikx = ikxs, ikxe
       DO iky = ikys, ikye
@@ -104,7 +105,8 @@ MODULE fourier
     ENDDO
     call fftw_mpi_execute_dft_c2r(planb, cmpx_data_f, real_data_f)
     call fftw_mpi_execute_dft_c2r(planb, cmpx_data_g, real_data_g)
-    bracket_sum_r = bracket_sum_r - real_data_f*inv_Ny*inv_Nx  * real_data_g*inv_Ny*inv_Nx
+    ! bracket_sum_r = bracket_sum_r - real_data_f*inv_Ny*inv_Nx  * real_data_g*inv_Ny*inv_Nx
+    bracket_sum_r = bracket_sum_r - real_data_f  * real_data_g*inv_Ny*inv_Nx
 END SUBROUTINE poisson_bracket_and_sum
 
 !!! Compute the poisson bracket of [F,G] to real space
@@ -122,7 +124,8 @@ SUBROUTINE convolve_and_add( F_, G_)
   ENDDO
   call fftw_mpi_execute_dft_c2r(planb, cmpx_data_f, real_data_f)
   call fftw_mpi_execute_dft_c2r(planb, cmpx_data_g, real_data_g)
-  bracket_sum_r = bracket_sum_r + real_data_f*inv_Ny*inv_Nx  * real_data_g*inv_Ny*inv_Nx
+  ! bracket_sum_r = bracket_sum_r + real_data_f*inv_Ny*inv_Nx  * real_data_g*inv_Ny*inv_Nx
+  bracket_sum_r = bracket_sum_r + real_data_f  * real_data_g*inv_Ny*inv_Nx
 END SUBROUTINE convolve_and_add