diff --git a/src/modules/mesh/2DCart/moduleMesh2DCart.f90 b/src/modules/mesh/2DCart/moduleMesh2DCart.f90
index c7308bf..eb834e5 100644
--- a/src/modules/mesh/2DCart/moduleMesh2DCart.f90
+++ b/src/modules/mesh/2DCart/moduleMesh2DCart.f90
@@ -520,7 +520,7 @@ MODULE moduleMesh2DCart
         fPsi = self%fPsi(XiO, 4)
         f(1:2) = (/ DOT_PRODUCT(fPsi,self%x), &
                     DOT_PRODUCT(fPsi,self%y) /) - r(1:2)
-        Xi = XiO - MATMUL(invJ, f)/detJ
+        Xi = XiO - MATMUL(transpose(invJ), f)/detJ
         conv = MAXVAL(DABS(Xi-XiO),1)
         XiO = Xi
 
@@ -831,14 +831,14 @@ MODULE moduleMesh2DCart
       REAL(8):: deltaR(1:3)
       REAL(8):: dPsi(1:3, 1:3)
       REAL(8):: pDer(1:3, 1:3)
-      REAL(8):: invJ(1:3,  1:3), detJ
+      REAL(8):: invJ(1:3, 1:3), detJ
 
       !Direct method to convert coordinates
       Xi     = 0.D0
       deltaR = (/ r(1) - self%x(1), r(2) - self%y(1), 0.D0 /)
       dPsi   = self%dPsi(Xi, 3)
       pDer   = self%partialDer(3, dPsi)
-      invJ   = self%invJac(pDer)
+      invJ   = transpose(self%invJac(pDer))
       detJ   = self%detJac(pDer)
       Xi     = MATMUL(invJ,deltaR)/detJ
 
@@ -915,8 +915,8 @@ MODULE moduleMesh2DCart
 
       invJ = 0.D0
 
-      invJ(1, 1:2) = (/  pDer(2,2), -pDer(1,2) /)
-      invJ(2, 1:2) = (/ -pDer(2,1),  pDer(1,1) /)
+      invJ(1, 1:2) = (/  pDer(2,2), -pDer(2,1) /)
+      invJ(2, 1:2) = (/ -pDer(1,2),  pDer(1,1) /)
       invJ(3, 3)   = 1.D0
 
     END FUNCTION invJ2DCart