[cig-commits] r14219 - seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Wed Mar 4 06:47:17 PST 2009
Author: dkomati1
Date: 2009-03-04 06:47:17 -0800 (Wed, 04 Mar 2009)
New Revision: 14219
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 a new 2D asteroid model from Celine Blitz
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-03-03 23:53:15 UTC (rev 14218)
+++ seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/check_mesh_quality_CUBIT_Abaqus.f90 2009-03-04 14:47:17 UTC (rev 14219)
@@ -64,10 +64,15 @@
! 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_modifie_in_meters.inp'
- integer, parameter :: NPOIN = 57807, NSPEC = 56983, NGNOD = 4
+! 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
+
+ character(len=100), parameter :: cubit_mesh_file = 'REGOLITE_only_no_fractures_2D_in_meters.inp'
+ integer, parameter :: NPOIN = 32536, NSPEC = 31695, NGNOD = 4
logical, parameter :: IGNORE_OTHER_HEADERS = .true.
- double precision, parameter :: delta_t = 1.5d-4
+ double precision, parameter :: delta_t = 0.15e-3
double precision, parameter :: VP_MAX = 3000.d0
@@ -160,8 +165,12 @@
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 == 5709) read(10,*)
+ if(IGNORE_OTHER_HEADERS .and. cubit_mesh_file == 'eros_complexe_2d_regolite_fractures_modifie_in_meters.inp' &
+ .and. i == 5709) read(10,*)
+ if(IGNORE_OTHER_HEADERS .and. cubit_mesh_file == 'REGOLITE_only_no_fractures_2D_in_meters.inp' &
+ .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
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)
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-03-03 23:53:15 UTC (rev 14218)
+++ seismo/3D/SPECFEM3D_SESAME/trunk/UTILS/external_mesh/analyze_CUBIT_Abaqus_mesh/multiply_CUBIT_Abaqus_mesh_by_1000.f90 2009-03-04 14:47:17 UTC (rev 14219)
@@ -17,9 +17,12 @@
! 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_modifie_in_meters.inp'
- integer, parameter :: NPOIN = 57807
+! 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
+
real, dimension(NPOIN) :: x,y,z
integer :: i,iread
More information about the CIG-COMMITS
mailing list