diff --git a/src/ExB_shear_flow_mod.F90 b/src/ExB_shear_flow_mod.F90
index ec34ae7d0017c4cac0f7d30548546af0e6686338..dadf1aad8e7718cb4613536f01d9bbcf5344ecc2 100644
--- a/src/ExB_shear_flow_mod.F90
+++ b/src/ExB_shear_flow_mod.F90
@@ -128,6 +128,7 @@ CONTAINS
                     sky_ExB_full(iky) = sky_ExB_full(iky) - REAL(iky-1,xp)*deltaky*gamma_E*dt_sub
                 ENDDO
                 ! Shift the arrays if the shear value sky is too high
+                IF(LINEARITY .EQ. 'nonlinear') &
                 CALL Array_shift_ExB_shear_flow
 
                 ! We update the operators and grids