diff --git a/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90 b/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90
index f6f2b26..f8e41b3 100644
--- a/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90
+++ b/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90
@@ -524,7 +524,7 @@ MODULE moduleMesh2DCyl
       conv = 1.D0
       XiO  = 0.D0
 
-      DO WHILE(conv > 1.D-2)
+      DO WHILE(conv > 1.D-4)
         dPsi = self%dPsi(XiO, 4)
         pDer = self%partialDer(4, dPsi)
         detJ = self%detJac(pDer)
diff --git a/src/modules/mesh/moduleMesh.f90 b/src/modules/mesh/moduleMesh.f90
index cb252e8..0804059 100644
--- a/src/modules/mesh/moduleMesh.f90
+++ b/src/modules/mesh/moduleMesh.f90
@@ -639,8 +639,8 @@ MODULE moduleMesh
       CLASS(meshCell), INTENT(inout):: self
       CLASS(particle), INTENT(inout), TARGET:: part
       CLASS(meshCell), OPTIONAL, INTENT(in):: oldCell
-      REAL(8):: Xi(1:3)
-      CLASS(meshElement), POINTER:: neighbourElement
+      REAL(8):: Xi(1:3) = 0.D0
+      CLASS(meshElement), POINTER:: neighbourElement => NULL()
       INTEGER:: sp
 
       Xi = self%phy2log(part%r)