[cig-commits] r14155 - seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2009-02-26 05:51:24 -0800 (Thu, 26 Feb 2009)
New Revision: 14155

Modified:
   seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/check_mesh_quality_CUBIT_Abaqus.f90
   seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/multiply_CUBIT_Abaqus_mesh_by_1000.f90
Log:
added support for Celine's new complex 2D mesh with fractures


Modified: seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/check_mesh_quality_CUBIT_Abaqus.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/check_mesh_quality_CUBIT_Abaqus.f90	2009-02-26 13:45:48 UTC (rev 14154)
+++ seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/check_mesh_quality_CUBIT_Abaqus.f90	2009-02-26 13:51:24 UTC (rev 14155)
@@ -64,7 +64,7 @@
 ! logical, parameter :: IGNORE_OTHER_HEADERS = .false.
 ! double precision, parameter :: delta_t = 5.d-3
 
-  character(len=100), parameter :: cubit_mesh_file = 'eros_complexe_2d_regolite_fractures_in_meters.inp'
+  character(len=100), parameter :: cubit_mesh_file = 'eros_complexe_2d_regolite_fractures_modifie_in_meters.inp'
   integer, parameter :: NPOIN = 57807, NSPEC = 56983, NGNOD = 4
   logical, parameter :: IGNORE_OTHER_HEADERS = .true.
   double precision, parameter :: delta_t = 1.5d-4
@@ -160,7 +160,7 @@
     if(NGNOD == 4) then
 
 !! DK DK ignore other headers for 2D mesh of Eros with fractures, which has multiple material sets
-      if(IGNORE_OTHER_HEADERS .and. (i == 51276 .or. i == 53917)) read(10,*)
+      if(IGNORE_OTHER_HEADERS .and. i == 5709) read(10,*)
 
       read(10,*) iread,ibool(1,i),ibool(2,i),ibool(3,i),ibool(4,i)
     else if(NGNOD == 8) then

Modified: seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/multiply_CUBIT_Abaqus_mesh_by_1000.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/multiply_CUBIT_Abaqus_mesh_by_1000.f90	2009-02-26 13:45:48 UTC (rev 14154)
+++ seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/multiply_CUBIT_Abaqus_mesh_by_1000.f90	2009-02-26 13:51:24 UTC (rev 14155)
@@ -17,7 +17,7 @@
 ! character(len=100), parameter :: cubit_mesh_file = 'HOMOGENE_2D_in_meters.inp'
 ! integer, parameter :: NPOIN = 3882
 
-  character(len=100), parameter :: cubit_mesh_file = 'eros_complexe_2d_regolite_fractures_in_meters.inp'
+  character(len=100), parameter :: cubit_mesh_file = 'eros_complexe_2d_regolite_fractures_modifie_in_meters.inp'
   integer, parameter :: NPOIN = 57807
 
   real, dimension(NPOIN) :: x,y,z