[cig-commits] [commit] devel, master: updated homogeneous_halfspace example (09cae29)

[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 09cae292cabcb96e7734b35f2b5bae6df56ae124
Author: Carl Tape <carltape at gi.alaska.edu>
Date:   Thu Oct 3 22:47:16 2013 +0000

    updated homogeneous_halfspace example


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

09cae292cabcb96e7734b35f2b5bae6df56ae124
 .../MESH-default/nummaterial_velocity_file         |  2 +-
 .../nummaterial_velocity_file.original             |  1 +
 .../nummaterial_velocity_file.reference            |  1 -
 .../README                                         | 32 ++++++++++++++++------
 .../make_mesh.sh                                   | 16 +----------
 5 files changed, 27 insertions(+), 25 deletions(-)

diff --git a/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file b/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file
index bcb1ce4..67fba13 100644
--- a/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file
+++ b/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file
@@ -1 +1 @@
-1   1          2300.000000          2800.000000          0.000000    9999.00000  9999.00000  0
\ No newline at end of file
+2   1          2300.000000          2800.000000          1500.000000          9999.000000          9999.000000  0
\ No newline at end of file
diff --git a/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.original b/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.original
new file mode 100644
index 0000000..5c872f6
--- /dev/null
+++ b/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.original
@@ -0,0 +1 @@
+0   1 vol4  -->       sintax: #material_domain_id #material_id #rho #vp #vs #Q_kappa #Q_mu #anisotropy
diff --git a/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.reference b/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.reference
deleted file mode 100644
index 67fba13..0000000
--- a/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.reference
+++ /dev/null
@@ -1 +0,0 @@
-2   1          2300.000000          2800.000000          1500.000000          9999.000000          9999.000000  0
\ No newline at end of file
diff --git a/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/README b/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/README
index 8d16a4f..53e16e1 100644
--- a/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/README
+++ b/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/README
@@ -2,16 +2,12 @@
 README
 ----------------------------------------------------------------------
 
+This example creates a homogeneous halfspace, i.e. a single volume block with a
+constant elastic material property, using CUBIT, and runs a forward simulation.
 
-Note: to create the mesh again using CUBIT, just type " claro ./block_mesh.py "
-("claro" is the way of calling CUBIT from the command line)
-
-
-
-This examples creates a homogeneous halfspace, i.e. a single volume block with a
-constant elastic material property, using CUBIT and runs a forward simulation.
+----------------------------------------------------------------------
 
-step-by-step tutorial:
+step-by-step tutorial (short version):
 
 0. check that all software is available (or that modules are loaded):
 	openmpi: > which mpirun
@@ -22,6 +18,20 @@ step-by-step tutorial:
         gnuplot: > which gnuplot
         xmgrace: > which xmgrace
 
+1. check paths and set paths for meshing with GEOCUBIT
+
+    - From the SPECFEM3D 'base' directory SPECFEM3D/
+      > cd CUBIT_GEOCUBIT
+      > source setpaths.sh
+
+2. run process.sh
+
+----------------------------------------------------------------------
+
+step-by-step tutorial (long version):
+
+0. perform steps 0. and 1. above
+
 1. configure package:
 
    - determine the path for your compiled version of SCOTCH
@@ -54,6 +64,10 @@ step-by-step tutorial:
 
 3. create mesh:
 
+   - check path, then make mesh
+     > which GEOCUBIT.py
+     > make_mesh.sh
+
    - see instructions in make_mesh.sh
    
      if everything goes fine, this creates the ten mesh files in subdirectory MESH/:
@@ -68,10 +82,12 @@ step-by-step tutorial:
         nodes_coords_file
         nummaterial_velocity_file
 
+     (there will also be several local files generated)
      the mesh can be viewed using either cubit or paraview:
 
      > claro &
      then open TOTALMESH_MERGED.cub
+     optional: type "draw block all" in the cubit command line window
 
      > paraview &
      then open TOTALMESH_MERGED.e
diff --git a/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/make_mesh.sh b/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/make_mesh.sh
index 4c28b8d..fb4f585 100755
--- a/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/make_mesh.sh
+++ b/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/make_mesh.sh
@@ -1,19 +1,5 @@
 #!/bin/bash
 
-#
-# Temporary instructions
-#
-# 1. set path to GEOCUBIT base directory (../../CUBIT_GEOCUBIT/), for example:
-#    export PYTHONPATH=$PYTHONPATH:/import/c/w/tape/3D/SPECFEM3D/CUBIT_GEOCUBIT
-#    export PATH=$PATH:/import/c/w/tape/3D/SPECFEM3D/CUBIT_GEOCUBIT
-#
-#    check path:
-#    which GEOCUBIT.py
-#
-# 2. run this script to generate mesh
-#    ./make_mesh.sh
-#
-
 GEOCUBIT.py --build_volume --mesh --cfg=homogeneous_halfspace.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-default/nummaterial_velocity_file MESH/
\ No newline at end of file