diff --git a/sources/integrals.F90 b/sources/integrals.F90
index 5a99fbf..e3d2227 100644
--- a/sources/integrals.F90
+++ b/sources/integrals.F90
@@ -743,7 +743,7 @@ module integrals
             tmp(:,1,:) = 0.5d+00 * (pdata%q(iby,nbu:neu,nb ,nb:ne)             &
                                  -  pdata%q(iby,nbl:nel,nb ,nb:ne)) / dh(1)    &
                                  - (pdata%q(ibx,nb :ne ,nb ,nb:ne)             &
-                                  - pdata%quibx,nb :ne ,nbl,nb:ne)) / dh(2)
+                                  - pdata%q(ibx,nb :ne ,nbl,nb:ne)) / dh(2)
 
             inarr(15) = inarr(15) - resistivity * sum(sign(tmp(:,1,:),         &
                                           pdata%q(ibx,nb:ne,nb ,nb:ne))) * dxz