[cig-commits] r14414 - seismo/3D/SPECFEM3D_SESAME/trunk/check_mesh_quality_CUBIT_Abaqus

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2009-03-22 08:36:18 -0700 (Sun, 22 Mar 2009)
New Revision: 14414

Modified:
   seismo/3D/SPECFEM3D_SESAME/trunk/check_mesh_quality_CUBIT_Abaqus/check_mesh_quality_CUBIT_Abaqus.f90
   seismo/3D/SPECFEM3D_SESAME/trunk/check_mesh_quality_CUBIT_Abaqus/multiply_CUBIT_Abaqus_mesh_by_1000.f90
Log:
added support for new 3D model from Celine Blitz


Modified: seismo/3D/SPECFEM3D_SESAME/trunk/check_mesh_quality_CUBIT_Abaqus/check_mesh_quality_CUBIT_Abaqus.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/trunk/check_mesh_quality_CUBIT_Abaqus/check_mesh_quality_CUBIT_Abaqus.f90	2009-03-21 16:57:46 UTC (rev 14413)
+++ seismo/3D/SPECFEM3D_SESAME/trunk/check_mesh_quality_CUBIT_Abaqus/check_mesh_quality_CUBIT_Abaqus.f90	2009-03-22 15:36:18 UTC (rev 14414)
@@ -47,11 +47,17 @@
 ! double precision, parameter :: delta_t = 1.d-3
 ! double precision, parameter :: VP_MAX = 3000.d0
 
-  character(len=100), parameter :: cubit_mesh_file = 'regolite_3D_rego3d_70m_in_meters.inp'
-  integer, parameter :: NPOIN = 4050696, NSPEC = 3410265, NGNOD = 8
-  logical, parameter :: IGNORE_OTHER_HEADERS = .false.
+! character(len=100), parameter :: cubit_mesh_file = 'regolite_3D_rego3d_70m_in_meters.inp'
+! integer, parameter :: NPOIN = 4050696, NSPEC = 3410265, NGNOD = 8
+! logical, parameter :: IGNORE_OTHER_HEADERS = .false.
+! double precision, parameter :: delta_t = 3.d-4
+! double precision, parameter :: VP_MAX = 900.d0 ! because the smallest element is in the regolith layer, not in the bedrock
+
+  character(len=100), parameter :: cubit_mesh_file = 'rego3d_70_disp.inp'
+  integer, parameter :: NPOIN = 5924713, NSPEC = 5797440, NGNOD = 8
+  logical, parameter :: IGNORE_OTHER_HEADERS = .true.
   double precision, parameter :: delta_t = 3.d-4
-  double precision, parameter :: VP_MAX = 900.d0 ! because the smallest element is in the regolith layer, not in the bedrock
+  double precision, parameter :: VP_MAX = 3000.d0 ! 900.d0 ! because the smallest element is in the regolith layer, not in the bedrock
 
 ! character(len=100), parameter :: cubit_mesh_file = 'HOMOGENE_2D_in_meters.inp'
 ! integer, parameter :: NPOIN = 3882, NSPEC = 3744, NGNOD = 4
@@ -183,8 +189,14 @@
                  .and. i == 28429) read(10,*)
 
       read(10,*) iread,ibool(1,i),ibool(2,i),ibool(3,i),ibool(4,i)
+
     else if(NGNOD == 8) then
+
+      if(IGNORE_OTHER_HEADERS .and. cubit_mesh_file == 'rego3d_70_disp.inp' &
+                 .and. i == 252929) read(10,*)
+
       read(10,*) iread,ibool(1,i),ibool(2,i),ibool(3,i),ibool(4,i),ibool(5,i),ibool(6,i),ibool(7,i),ibool(8,i)
+
     endif
 
     if(iread /= i) then

Modified: seismo/3D/SPECFEM3D_SESAME/trunk/check_mesh_quality_CUBIT_Abaqus/multiply_CUBIT_Abaqus_mesh_by_1000.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/trunk/check_mesh_quality_CUBIT_Abaqus/multiply_CUBIT_Abaqus_mesh_by_1000.f90	2009-03-21 16:57:46 UTC (rev 14413)
+++ seismo/3D/SPECFEM3D_SESAME/trunk/check_mesh_quality_CUBIT_Abaqus/multiply_CUBIT_Abaqus_mesh_by_1000.f90	2009-03-22 15:36:18 UTC (rev 14414)
@@ -20,9 +20,12 @@
 ! character(len=100), parameter :: cubit_mesh_file = 'eros_complexe_2d_regolite_fractures_modifie_in_meters.inp'
 ! integer, parameter :: NPOIN = 57807
 
-  character(len=100), parameter :: cubit_mesh_file = 'REGOLITE_only_no_fractures_2D_in_meters.inp'
-  integer, parameter :: NPOIN = 32536
+! character(len=100), parameter :: cubit_mesh_file = 'REGOLITE_only_no_fractures_2D_in_meters.inp'
+! integer, parameter :: NPOIN = 32536
 
+  character(len=100), parameter :: cubit_mesh_file = 'rego3d_70_disp.inp'
+  integer, parameter :: NPOIN = 5924713
+
   real, dimension(NPOIN) :: x,y,z
 
   integer :: i,iread