[cig-commits] [commit] devel,master: updated waterlayered_halfspace example; did not delete the old .py scripts (ffb4675)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem3d

On branches: devel,master
Link       : https://github.com/geodynamics/specfem3d/compare/6026e367984905ab133865f62fa6293b343759b9...47f703851338234f96397e7da9fbff63d8178b8a

>---------------------------------------------------------------

commit ffb4675f57099fe244afade5939e84a662f28f99
Author: Carl Tape <carltape at gi.alaska.edu>
Date:   Fri Sep 27 07:17:41 2013 +0000

    updated waterlayered_halfspace example; did not delete the old .py scripts


>---------------------------------------------------------------

ffb4675f57099fe244afade5939e84a662f28f99
 .../MESH/nummaterial_velocity_file.reference       |  3 ++
 ...terial_velocity_file.reference.verticalboundary |  2 ++
 waterlayered_halfspace/boundary_definition.py      |  1 -
 waterlayered_halfspace/cubit2specfem3d.py          |  1 -
 .../make_mesh_waterlater_verticalboundary.py       | 33 ++++++++++++++++++++++
 .../make_mesh_waterlayer.sh                        |  6 ++--
 .../waterlayer_only.cfg                            |  0
 7 files changed, 41 insertions(+), 5 deletions(-)

diff --git a/waterlayered_halfspace/MESH/nummaterial_velocity_file.reference b/waterlayered_halfspace/MESH/nummaterial_velocity_file.reference
new file mode 100644
index 0000000..2e08bc5
--- /dev/null
+++ b/waterlayered_halfspace/MESH/nummaterial_velocity_file.reference
@@ -0,0 +1,3 @@
+1   1 1028 1480 0 9999.0 9999.0 0  
+1   2 1028 1480 0 9999.0 9999.0 0 
+1   3 1028 1480 0 9999.0 9999.0 0
\ No newline at end of file
diff --git a/waterlayered_halfspace/MESH/nummaterial_velocity_file.reference.verticalboundary b/waterlayered_halfspace/MESH/nummaterial_velocity_file.reference.verticalboundary
new file mode 100644
index 0000000..cee0aa8
--- /dev/null
+++ b/waterlayered_halfspace/MESH/nummaterial_velocity_file.reference.verticalboundary
@@ -0,0 +1,2 @@
+1   1 2300 3000 0 9999.0 9999.0 0
+2   2 2300 3000 0.1 9999.0 9999.0 0
\ No newline at end of file
diff --git a/waterlayered_halfspace/boundary_definition.py b/waterlayered_halfspace/boundary_definition.py
deleted file mode 120000
index d4dcc5f..0000000
--- a/waterlayered_halfspace/boundary_definition.py
+++ /dev/null
@@ -1 +0,0 @@
-../../CUBIT_GEOCUBIT/geocubitlib/boundary_definition.py
\ No newline at end of file
diff --git a/waterlayered_halfspace/cubit2specfem3d.py b/waterlayered_halfspace/cubit2specfem3d.py
deleted file mode 120000
index f70980c..0000000
--- a/waterlayered_halfspace/cubit2specfem3d.py
+++ /dev/null
@@ -1 +0,0 @@
-../../CUBIT_GEOCUBIT/geocubitlib/cubit2specfem3d.py
\ No newline at end of file
diff --git a/waterlayered_halfspace/make_mesh_waterlater_verticalboundary.py b/waterlayered_halfspace/make_mesh_waterlater_verticalboundary.py
new file mode 100755
index 0000000..cd4c2c0
--- /dev/null
+++ b/waterlayered_halfspace/make_mesh_waterlater_verticalboundary.py
@@ -0,0 +1,33 @@
+#!/usr/bin/env python
+
+SEMoutput='MESH'
+CUBIToutput='MESH_GEOCUBIT'
+
+import cubit
+try:
+    cubit.init([""])
+except:
+    pass
+    
+
+
+cubit.cmd('reset')
+cubit.cmd('brick x 67000 y 134000 z 60000')
+cubit.cmd('volume 1 move x 33500 y 67000 z -30000')
+cubit.cmd('brick x 67000 y 134000 z 60000')
+cubit.cmd('volume 2 move x 100500 y 67000 z -30000')
+cubit.cmd('merge all')
+
+# Meshing the volumes
+elementsize = 3000.0
+
+cubit.cmd('volume 1 size '+str(elementsize))
+cubit.cmd('volume 2 size '+str(elementsize))
+cubit.cmd('mesh volume 1 2')
+
+
+from geocubitlib import boundary_definition,exportlib
+
+boundary_definition.define_bc(parallel=True)
+exportlib.collect(outdir=CUBIToutput)
+exportlib.e2SEM(outdir=SEMoutput)
diff --git a/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/make_mesh.sh b/waterlayered_halfspace/make_mesh_waterlayer.sh
similarity index 72%
copy from homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/make_mesh.sh
copy to waterlayered_halfspace/make_mesh_waterlayer.sh
index 4c28b8d..500ee7f 100755
--- a/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/make_mesh.sh
+++ b/waterlayered_halfspace/make_mesh_waterlayer.sh
@@ -11,9 +11,9 @@
 #    which GEOCUBIT.py
 #
 # 2. run this script to generate mesh
-#    ./make_mesh.sh
+#    ./make_mesh_waterlayer.sh
 #
 
-GEOCUBIT.py --build_volume --mesh --cfg=homogeneous_halfspace.cfg
+GEOCUBIT.py --build_volume --mesh --cfg=waterlayer_only.cfg
 GEOCUBIT.py --collect --meshfiles=MESH_GEOCUBIT/mesh_vol_0.e --export2SPECFEM3D --SEMoutput=MESH
-cp MESH-default/nummaterial_velocity_file.reference MESH/nummaterial_velocity_file
\ No newline at end of file
+cp MESH/nummaterial_velocity_file.reference MESH/nummaterial_velocity_file
\ No newline at end of file
diff --git a/layered_halfspace/layered_halfspace_notripling.cfg b/waterlayered_halfspace/waterlayer_only.cfg
similarity index 100%
copy from layered_halfspace/layered_halfspace_notripling.cfg
copy to waterlayered_halfspace/waterlayer_only.cfg