diff --git a/src/mesh_marc.f90 b/src/mesh_marc.f90
index 3464f616e..4a5f17674 100644
--- a/src/mesh_marc.f90
+++ b/src/mesh_marc.f90
@@ -888,11 +888,11 @@ subroutine buildCells(thisMesh,elem,connectivity_elem)
   call results_closeJobFile
 #endif  
 
-contains
+  contains
 
-  !--------------------------------------------------------------------------------------------------
-  !> @brief count unique rows (same rows need to be stored consequtively)
-  !--------------------------------------------------------------------------------------------------
+  !------------------------------------------------------------------------------------------------
+  !> @brief count unique rows (same rows need to be stored consecutively)
+  !------------------------------------------------------------------------------------------------
   pure function uniqueRows(A) result(u)
     
     integer, dimension(:,:), intent(in) :: A                                                        !< array, rows need to be sorted
@@ -1223,7 +1223,7 @@ subroutine mesh_build_ipAreas
              forall(n = 1:FE_NcellnodesPerCellface(c)) &
                nodePos(1:3,n) = mesh_cellnode(1:3,mesh_cell(theMesh%elem%cellface(n,f),i,e))
              normal = math_cross(nodePos(1:3,2) - nodePos(1:3,1), &
-                                         nodePos(1:3,3) - nodePos(1:3,1))
+                                 nodePos(1:3,3) - nodePos(1:3,1))
              mesh_ipArea(f,i,e) = norm2(normal)
              mesh_ipAreaNormal(1:3,f,i,e) = normal / norm2(normal)                             ! ensure unit length of area normal
            enddo
@@ -1242,7 +1242,7 @@ subroutine mesh_build_ipAreas
              forall(n = 1:FE_NcellnodesPerCellface(c)) &
                normals(1:3,n) = 0.5_pReal &
                               * math_cross(nodePos(1:3,1+mod(n  ,m)) - nodePos(1:3,n), &
-                                                   nodePos(1:3,1+mod(n+1,m)) - nodePos(1:3,n))
+                                           nodePos(1:3,1+mod(n+1,m)) - nodePos(1:3,n))
              normal = 0.5_pReal * sum(normals,2)
              mesh_ipArea(f,i,e) = norm2(normal)
              mesh_ipAreaNormal(1:3,f,i,e) = normal / norm2(normal)