diff --git a/src/modules/moduleInject.f90 b/src/modules/moduleInject.f90 index 52e26d3..6e7deee 100644 --- a/src/modules/moduleInject.f90 +++ b/src/modules/moduleInject.f90 @@ -273,8 +273,6 @@ MODULE moduleInject do e = 1, self%nEdges edge => self%edges(e)%obj - Xi = edge%centerXi() - if (associated(edge%e1)) then cell => edge%e1 @@ -283,6 +281,8 @@ MODULE moduleInject end if + Xi = edge%centerXi() + ! Projection of EF on the edge normal vector EF_normal = dot_product(cell%gatherElectricField(Xi), edge%normal)