debug multiple fault partition

elifo · Sep 30, 2023 · 54e2d1f · 54e2d1f
1 parent d2105bb
commit 54e2d1f
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/src/decompose_mesh/fault_scotch.f90 b/src/decompose_mesh/fault_scotch.f90
@@ -69,7 +69,8 @@ module fault_scotch
   ! parallel faults
 
   ! for parallel faults, chooses a greedy scheme (.true.) or a heuristic balancing (.false.)
-  logical,parameter :: FAULT_BALANCE_GREEDY = .false.
+  !logical,parameter :: FAULT_BALANCE_GREEDY = .false.
+  logical,parameter :: FAULT_BALANCE_GREEDY = .true.
 
   ! adds elements only touching fault by a node/edge to fault partitions
   logical,parameter :: ADD_FAULT_TOUCH_ELEMENTS = .true.
@@ -1484,6 +1485,7 @@ subroutine write_fault_database(IIN_database, iproc, nelmnts, glob2loc_elmnts, &
       print *, '  ispec1 : ', nspec_fault_1
       print *, '  ispec2 : ', nspec_fault_2
       print *, 'Fatal error: Number of fault elements do not coincide. Abort.'
+      print *, 'If the error originates from repartition, try FAULT_BALANCE_GREEDY=.true.' 
       stop
     endif
     write(IIN_database) nspec_fault_1