diff --git a/src/solve_EM_fields.F90 b/src/solve_EM_fields.F90
index d28a4aca89bd83c3e024baa31cc9e580e563d9a6..dfbc9d2dedfa97ca32b572316832535855cf0fec 100644
--- a/src/solve_EM_fields.F90
+++ b/src/solve_EM_fields.F90
@@ -69,10 +69,7 @@ CONTAINS
             !!!!!!!!!!!!!!! adiabatic ion model
             ! Candy et al. 2007, rho_i = -q_i/tau_i phi
             IF (ADIAB_I) THEN
-              DO iz = 1,local_nz
-                izi = iz+ngz/2
-                rho(iz) = rho(iz) - 0*q_i/tau_i * phi(iky,ikx,izi)
-              ENDDO
+              ! No ion contribution
             ENDIF
 
             !!!!!!!!!!!!!!! Inverting the poisson equation