Abort simulation when PARTICLES=OFF and load balance occurs as this f…

…eature is not yet implemented.
piclas-framework · Dec 5, 2022 · 9e74c81 · 9e74c81
1 parent 388129a
commit 9e74c81
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 2 deletions.
diff --git a/src/loadbalance/loadbalance_tools.f90 b/src/loadbalance/loadbalance_tools.f90
@@ -96,6 +96,9 @@ SUBROUTINE DomainDecomposition()
 IF (PerformLoadBalance) THEN
   nElemsOld     = nElems
   offsetElemOld = offsetElem
+#if ! defined(PARTICLES)
+  CALL CollectiveStop(__STAMP__,'Load balance not implemented for PARTICLES=OFF')
+#endif /*defined(PARTICLES)*/
   IF (myComputeNodeRank.EQ.0) &
     ElemInfoRank_Shared  = ElemInfo_Shared(ELEM_RANK,:)
   CALL BARRIER_AND_SYNC(ElemInfoRank_Shared_Win,MPI_COMM_SHARED)

diff --git a/src/particles/particle_mesh/particle_mesh_readin.f90 b/src/particles/particle_mesh/particle_mesh_readin.f90
@@ -343,7 +343,7 @@ SUBROUTINE ReadMeshNodes()
 #endif  /*USE_MPI*/
 
 ! Reduce NodeCoords if no curved elements are to be used
-ELSE
+ELSE ! .NOT. (useCurveds.OR.NGeo.EQ.1)
   ! every proc needs to allocate the array
   ALLOCATE(NodeInfo(1:nNonUniqueGlobalNodes))
 
@@ -431,7 +431,7 @@ SUBROUTINE ReadMeshNodes()
 
   END ASSOCIATE
 
-END IF
+END IF ! useCurveds.OR.NGeo.EQ.1
 
 ! Update node counters
 offsetNodeID = ElemInfo_Shared(ELEM_FIRSTNODEIND,FirstElemInd) ! hdf5 array starts at 0-> -1