[cig-commits] r18532 - seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH

[ampuero at geodynamics.org]

Author: ampuero
Date: 2011-06-02 17:53:58 -0700 (Thu, 02 Jun 2011)
New Revision: 18532

Modified:
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90
Log:
fixed indices of array part in scotch

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90	2011-06-03 00:23:21 UTC (rev 18531)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90	2011-06-03 00:53:58 UTC (rev 18532)
@@ -275,7 +275,7 @@
     allocate(elem_proc_null(nproc_null))
    ! Filling up proc = 0 elements
     nproc_null = 0
-    do i = 1,nelmnts
+    do i = 0,nelmnts-1
       if ( part(i) == 0 ) then
         nproc_null = nproc_null + 1
         elem_proc_null(nproc_null) = i
@@ -294,6 +294,7 @@
       part(elem_proc_null(i)) = ipart
     end do
 
+    deallocate(elem_proc_null)
   endif
  ! call mesh2dual_ncommonnodes_fault(nelmnts, nnodes, nsize, sup_neighbour, &
  !                               elmnts, xadj, adjncy, nnodes_elmnts, &
@@ -330,7 +331,7 @@
 
  ! list of elements per node
  !   nnodes_elmnts(i) = number of elements containing node #i (i=0:nnodes-1)     
- !   nodes_elmnts(nsize*i:nsize*i+nnodes_elmnts(i)-1) = index of elements containing node #i  (k=0:nsize*nnodes-1)
+ !   nodes_elmnts(nsize*i:nsize*i+nnodes_elmnts(i)-1) = index of elements (starting at 0) containing node #i  
   nnodes_elmnts(:) = 0
   nodes_elmnts(:) = 0
   do i = 0, esize*nelmnts-1