diff --git a/src/modules/mesh/2DCart/moduleMesh2DCart.f90 b/src/modules/mesh/2DCart/moduleMesh2DCart.f90
index 7bc94ce..4bf1f32 100644
--- a/src/modules/mesh/2DCart/moduleMesh2DCart.f90
+++ b/src/modules/mesh/2DCart/moduleMesh2DCart.f90
@@ -556,6 +556,7 @@ MODULE moduleMesh2DCart
 
     !Compute element volume
     PURE SUBROUTINE volumeQuad(self)
+      USE moduleRefParam, ONLY: L_ref
       IMPLICIT NONE
 
       CLASS(meshCell2DCartQuad), INTENT(inout):: self
@@ -573,7 +574,7 @@ MODULE moduleMesh2DCart
       fPsi = self%fPsi(Xi, 4)
 
       !Compute total volume of the cell
-      self%volume  = detJ*4.D0
+      self%volume  = detJ*4.D0/L_ref
       !Compute volume per node
       self%n1%v = self%n1%v + fPsi(1)*self%volume
       self%n2%v = self%n2%v + fPsi(2)*self%volume
diff --git a/src/modules/mesh/inout/gmsh2/moduleMeshInputGmsh2.f90 b/src/modules/mesh/inout/gmsh2/moduleMeshInputGmsh2.f90
index 84b8535..2c02358 100644
--- a/src/modules/mesh/inout/gmsh2/moduleMeshInputGmsh2.f90
+++ b/src/modules/mesh/inout/gmsh2/moduleMeshInputGmsh2.f90
@@ -299,7 +299,7 @@ MODULE moduleMeshInputGmsh2
       ! Adjust node volume at axis
       SELECT CASE(self%geometry)
       CASE("Cyl")
-        DO n  = 1, numNodes
+        DO n  = 1, self%numNodes
           r = self%nodes(n)%obj%getCoordinates()
           IF (r(2) == 0.D0) THEN
             self%nodes(n)%obj%v = self%nodes(n)%obj%v * 2.0D0!2.0D0/3.0D0