[cig-commits] [commit] devel, master: converted block_mesh.py to symbolic links to the GEOCUBIT library (c8fd15a)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 16:54:51 PDT 2014
Repository : https://github.com/geodynamics/specfem3d
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d/compare/6026e367984905ab133865f62fa6293b343759b9...47f703851338234f96397e7da9fbff63d8178b8a
>---------------------------------------------------------------
commit c8fd15ab35d05fe93497fa99ab27e6474dcaa384
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date: Tue Jun 3 10:11:42 2014 +0200
converted block_mesh.py to symbolic links to the GEOCUBIT library
>---------------------------------------------------------------
c8fd15ab35d05fe93497fa99ab27e6474dcaa384
.../block_mesh.py | 60 +--------------------
.../block_mesh.py | 60 +--------------------
.../block_mesh.py | 60 +--------------------
.../block_mesh.py | 60 +--------------------
.../block_mesh.py | 60 +--------------------
.../block_mesh.py | 60 +--------------------
.../block_mesh.py | 61 +---------------------
.../block_mesh.py | 60 +--------------------
homogeneous_poroelastic/block_mesh.py | 49 +----------------
9 files changed, 9 insertions(+), 521 deletions(-)
diff --git a/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/block_mesh.py b/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/block_mesh.py
deleted file mode 100755
index f18426a..0000000
--- a/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/block_mesh.py
+++ /dev/null
@@ -1,59 +0,0 @@
-#!/usr/bin/env python
-
-import cubit
-import boundary_definition
-import cubit2specfem3d
-
-import os
-import sys
-
-cubit.cmd('reset')
-cubit.cmd('brick x 134000 y 134000 z 60000')
-cubit.cmd('volume 1 move x 67000 y 67000 z -30000')
-
-
-# Meshing the volumes
-elementsize = 3750.0
-
-cubit.cmd('volume 1 size '+str(elementsize))
-cubit.cmd('mesh volume 1')
-
-
-#### End of meshing
-
-###### This is boundary_definition.py of GEOCUBIT
-#..... which extracts the bounding faces and defines them into blocks
-boundary_definition.entities=['face']
-boundary_definition.define_bc(boundary_definition.entities,parallel=True)
-
-#### Define material properties for the 3 volumes ################
-cubit.cmd('#### DEFINE MATERIAL PROPERTIES #######################')
-
-
-cubit.cmd('block 1 name "elastic 1" ') # elastic material region
-cubit.cmd('block 1 attribute count 6')
-cubit.cmd('block 1 attribute index 1 1') # flag for material: 1 for 1. material
-cubit.cmd('block 1 attribute index 2 2800') # vp
-cubit.cmd('block 1 attribute index 3 1500') # vs
-cubit.cmd('block 1 attribute index 4 2300') # rho
-cubit.cmd('block 1 attribute index 5 9000.0') # Qmu
-cubit.cmd('block 1 attribute index 6 0 ') # anisotropy_flag
-
-
-#cubit.cmd('block 1 name "acoustic 1" ') # acoustic material region
-#cubit.cmd('block 1 attribute count 4')
-#cubit.cmd('block 1 attribute index 1 1 ') # material 1
-#cubit.cmd('block 1 attribute index 2 1480 ') # vp
-#cubit.cmd('block 1 attribute index 3 0 ') # vs
-#cubit.cmd('block 1 attribute index 4 1028 ') # rho (ocean salt water density:
-
-
-cubit.cmd('export mesh "top.e" dimension 3 overwrite')
-cubit.cmd('save as "meshing.cub" overwrite')
-
-#### Export to SPECFEM3D format using cubit2specfem3d.py of GEOCUBIT
-
-os.system('mkdir -p MESH')
-cubit2specfem3d.export2SPECFEM3D('MESH')
-
-# all files needed by SCOTCH are now in directory MESH
diff --git a/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/block_mesh.py b/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/block_mesh.py
new file mode 120000
index 0000000..b7f9fdf
--- /dev/null
+++ b/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/block_mesh.py
@@ -0,0 +1 @@
+../CUBIT_GEOCUBIT/geocubitlib/block_mesh.py
\ No newline at end of file
diff --git a/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/block_mesh.py b/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/block_mesh.py
deleted file mode 100755
index f18426a..0000000
--- a/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/block_mesh.py
+++ /dev/null
@@ -1,59 +0,0 @@
-#!/usr/bin/env python
-
-import cubit
-import boundary_definition
-import cubit2specfem3d
-
-import os
-import sys
-
-cubit.cmd('reset')
-cubit.cmd('brick x 134000 y 134000 z 60000')
-cubit.cmd('volume 1 move x 67000 y 67000 z -30000')
-
-
-# Meshing the volumes
-elementsize = 3750.0
-
-cubit.cmd('volume 1 size '+str(elementsize))
-cubit.cmd('mesh volume 1')
-
-
-#### End of meshing
-
-###### This is boundary_definition.py of GEOCUBIT
-#..... which extracts the bounding faces and defines them into blocks
-boundary_definition.entities=['face']
-boundary_definition.define_bc(boundary_definition.entities,parallel=True)
-
-#### Define material properties for the 3 volumes ################
-cubit.cmd('#### DEFINE MATERIAL PROPERTIES #######################')
-
-
-cubit.cmd('block 1 name "elastic 1" ') # elastic material region
-cubit.cmd('block 1 attribute count 6')
-cubit.cmd('block 1 attribute index 1 1') # flag for material: 1 for 1. material
-cubit.cmd('block 1 attribute index 2 2800') # vp
-cubit.cmd('block 1 attribute index 3 1500') # vs
-cubit.cmd('block 1 attribute index 4 2300') # rho
-cubit.cmd('block 1 attribute index 5 9000.0') # Qmu
-cubit.cmd('block 1 attribute index 6 0 ') # anisotropy_flag
-
-
-#cubit.cmd('block 1 name "acoustic 1" ') # acoustic material region
-#cubit.cmd('block 1 attribute count 4')
-#cubit.cmd('block 1 attribute index 1 1 ') # material 1
-#cubit.cmd('block 1 attribute index 2 1480 ') # vp
-#cubit.cmd('block 1 attribute index 3 0 ') # vs
-#cubit.cmd('block 1 attribute index 4 1028 ') # rho (ocean salt water density:
-
-
-cubit.cmd('export mesh "top.e" dimension 3 overwrite')
-cubit.cmd('save as "meshing.cub" overwrite')
-
-#### Export to SPECFEM3D format using cubit2specfem3d.py of GEOCUBIT
-
-os.system('mkdir -p MESH')
-cubit2specfem3d.export2SPECFEM3D('MESH')
-
-# all files needed by SCOTCH are now in directory MESH
diff --git a/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/block_mesh.py b/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/block_mesh.py
new file mode 120000
index 0000000..b7f9fdf
--- /dev/null
+++ b/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/block_mesh.py
@@ -0,0 +1 @@
+../CUBIT_GEOCUBIT/geocubitlib/block_mesh.py
\ No newline at end of file
diff --git a/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_5sides/block_mesh.py b/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_5sides/block_mesh.py
deleted file mode 100755
index f18426a..0000000
--- a/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_5sides/block_mesh.py
+++ /dev/null
@@ -1,59 +0,0 @@
-#!/usr/bin/env python
-
-import cubit
-import boundary_definition
-import cubit2specfem3d
-
-import os
-import sys
-
-cubit.cmd('reset')
-cubit.cmd('brick x 134000 y 134000 z 60000')
-cubit.cmd('volume 1 move x 67000 y 67000 z -30000')
-
-
-# Meshing the volumes
-elementsize = 3750.0
-
-cubit.cmd('volume 1 size '+str(elementsize))
-cubit.cmd('mesh volume 1')
-
-
-#### End of meshing
-
-###### This is boundary_definition.py of GEOCUBIT
-#..... which extracts the bounding faces and defines them into blocks
-boundary_definition.entities=['face']
-boundary_definition.define_bc(boundary_definition.entities,parallel=True)
-
-#### Define material properties for the 3 volumes ################
-cubit.cmd('#### DEFINE MATERIAL PROPERTIES #######################')
-
-
-cubit.cmd('block 1 name "elastic 1" ') # elastic material region
-cubit.cmd('block 1 attribute count 6')
-cubit.cmd('block 1 attribute index 1 1') # flag for material: 1 for 1. material
-cubit.cmd('block 1 attribute index 2 2800') # vp
-cubit.cmd('block 1 attribute index 3 1500') # vs
-cubit.cmd('block 1 attribute index 4 2300') # rho
-cubit.cmd('block 1 attribute index 5 9000.0') # Qmu
-cubit.cmd('block 1 attribute index 6 0 ') # anisotropy_flag
-
-
-#cubit.cmd('block 1 name "acoustic 1" ') # acoustic material region
-#cubit.cmd('block 1 attribute count 4')
-#cubit.cmd('block 1 attribute index 1 1 ') # material 1
-#cubit.cmd('block 1 attribute index 2 1480 ') # vp
-#cubit.cmd('block 1 attribute index 3 0 ') # vs
-#cubit.cmd('block 1 attribute index 4 1028 ') # rho (ocean salt water density:
-
-
-cubit.cmd('export mesh "top.e" dimension 3 overwrite')
-cubit.cmd('save as "meshing.cub" overwrite')
-
-#### Export to SPECFEM3D format using cubit2specfem3d.py of GEOCUBIT
-
-os.system('mkdir -p MESH')
-cubit2specfem3d.export2SPECFEM3D('MESH')
-
-# all files needed by SCOTCH are now in directory MESH
diff --git a/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_5sides/block_mesh.py b/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_5sides/block_mesh.py
new file mode 120000
index 0000000..b7f9fdf
--- /dev/null
+++ b/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_5sides/block_mesh.py
@@ -0,0 +1 @@
+../CUBIT_GEOCUBIT/geocubitlib/block_mesh.py
\ No newline at end of file
diff --git a/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_6sides/block_mesh.py b/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_6sides/block_mesh.py
deleted file mode 100755
index f18426a..0000000
--- a/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_6sides/block_mesh.py
+++ /dev/null
@@ -1,59 +0,0 @@
-#!/usr/bin/env python
-
-import cubit
-import boundary_definition
-import cubit2specfem3d
-
-import os
-import sys
-
-cubit.cmd('reset')
-cubit.cmd('brick x 134000 y 134000 z 60000')
-cubit.cmd('volume 1 move x 67000 y 67000 z -30000')
-
-
-# Meshing the volumes
-elementsize = 3750.0
-
-cubit.cmd('volume 1 size '+str(elementsize))
-cubit.cmd('mesh volume 1')
-
-
-#### End of meshing
-
-###### This is boundary_definition.py of GEOCUBIT
-#..... which extracts the bounding faces and defines them into blocks
-boundary_definition.entities=['face']
-boundary_definition.define_bc(boundary_definition.entities,parallel=True)
-
-#### Define material properties for the 3 volumes ################
-cubit.cmd('#### DEFINE MATERIAL PROPERTIES #######################')
-
-
-cubit.cmd('block 1 name "elastic 1" ') # elastic material region
-cubit.cmd('block 1 attribute count 6')
-cubit.cmd('block 1 attribute index 1 1') # flag for material: 1 for 1. material
-cubit.cmd('block 1 attribute index 2 2800') # vp
-cubit.cmd('block 1 attribute index 3 1500') # vs
-cubit.cmd('block 1 attribute index 4 2300') # rho
-cubit.cmd('block 1 attribute index 5 9000.0') # Qmu
-cubit.cmd('block 1 attribute index 6 0 ') # anisotropy_flag
-
-
-#cubit.cmd('block 1 name "acoustic 1" ') # acoustic material region
-#cubit.cmd('block 1 attribute count 4')
-#cubit.cmd('block 1 attribute index 1 1 ') # material 1
-#cubit.cmd('block 1 attribute index 2 1480 ') # vp
-#cubit.cmd('block 1 attribute index 3 0 ') # vs
-#cubit.cmd('block 1 attribute index 4 1028 ') # rho (ocean salt water density:
-
-
-cubit.cmd('export mesh "top.e" dimension 3 overwrite')
-cubit.cmd('save as "meshing.cub" overwrite')
-
-#### Export to SPECFEM3D format using cubit2specfem3d.py of GEOCUBIT
-
-os.system('mkdir -p MESH')
-cubit2specfem3d.export2SPECFEM3D('MESH')
-
-# all files needed by SCOTCH are now in directory MESH
diff --git a/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_6sides/block_mesh.py b/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_6sides/block_mesh.py
new file mode 120000
index 0000000..b7f9fdf
--- /dev/null
+++ b/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_6sides/block_mesh.py
@@ -0,0 +1 @@
+../CUBIT_GEOCUBIT/geocubitlib/block_mesh.py
\ No newline at end of file
diff --git a/CPML_examples/homogeneous_halfspace_HEX8_elastic_absorbing_CPML_5sides/block_mesh.py b/CPML_examples/homogeneous_halfspace_HEX8_elastic_absorbing_CPML_5sides/block_mesh.py
deleted file mode 100755
index f18426a..0000000
--- a/CPML_examples/homogeneous_halfspace_HEX8_elastic_absorbing_CPML_5sides/block_mesh.py
+++ /dev/null
@@ -1,59 +0,0 @@
-#!/usr/bin/env python
-
-import cubit
-import boundary_definition
-import cubit2specfem3d
-
-import os
-import sys
-
-cubit.cmd('reset')
-cubit.cmd('brick x 134000 y 134000 z 60000')
-cubit.cmd('volume 1 move x 67000 y 67000 z -30000')
-
-
-# Meshing the volumes
-elementsize = 3750.0
-
-cubit.cmd('volume 1 size '+str(elementsize))
-cubit.cmd('mesh volume 1')
-
-
-#### End of meshing
-
-###### This is boundary_definition.py of GEOCUBIT
-#..... which extracts the bounding faces and defines them into blocks
-boundary_definition.entities=['face']
-boundary_definition.define_bc(boundary_definition.entities,parallel=True)
-
-#### Define material properties for the 3 volumes ################
-cubit.cmd('#### DEFINE MATERIAL PROPERTIES #######################')
-
-
-cubit.cmd('block 1 name "elastic 1" ') # elastic material region
-cubit.cmd('block 1 attribute count 6')
-cubit.cmd('block 1 attribute index 1 1') # flag for material: 1 for 1. material
-cubit.cmd('block 1 attribute index 2 2800') # vp
-cubit.cmd('block 1 attribute index 3 1500') # vs
-cubit.cmd('block 1 attribute index 4 2300') # rho
-cubit.cmd('block 1 attribute index 5 9000.0') # Qmu
-cubit.cmd('block 1 attribute index 6 0 ') # anisotropy_flag
-
-
-#cubit.cmd('block 1 name "acoustic 1" ') # acoustic material region
-#cubit.cmd('block 1 attribute count 4')
-#cubit.cmd('block 1 attribute index 1 1 ') # material 1
-#cubit.cmd('block 1 attribute index 2 1480 ') # vp
-#cubit.cmd('block 1 attribute index 3 0 ') # vs
-#cubit.cmd('block 1 attribute index 4 1028 ') # rho (ocean salt water density:
-
-
-cubit.cmd('export mesh "top.e" dimension 3 overwrite')
-cubit.cmd('save as "meshing.cub" overwrite')
-
-#### Export to SPECFEM3D format using cubit2specfem3d.py of GEOCUBIT
-
-os.system('mkdir -p MESH')
-cubit2specfem3d.export2SPECFEM3D('MESH')
-
-# all files needed by SCOTCH are now in directory MESH
diff --git a/CPML_examples/homogeneous_halfspace_HEX8_elastic_absorbing_CPML_5sides/block_mesh.py b/CPML_examples/homogeneous_halfspace_HEX8_elastic_absorbing_CPML_5sides/block_mesh.py
new file mode 120000
index 0000000..b7f9fdf
--- /dev/null
+++ b/CPML_examples/homogeneous_halfspace_HEX8_elastic_absorbing_CPML_5sides/block_mesh.py
@@ -0,0 +1 @@
+../CUBIT_GEOCUBIT/geocubitlib/block_mesh.py
\ No newline at end of file
diff --git a/CPML_examples/homogeneous_halfspace_HEX8_elastic_absorbing_CPML_6sides/block_mesh.py b/CPML_examples/homogeneous_halfspace_HEX8_elastic_absorbing_CPML_6sides/block_mesh.py
deleted file mode 100755
index f18426a..0000000
--- a/CPML_examples/homogeneous_halfspace_HEX8_elastic_absorbing_CPML_6sides/block_mesh.py
+++ /dev/null
@@ -1,59 +0,0 @@
-#!/usr/bin/env python
-
-import cubit
-import boundary_definition
-import cubit2specfem3d
-
-import os
-import sys
-
-cubit.cmd('reset')
-cubit.cmd('brick x 134000 y 134000 z 60000')
-cubit.cmd('volume 1 move x 67000 y 67000 z -30000')
-
-
-# Meshing the volumes
-elementsize = 3750.0
-
-cubit.cmd('volume 1 size '+str(elementsize))
-cubit.cmd('mesh volume 1')
-
-
-#### End of meshing
-
-###### This is boundary_definition.py of GEOCUBIT
-#..... which extracts the bounding faces and defines them into blocks
-boundary_definition.entities=['face']
-boundary_definition.define_bc(boundary_definition.entities,parallel=True)
-
-#### Define material properties for the 3 volumes ################
-cubit.cmd('#### DEFINE MATERIAL PROPERTIES #######################')
-
-
-cubit.cmd('block 1 name "elastic 1" ') # elastic material region
-cubit.cmd('block 1 attribute count 6')
-cubit.cmd('block 1 attribute index 1 1') # flag for material: 1 for 1. material
-cubit.cmd('block 1 attribute index 2 2800') # vp
-cubit.cmd('block 1 attribute index 3 1500') # vs
-cubit.cmd('block 1 attribute index 4 2300') # rho
-cubit.cmd('block 1 attribute index 5 9000.0') # Qmu
-cubit.cmd('block 1 attribute index 6 0 ') # anisotropy_flag
-
-
-#cubit.cmd('block 1 name "acoustic 1" ') # acoustic material region
-#cubit.cmd('block 1 attribute count 4')
-#cubit.cmd('block 1 attribute index 1 1 ') # material 1
-#cubit.cmd('block 1 attribute index 2 1480 ') # vp
-#cubit.cmd('block 1 attribute index 3 0 ') # vs
-#cubit.cmd('block 1 attribute index 4 1028 ') # rho (ocean salt water density:
-
-
-cubit.cmd('export mesh "top.e" dimension 3 overwrite')
-cubit.cmd('save as "meshing.cub" overwrite')
-
-#### Export to SPECFEM3D format using cubit2specfem3d.py of GEOCUBIT
-
-os.system('mkdir -p MESH')
-cubit2specfem3d.export2SPECFEM3D('MESH')
-
-# all files needed by SCOTCH are now in directory MESH
diff --git a/CPML_examples/homogeneous_halfspace_HEX8_elastic_absorbing_CPML_6sides/block_mesh.py b/CPML_examples/homogeneous_halfspace_HEX8_elastic_absorbing_CPML_6sides/block_mesh.py
new file mode 120000
index 0000000..b7f9fdf
--- /dev/null
+++ b/CPML_examples/homogeneous_halfspace_HEX8_elastic_absorbing_CPML_6sides/block_mesh.py
@@ -0,0 +1 @@
+../CUBIT_GEOCUBIT/geocubitlib/block_mesh.py
\ No newline at end of file
diff --git a/homogeneous_halfspace_HEX27_elastic_no_absorbing/block_mesh.py b/homogeneous_halfspace_HEX27_elastic_no_absorbing/block_mesh.py
deleted file mode 100644
index 6e849c5..0000000
--- a/homogeneous_halfspace_HEX27_elastic_no_absorbing/block_mesh.py
+++ /dev/null
@@ -1,60 +0,0 @@
-#!/usr/bin/env python
-
-import cubit
-import boundary_definition
-import cubit2specfem3d
-
-import os
-import sys
-
-cubit.cmd('reset')
-cubit.cmd('brick x 134000 y 134000 z 60000')
-cubit.cmd('volume 1 move x 67000 y 67000 z -30000')
-
-
-# Meshing the volumes
-elementsize = 3750.0
-
-cubit.cmd('volume 1 size '+str(elementsize))
-cubit.cmd('mesh volume 1')
-
-
-#### End of meshing
-
-###### This is boundary_definition.py of GEOCUBIT
-#..... which extracts the bounding faces and defines them into blocks
-boundary_definition.entities=['face']
-boundary_definition.define_bc(boundary_definition.entities,parallel=True)
-
-#### Define material properties for the 3 volumes ################
-cubit.cmd('#### DEFINE MATERIAL PROPERTIES #######################')
-
-
-cubit.cmd('block 1 name "elastic 1" ') # elastic material region
-cubit.cmd('block 1 attribute count 6')
-cubit.cmd('block 1 attribute index 1 1') # flag for material: 1 for 1. material
-cubit.cmd('block 1 attribute index 2 2800') # vp
-cubit.cmd('block 1 attribute index 3 1500') # vs
-cubit.cmd('block 1 attribute index 4 2300') # rho
-cubit.cmd('block 1 attribute index 5 9000.0') # Qmu
-cubit.cmd('block 1 attribute index 6 0 ') # anisotropy_flag
-
-
-#cubit.cmd('block 1 name "acoustic 1" ') # acoustic material region
-#cubit.cmd('block 1 attribute count 4')
-#cubit.cmd('block 1 attribute index 1 1 ') # material 1
-#cubit.cmd('block 1 attribute index 2 1480 ') # vp
-#cubit.cmd('block 1 attribute index 3 0 ') # vs
-#cubit.cmd('block 1 attribute index 4 1028 ') # rho (ocean salt water density:
-
-
-cubit.cmd('export mesh "top.e" dimension 3 overwrite')
-cubit.cmd('save as "meshing.cub" overwrite')
-
-#### Export to SPECFEM3D format using cubit2specfem3d.py of GEOCUBIT
-
-os.system('mkdir -p MESH')
-cubit2specfem3d.export2SPECFEM3D('MESH',hex27=True)
-
-# all files needed by SCOTCH are now in directory MESH
-print os.getcwd()
diff --git a/homogeneous_halfspace_HEX27_elastic_no_absorbing/block_mesh.py b/homogeneous_halfspace_HEX27_elastic_no_absorbing/block_mesh.py
new file mode 120000
index 0000000..b7f9fdf
--- /dev/null
+++ b/homogeneous_halfspace_HEX27_elastic_no_absorbing/block_mesh.py
@@ -0,0 +1 @@
+../CUBIT_GEOCUBIT/geocubitlib/block_mesh.py
\ No newline at end of file
diff --git a/homogeneous_halfspace_HEX8_elastic_no_absorbing/block_mesh.py b/homogeneous_halfspace_HEX8_elastic_no_absorbing/block_mesh.py
deleted file mode 100755
index f18426a..0000000
--- a/homogeneous_halfspace_HEX8_elastic_no_absorbing/block_mesh.py
+++ /dev/null
@@ -1,59 +0,0 @@
-#!/usr/bin/env python
-
-import cubit
-import boundary_definition
-import cubit2specfem3d
-
-import os
-import sys
-
-cubit.cmd('reset')
-cubit.cmd('brick x 134000 y 134000 z 60000')
-cubit.cmd('volume 1 move x 67000 y 67000 z -30000')
-
-
-# Meshing the volumes
-elementsize = 3750.0
-
-cubit.cmd('volume 1 size '+str(elementsize))
-cubit.cmd('mesh volume 1')
-
-
-#### End of meshing
-
-###### This is boundary_definition.py of GEOCUBIT
-#..... which extracts the bounding faces and defines them into blocks
-boundary_definition.entities=['face']
-boundary_definition.define_bc(boundary_definition.entities,parallel=True)
-
-#### Define material properties for the 3 volumes ################
-cubit.cmd('#### DEFINE MATERIAL PROPERTIES #######################')
-
-
-cubit.cmd('block 1 name "elastic 1" ') # elastic material region
-cubit.cmd('block 1 attribute count 6')
-cubit.cmd('block 1 attribute index 1 1') # flag for material: 1 for 1. material
-cubit.cmd('block 1 attribute index 2 2800') # vp
-cubit.cmd('block 1 attribute index 3 1500') # vs
-cubit.cmd('block 1 attribute index 4 2300') # rho
-cubit.cmd('block 1 attribute index 5 9000.0') # Qmu
-cubit.cmd('block 1 attribute index 6 0 ') # anisotropy_flag
-
-
-#cubit.cmd('block 1 name "acoustic 1" ') # acoustic material region
-#cubit.cmd('block 1 attribute count 4')
-#cubit.cmd('block 1 attribute index 1 1 ') # material 1
-#cubit.cmd('block 1 attribute index 2 1480 ') # vp
-#cubit.cmd('block 1 attribute index 3 0 ') # vs
-#cubit.cmd('block 1 attribute index 4 1028 ') # rho (ocean salt water density:
-
-
-cubit.cmd('export mesh "top.e" dimension 3 overwrite')
-cubit.cmd('save as "meshing.cub" overwrite')
-
-#### Export to SPECFEM3D format using cubit2specfem3d.py of GEOCUBIT
-
-os.system('mkdir -p MESH')
-cubit2specfem3d.export2SPECFEM3D('MESH')
-
-# all files needed by SCOTCH are now in directory MESH
diff --git a/homogeneous_halfspace_HEX8_elastic_no_absorbing/block_mesh.py b/homogeneous_halfspace_HEX8_elastic_no_absorbing/block_mesh.py
new file mode 120000
index 0000000..b7f9fdf
--- /dev/null
+++ b/homogeneous_halfspace_HEX8_elastic_no_absorbing/block_mesh.py
@@ -0,0 +1 @@
+../CUBIT_GEOCUBIT/geocubitlib/block_mesh.py
\ No newline at end of file
diff --git a/homogeneous_poroelastic/block_mesh.py b/homogeneous_poroelastic/block_mesh.py
deleted file mode 100755
index 8945e17..0000000
--- a/homogeneous_poroelastic/block_mesh.py
+++ /dev/null
@@ -1,48 +0,0 @@
-#!/usr/bin/env python
-
-import cubit
-import boundary_definition
-import cubit2specfem3d
-
-import os
-import sys
-
-cubit.cmd('reset')
-cubit.cmd('brick x 134000 y 134000 z 60000')
-cubit.cmd('volume 1 move x 67000 y 67000 z -30000')
-
-
-# Meshing the volumes
-elementsize = 3750.0
-
-cubit.cmd('volume 1 size '+str(elementsize))
-cubit.cmd('mesh volume 1')
-
-
-#### End of meshing
-
-###### This is boundary_definition.py of GEOCUBIT
-#..... which extracts the bounding faces and defines them into blocks
-boundary_definition.entities=['face']
-boundary_definition.define_bc(boundary_definition.entities,parallel=True)
-
-#### Define material properties for the volume ################
-cubit.cmd('#### DEFINE MATERIAL PROPERTIES #######################')
-
-
-cubit.cmd('block 1 name "poroelastic 1" ') # porous material region
-cubit.cmd('block 1 attribute count 4')
-cubit.cmd('block 1 attribute index 1 1') # flag for material: 1 for 1. material
-cubit.cmd('block 1 attribute index 2 3677 ') # vp-dummy
-cubit.cmd('block 1 attribute index 3 2378 ') # vs-dummy
-cubit.cmd('block 1 attribute index 4 1908 ') # rho-dummy
-
-cubit.cmd('export mesh "top.e" dimension 3 overwrite')
-cubit.cmd('save as "meshing.cub" overwrite')
-
-#### Export to SPECFEM3D format using cubit2specfem3d.py of GEOCUBIT
-
-os.system('mkdir -p MESH')
-cubit2specfem3d.export2SPECFEM3D('MESH')
-
-# all files needed by SCOTCH are now in directory MESH
diff --git a/homogeneous_poroelastic/block_mesh.py b/homogeneous_poroelastic/block_mesh.py
new file mode 120000
index 0000000..b7f9fdf
--- /dev/null
+++ b/homogeneous_poroelastic/block_mesh.py
@@ -0,0 +1 @@
+../CUBIT_GEOCUBIT/geocubitlib/block_mesh.py
\ No newline at end of file
More information about the CIG-COMMITS
mailing list