diff --git a/src/modules/init/moduleInput.f90 b/src/modules/init/moduleInput.f90
index b09992a..a54ef24 100644
--- a/src/modules/init/moduleInput.f90
+++ b/src/modules/init/moduleInput.f90
@@ -1012,9 +1012,6 @@ MODULE moduleInput
       if (found) then
         allocate(physicalSurfaces(1:nPhysicalSurfaces))
 
-      else
-        call criticalError('No physical surfaces found in the input file', 'readPhysicalSurfaces')
-
       end if
 
       do ps = 1, nPhysicalSurfaces
@@ -1339,8 +1336,11 @@ MODULE moduleInput
 
       end do
 
-      ! Builds the K matrix for the Particles mesh
-      call mesh%constructGlobalK()
+      if (mesh%dimen > 0) then
+        ! Builds the K matrix for the Particles mesh
+        call mesh%constructGlobalK()
+
+      end if
 
     END SUBROUTINE readGeometry
 
diff --git a/src/modules/mesh/moduleMesh@elements.f90 b/src/modules/mesh/moduleMesh@elements.f90
index b96d737..93ea8ff 100644
--- a/src/modules/mesh/moduleMesh@elements.f90
+++ b/src/modules/mesh/moduleMesh@elements.f90
@@ -183,7 +183,7 @@ submodule(moduleMesh) elements
       !Compute the PLU factorization of K once boundary conditions have been read
       CALL dgetrf(self%numNodes, self%numNodes, self%K, self%numNodes, self%IPIV, info)
       IF (info /= 0) THEN
-        CALL criticalError('Factorization of K matrix failed', 'readBoundaryEM')
+        CALL criticalError('Factorization of K matrix failed', 'constructGlobalK')
 
       END IF