diff --git a/src/nonlinear_mod.F90 b/src/nonlinear_mod.F90
index abd6b9b07f12224641b5429a4bd8628b57cb4aa2..4c7f1eafe7a8c48d3ab141ee2ccc69d05c5554ce 100644
--- a/src/nonlinear_mod.F90
+++ b/src/nonlinear_mod.F90
@@ -14,7 +14,7 @@ MODULE nonlinear
   USE prec_const,  ONLY : xp
   USE species,     ONLY : sqrt_tau_o_sigma
   USE time_integration, ONLY : updatetlevel
-  USE ExB_shear_flow,   ONLY : ExB_NL_factor, inv_ExB_NL_factor, sky_ExB
+  USE ExB_shear_flow,   ONLY : ExB_NL_factor, inv_ExB_NL_factor
   use, intrinsic :: iso_c_binding
 
   IMPLICIT NONE
@@ -101,7 +101,7 @@ SUBROUTINE compute_nonlinear
               ! this function adds its result to bracket_sum_r
                 CALL poisson_bracket_and_sum( kyarray,kxarray,inv_Ny,inv_Nx,AA_y,AA_x,&
                                               local_nky,total_nkx,F_cmpx,G_cmpx,&
-                                              ExB_NL_CORRECTION, ExB_NL_factor, sky_ExB, bracket_sum_r)
+                                              ExB_NL_CORRECTION, ExB_NL_factor, bracket_sum_r)
   !-----------!! ELECTROMAGNETIC CONTRIBUTION -sqrt(tau)/sigma*{Sum_s dnjs [sqrt(p+1)Nap+1s + sqrt(p)Nap-1s], Kernel psi}
               IF(EM) THEN
                 ! First convolution terms
@@ -117,7 +117,7 @@ SUBROUTINE compute_nonlinear
                 ! this function adds its result to bracket_sum_r
                 CALL poisson_bracket_and_sum( kyarray,kxarray,inv_Ny,inv_Nx,AA_y,AA_x,&
                                               local_nky,total_nkx,F_cmpx,G_cmpx,&
-                                              ExB_NL_CORRECTION, ExB_NL_factor, sky_ExB, bracket_sum_r)
+                                              ExB_NL_CORRECTION, ExB_NL_factor, bracket_sum_r)
               ENDIF
             ENDDO n
             ! Apply the ExB shearing rate factor before going back to k-space