diff --git a/srcs/BEM/postProcessing.cpp b/srcs/BEM/postProcessing.cpp
index 408d57b0bd40473218dcbc3576f4fcf8a900c28b..7a45cc63409d842bee3bc799a04536f7a788eb64 100644
--- a/srcs/BEM/postProcessing.cpp
+++ b/srcs/BEM/postProcessing.cpp
@@ -428,6 +428,10 @@ void computeElementData(std::vector<std::vector<std::vector<std::vector<double>>
                     }
                 }
 
+                currentPhi += localCurrentPhi;
+                gradPhiX += localGradPhiX;
+                gradPhiY += localGradPhiY;
+
                 phi[i][j][k].resize(1, currentPhi);
                 electricField[i][j][k].resize(3);
                 electricField[i][j][k][0] = -gradPhiX;