[cig-commits] r22920 - in seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides: . MESH-default
carltape at geodynamics.org
carltape at geodynamics.org
Thu Oct 3 15:47:16 PDT 2013
Author: carltape
Date: 2013-10-03 15:47:16 -0700 (Thu, 03 Oct 2013)
New Revision: 22920
Added:
seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file
seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.original
Removed:
seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file
seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.reference
Modified:
seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/README
seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/make_mesh.sh
Log:
updated homogeneous_halfspace example
Deleted: seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file 2013-10-03 09:53:40 UTC (rev 22919)
+++ seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file 2013-10-03 22:47:16 UTC (rev 22920)
@@ -1 +0,0 @@
-1 1 2300.000000 2800.000000 0.000000 9999.00000 9999.00000 0
\ No newline at end of file
Copied: seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file (from rev 22919, seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.reference)
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file 2013-10-03 22:47:16 UTC (rev 22920)
@@ -0,0 +1 @@
+2 1 2300.000000 2800.000000 1500.000000 9999.000000 9999.000000 0
\ No newline at end of file
Copied: seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.original (from rev 22919, seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file)
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.original (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.original 2013-10-03 22:47:16 UTC (rev 22920)
@@ -0,0 +1 @@
+0 1 vol4 --> sintax: #material_domain_id #material_id #rho #vp #vs #Q_kappa #Q_mu #anisotropy
Deleted: seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.reference
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.reference 2013-10-03 09:53:40 UTC (rev 22919)
+++ seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/MESH-default/nummaterial_velocity_file.reference 2013-10-03 22:47:16 UTC (rev 22920)
@@ -1 +0,0 @@
-2 1 2300.000000 2800.000000 1500.000000 9999.000000 9999.000000 0
\ No newline at end of file
Modified: seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/README
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/README 2013-10-03 09:53:40 UTC (rev 22919)
+++ seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/README 2013-10-03 22:47:16 UTC (rev 22920)
@@ -2,17 +2,13 @@
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)
+----------------------------------------------------------------------
+step-by-step tutorial (short version):
-
-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:
-
0. check that all software is available (or that modules are loaded):
openmpi: > which mpirun
cubit: > which cubit
@@ -22,6 +18,20 @@
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 @@
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 @@
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
Modified: seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/make_mesh.sh
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/make_mesh.sh 2013-10-03 09:53:40 UTC (rev 22919)
+++ seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8_elastic_absorbing_Stacey_5sides/make_mesh.sh 2013-10-03 22:47:16 UTC (rev 22920)
@@ -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
More information about the CIG-COMMITS
mailing list