diff --git a/src/numerics_mod.F90 b/src/numerics_mod.F90
index 1c769cae61318c30dd5f275f944b94934a8e26a8..5a29eda09542cbf67181fcdbf693dba5dc148d67 100644
--- a/src/numerics_mod.F90
+++ b/src/numerics_mod.F90
@@ -158,7 +158,7 @@ SUBROUTINE evaluate_poisson_op
     ELSE ! Adiabatic model
       pol_e = qe2_taue - 1._dp
     ENDIF
-    inv_poisson_op(iky, ikx, iz) =  1._dp/(qe2_taue + qi2_taui - pol_i - pol_e)
+    inv_poisson_op(iky, ikx, iz) =  1._dp/(qi2_taui - pol_i + qe2_taue - pol_e)
     inv_pol_ion   (iky, ikx, iz) =  1._dp/(qi2_taui - pol_i)
   ENDIF
   END DO zloop