diff --git a/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90 b/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90
index 9278a27..df01782 100644
--- a/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90
+++ b/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90
@@ -580,6 +580,7 @@ MODULE moduleMesh2DCyl
       REAL(8):: dPsi(1:3, 1:4), pDer(1:3, 1:3)
 
       self%volume  = 0.D0
+
       !2D 1 point Gauss Quad Integral
       Xi = 0.D0
       dPsi = self%dPsi(Xi, 4)
@@ -590,17 +591,9 @@ MODULE moduleMesh2DCyl
       !Computes total volume of the cell
       self%volume  = r*detJ*PI8 !4*2*pi
       !Computes volume per node
-      Xi = (/-5.D-1, -5.D-1, 0.D0/)
-      r = self%gatherF(Xi, 4, self%r)
       self%n1%v = self%n1%v + fPsi(1)*self%volume
-      Xi = (/ 5.D-1, -5.D-1, 0.D0/)
-      r = self%gatherF(Xi, 4, self%r)
       self%n2%v = self%n2%v + fPsi(2)*self%volume
-      Xi = (/ 5.D-1,  5.D-1, 0.D0/)
-      r = self%gatherF(Xi, 4, self%r)
       self%n3%v = self%n3%v + fPsi(3)*self%volume
-      Xi = (/-5.D-1,  5.D-1, 0.D0/)
-      r = self%gatherF(Xi, 4, self%r)
       self%n4%v = self%n4%v + fPsi(4)*self%volume
 
     END SUBROUTINE volumeQuad