diff --git a/src/modules/mesh/moduleMesh.f90 b/src/modules/mesh/moduleMesh.f90
index a62cfb2..7c3c4c3 100644
--- a/src/modules/mesh/moduleMesh.f90
+++ b/src/modules/mesh/moduleMesh.f90
@@ -1067,7 +1067,9 @@ MODULE moduleMesh
             END DO
 
             !Divide total momentum exchanged among all the particles of species j
-            normDeltaV = totalDeltaV_ij / REAL(cell%listPart_in(j)%amount) * (self%sp_i%weight*self%sp_i%m)/(self%sp_j%weight*self%sp_j%m)
+            normDeltaV = totalDeltaV_ij / REAL(cell%listPart_in(j)%amount) *       &
+                         (coulombMatrix(k)%sp_i%weight*coulombMatrix(k)%sp_i%m) /  &
+                         (coulombMatrix(k)%sp_j%weight*coulombMatrix(k)%sp_j%m)
             !Loop over particles of species_j
             partTemp => cell%listPart_in(j)%head
             DO WHILE(ASSOCIATED(partTemp))
@@ -1124,6 +1126,8 @@ MODULE moduleMesh
 
         END DO
 
+        DEALLOCATE(densityNodes, velocityNodes, temperatureNodes)
+
       END DO
       !$OMP END DO