diff --git a/src/modules/mesh/2DCart/moduleMesh2DCart.f90 b/src/modules/mesh/2DCart/moduleMesh2DCart.f90
index dbc8b25..c7308bf 100644
--- a/src/modules/mesh/2DCart/moduleMesh2DCart.f90
+++ b/src/modules/mesh/2DCart/moduleMesh2DCart.f90
@@ -701,6 +701,8 @@ MODULE moduleMesh2DCart
       pDer(2, 1:2) = (/ DOT_PRODUCT(dPsi(1,1:3),self%y(1:3)), &
                         DOT_PRODUCT(dPsi(2,1:3),self%y(1:3)) /)
 
+      pDer(3,3) = 1.D0
+
     END FUNCTION partialDerTria
 
     !Gather electric field at position Xi
diff --git a/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90 b/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90
index ae1eb92..bb94440 100644
--- a/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90
+++ b/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90
@@ -705,8 +705,8 @@ MODULE moduleMesh2DCyl
 
       dPsi = 0.D0
 
-      dPsi(1,:) = (/ -1.D0, 1.D0, 0.D0 /)
-      dPsi(2,:) = (/ -1.D0, 0.D0, 1.D0 /)
+      dPsi(1,1:3) = (/ -1.D0, 1.D0, 0.D0 /)
+      dPsi(2,1:3) = (/ -1.D0, 0.D0, 1.D0 /)
 
     END FUNCTION dPsiTria
 
diff --git a/src/modules/solver/moduleSolver.f90 b/src/modules/solver/moduleSolver.f90
index f85d812..df8d5e8 100644
--- a/src/modules/solver/moduleSolver.f90
+++ b/src/modules/solver/moduleSolver.f90
@@ -140,6 +140,7 @@ MODULE moduleSolver
 
       CASE('ElectrostaticBoltzmann')
         self%solveEM => solveElecFieldBoltzmann
+
       END SELECT
 
     END SUBROUTINE initEM