[cig-commits] r15642 - seismo/3D/SPECFEM3D_SESAME/trunk/decompose_mesh_SCOTCH

[pieyre at geodynamics.org]

Author: pieyre
Date: 2009-09-03 02:37:08 -0700 (Thu, 03 Sep 2009)
New Revision: 15642

Modified:
   seismo/3D/SPECFEM3D_SESAME/trunk/decompose_mesh_SCOTCH/decompose_mesh_SCOTCH.f90
Log:
Changed the order of nodes in elements read from the external file mesh_file, Stacey's conditions validated

Modified: seismo/3D/SPECFEM3D_SESAME/trunk/decompose_mesh_SCOTCH/decompose_mesh_SCOTCH.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/trunk/decompose_mesh_SCOTCH/decompose_mesh_SCOTCH.f90	2009-09-01 18:36:22 UTC (rev 15641)
+++ seismo/3D/SPECFEM3D_SESAME/trunk/decompose_mesh_SCOTCH/decompose_mesh_SCOTCH.f90	2009-09-03 09:37:08 UTC (rev 15642)
@@ -73,8 +73,8 @@
   read(98,*) nspec
   allocate(elmnts(esize,nspec))
    do ispec = 1, nspec
-     read(98,*) num_elmnt, elmnts(1,num_elmnt), elmnts(2,num_elmnt), elmnts(3,num_elmnt), elmnts(4,num_elmnt), &
-          elmnts(5,num_elmnt), elmnts(6,num_elmnt), elmnts(7,num_elmnt), elmnts(8,num_elmnt)
+     read(98,*) num_elmnt, elmnts(5,num_elmnt), elmnts(1,num_elmnt),elmnts(4,num_elmnt), elmnts(8,num_elmnt), &
+          elmnts(6,num_elmnt), elmnts(2,num_elmnt), elmnts(3,num_elmnt), elmnts(7,num_elmnt)
      if((num_elmnt > nspec) .or. (num_elmnt < 1) )  stop "ERROR : Invalid mesh file."
   end do
   close(98)