[cig-commits] r14196 - short/2.5D/benchmarks/savageprescott/meshes

willic3 at geodynamics.org willic3 at geodynamics.org
Mon Mar 2 18:19:01 PST 2009 

Author: willic3
Date: 2009-03-02 18:19:00 -0800 (Mon, 02 Mar 2009)
New Revision: 14196

Added:
   short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded1_20km.jou
   short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded2_20km.jou
   short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif1_20km.exo
   short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif1_20km.jou
   short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif2_20km.exo
   short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif2_20km.jou
Log:
Added Cubit journal files along with the resulting meshes for the case of
a uniform mesh at 20 km resolution and a graded mesh with one level of
refinement.
The two graded meshes are still being generated and will be added shortly.



Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded1_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded1_20km.jou	                        (rev 0)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded1_20km.jou	2009-03-03 02:19:00 UTC (rev 14196)
@@ -0,0 +1,172 @@
+## /Applications/Plotting/Cubit_10_2.app/Contents/MacOS/Cubit_10_2
+## Cubit Version 10.2
+## Cubit Build 8
+## Revised 10/17/2006 15:13:41 MST
+## Running 06/16/2008 02:47:47 PM
+## Command Options:
+## -warning = On
+## -information = On
+
+# Description:
+# This journal creates a mesh with dimensions of 2000 km (x), 1000 km (y),
+# and 200 km (z).  The mesh is divided into 3 horizontal layers:
+# a surface layer from 0 to 20 km, marking the bottom of the locking depth;
+# a middle layer from 20 to 40 km, marking the creeping section;
+# a bottom layer from 40 to 200 km, which approximates the viscoelastic
+# half-space.
+# The mesh is cut vertically by the fault surface, which lies along the
+# x-axis, and by a symmetry plane, which lies along the y-axis.
+# In this version of the mesh, the fault is restricted to the upper two
+# (elastic) layers.
+# The volume is meshed at an initial resolution of 20 km.
+# An inner volume with dimensions of 1000 km (x), 500 km (y), and 100 km (z)
+# is then put through 1 level of refinement.
+# The final mesh contains 247405 vertices and 235482 cells.
+# ----------------------------------------------------------------------
+# Create outer block
+# ----------------------------------------------------------------------
+
+# Block is 2000 km x 1000 km x 200 km
+# -1000 km <= x <= 1000 km
+# -500 km <= y <= 500 km
+# -200 km <= z <= 0 km
+reset
+brick x 2000000 y 1000000 z 200000
+
+# Translate block so the top is at z=0
+volume 1 move x 0 y 0 z -100000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create planar surface with plane xplane offset 0
+surface 7 name "fault_surface"
+create planar surface with plane zplane offset -40000
+surface 8 name "material_interface"
+create planar surface with plane yplane offset 0
+surface 9 name "symmetry_plane"
+create planar surface with plane zplane offset -20000
+surface 10 name "fault_bottom"
+
+# ----------------------------------------------------------------------
+# Cut volume with interfaces.
+# ----------------------------------------------------------------------
+compress all
+webcut volume all with plane fault_surface
+webcut volume all with plane material_interface
+webcut volume all with plane symmetry_plane
+webcut volume all with plane fault_bottom
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces.
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+
+# ----------------------------------------------------------------------
+# Remove surface bodies.
+# ----------------------------------------------------------------------
+delete body 2 3 4 5 7 8 9 11 12 14 16 17 18 20 21 23 26 27 29 31
+imprint all with volume all
+merge all
+
+# ----------------------------------------------------------------------
+# Set interval size for outer surfaces to 20 km
+# ----------------------------------------------------------------------
+volume all size 20000.0
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Refine region inside inner box.
+# ----------------------------------------------------------------------
+group "inner_box" add hex with x_coord > -500000.0 and with x_coord < 500000.0 and with y_coord > -250000.0 and with y_coord < 250000.0 and with z_coord > -100000.0
+refine hex in inner_box depth 0
+
+# ----------------------------------------------------------------------
+# Smooth mesh.
+# ----------------------------------------------------------------------
+volume all smooth scheme smart laplacian
+smooth volume all
+volume all smooth scheme condition number
+set debug 91 on
+smooth volume all
+volume all smooth scheme mean ratio
+smooth volume all
+volume all smooth scheme condition number
+smooth volume all
+
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1 6 15 19 25 28 30 32
+block 1 name "elastic"
+block 2 volume 10 13 22 24
+block 2 name "viscoelastic"
+
+# ----------------------------------------------------------------------
+# Create fault nodeset
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+group "fault" add node in fault_surface at A
+group "fault" add node in fault_surface at D
+group "fault" add node in fault_surface at E
+# Leave out the portion of the fault surface in the viscoelastic layer
+# group "fault" add node in fault_surface at B
+# group "fault" add node in fault_surface at C
+nodeset 10 fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 133 109 139 111 89 87
+nodeset 11 face_xpos
+nodeset 11 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 121 147 127 149 97 104
+nodeset 12 face_xneg
+nodeset 12 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 106 120 113 128 35 47
+nodeset 13 face_ypos
+nodeset 13 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 144 132 151 140 49 33
+nodeset 14 face_yneg
+nodeset 14 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 103 92 96 84
+nodeset 15 face_zneg
+nodeset 15 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face without fault
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 16 face_zneg_nofault
+nodeset 16 name "face_zneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 78 62 70 54
+nodeset 17 face_zpos
+nodeset 17 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "spbm_hex8_graded1_20km.exo" dimension 3 overwrite

Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded2_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded2_20km.jou	                        (rev 0)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded2_20km.jou	2009-03-03 02:19:00 UTC (rev 14196)
@@ -0,0 +1,171 @@
+## /Applications/Plotting/Cubit_10_2.app/Contents/MacOS/Cubit_10_2
+## Cubit Version 10.2
+## Cubit Build 8
+## Revised 10/17/2006 15:13:41 MST
+## Running 06/16/2008 02:47:47 PM
+## Command Options:
+## -warning = On
+## -information = On
+
+# Description:
+# This journal creates a mesh with dimensions of 2000 km (x), 1000 km (y),
+# and 200 km (z).  The mesh is divided into 3 horizontal layers:
+# a surface layer from 0 to 20 km, marking the bottom of the locking depth;
+# a middle layer from 20 to 40 km, marking the creeping section;
+# a bottom layer from 40 to 200 km, which approximates the viscoelastic
+# half-space.
+# The mesh is cut vertically by the fault surface, which lies along the
+# x-axis, and by a symmetry plane, which lies along the y-axis.
+# In this version of the mesh, the fault cuts through the entire mesh.
+# The volume is meshed at an initial resolution of 20 km.
+# An inner volume with dimensions of 1000 km (x), 500 km (y), and 100 km (z)
+# is then put through 1 level of refinement.
+# The final mesh contains 247405 vertices and 235482 cells.
+# ----------------------------------------------------------------------
+# Create outer block
+# ----------------------------------------------------------------------
+
+# Block is 2000 km x 1000 km x 200 km
+# -1000 km <= x <= 1000 km
+# -500 km <= y <= 500 km
+# -200 km <= z <= 0 km
+reset
+brick x 2000000 y 1000000 z 200000
+
+# Translate block so the top is at z=0
+volume 1 move x 0 y 0 z -100000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create planar surface with plane xplane offset 0
+surface 7 name "fault_surface"
+create planar surface with plane zplane offset -40000
+surface 8 name "material_interface"
+create planar surface with plane yplane offset 0
+surface 9 name "symmetry_plane"
+create planar surface with plane zplane offset -20000
+surface 10 name "fault_bottom"
+
+# ----------------------------------------------------------------------
+# Cut volume with interfaces.
+# ----------------------------------------------------------------------
+compress all
+webcut volume all with plane fault_surface
+webcut volume all with plane material_interface
+webcut volume all with plane symmetry_plane
+webcut volume all with plane fault_bottom
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces.
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+
+# ----------------------------------------------------------------------
+# Remove surface bodies.
+# ----------------------------------------------------------------------
+delete body 2 3 4 5 7 8 9 11 12 14 16 17 18 20 21 23 26 27 29 31
+imprint all with volume all
+merge all
+
+# ----------------------------------------------------------------------
+# Set interval size for outer surfaces to 20 km
+# ----------------------------------------------------------------------
+volume all size 20000.0
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Refine region inside inner box.
+# ----------------------------------------------------------------------
+group "inner_box" add hex with x_coord > -500000.0 and with x_coord < 500000.0 and with y_coord > -250000.0 and with y_coord < 250000.0 and with z_coord > -100000.0
+refine hex in inner_box depth 0
+
+# ----------------------------------------------------------------------
+# Smooth mesh.
+# ----------------------------------------------------------------------
+volume all smooth scheme smart laplacian
+smooth volume all
+volume all smooth scheme condition number
+set debug 91 on
+smooth volume all
+volume all smooth scheme mean ratio
+smooth volume all
+volume all smooth scheme condition number
+smooth volume all
+
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1 6 15 19 25 28 30 32
+block 1 name "elastic"
+block 2 volume 10 13 22 24
+block 2 name "viscoelastic"
+
+# ----------------------------------------------------------------------
+# Create fault nodeset
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+group "fault" add node in fault_surface at A
+group "fault" add node in fault_surface at D
+group "fault" add node in fault_surface at E
+# Include portion of fault mesh in the viscoelastic layer.
+group "fault" add node in fault_surface at B
+group "fault" add node in fault_surface at C
+nodeset 10 fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 133 109 139 111 89 87
+nodeset 11 face_xpos
+nodeset 11 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 121 147 127 149 97 104
+nodeset 12 face_xneg
+nodeset 12 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 106 120 113 128 35 47
+nodeset 13 face_ypos
+nodeset 13 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 144 132 151 140 49 33
+nodeset 14 face_yneg
+nodeset 14 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 103 92 96 84
+nodeset 15 face_zneg
+nodeset 15 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face without fault
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 16 face_zneg_nofault
+nodeset 16 name "face_zneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 78 62 70 54
+nodeset 17 face_zpos
+nodeset 17 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "spbm_hex8_graded2_20km.exo" dimension 3 overwrite

Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif1_20km.exo
===================================================================
(Binary files differ)


Property changes on: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif1_20km.exo
___________________________________________________________________
Name: svn:mime-type
   + application/octet-stream

Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif1_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif1_20km.jou	                        (rev 0)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif1_20km.jou	2009-03-03 02:19:00 UTC (rev 14196)
@@ -0,0 +1,164 @@
+## /Applications/Plotting/Cubit_10_2.app/Contents/MacOS/Cubit_10_2
+## Cubit Version 10.2
+## Cubit Build 8
+## Revised 10/17/2006 15:13:41 MST
+## Running 06/16/2008 02:47:47 PM
+## Command Options:
+## -warning = On
+## -information = On
+
+# Description:
+# This journal creates a mesh with dimensions of 2000 km (x), 1000 km (y),
+# and 200 km (z).  The mesh is divided into 3 horizontal layers:
+# a surface layer from 0 to 20 km, marking the bottom of the locking depth;
+# a middle layer from 20 to 40 km, marking the creeping section;
+# a bottom layer from 40 to 200 km, which approximates the viscoelastic
+# half-space.
+# The mesh is cut vertically by the fault surface, which lies along the
+# x-axis, and by a symmetry plane, which lies along the y-axis.
+# In this version of the mesh, the fault is restricted to the upper two
+# (elastic) layers.
+# The entire volume is meshed at a resolution of 20 km.
+# The final mesh contains 56661 vertices and 50000 cells.
+# ----------------------------------------------------------------------
+# Create outer block
+# ----------------------------------------------------------------------
+
+# Block is 2000 km x 1000 km x 200 km
+# -1000 km <= x <= 1000 km
+# -500 km <= y <= 500 km
+# -200 km <= z <= 0 km
+reset
+brick x 2000000 y 1000000 z 200000
+
+# Translate block so the top is at z=0
+volume 1 move x 0 y 0 z -100000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create planar surface with plane xplane offset 0
+surface 7 name "fault_surface"
+create planar surface with plane zplane offset -40000
+surface 8 name "material_interface"
+create planar surface with plane yplane offset 0
+surface 9 name "symmetry_plane"
+create planar surface with plane zplane offset -20000
+surface 10 name "fault_bottom"
+
+# ----------------------------------------------------------------------
+# Cut volume with interfaces.
+# ----------------------------------------------------------------------
+compress all
+webcut volume all with plane fault_surface
+webcut volume all with plane material_interface
+webcut volume all with plane symmetry_plane
+webcut volume all with plane fault_bottom
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces.
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+
+# ----------------------------------------------------------------------
+# Remove surface bodies.
+# ----------------------------------------------------------------------
+delete body 2 3 4 5 7 8 9 11 12 14 16 17 18 20 21 23 26 27 29 31
+imprint all with volume all
+merge all
+
+# ----------------------------------------------------------------------
+# Set interval size for outer surfaces to 20 km
+# ----------------------------------------------------------------------
+volume all size 20000.0
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Smooth mesh.
+# ----------------------------------------------------------------------
+volume all smooth scheme smart laplacian
+smooth volume all
+volume all smooth scheme condition number
+set debug 91 on
+smooth volume all
+volume all smooth scheme mean ratio
+smooth volume all
+volume all smooth scheme condition number
+smooth volume all
+
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1 6 15 19 25 28 30 32
+block 1 name "elastic"
+block 2 volume 10 13 22 24
+block 2 name "viscoelastic"
+
+# ----------------------------------------------------------------------
+# Create fault nodeset
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+group "fault" add node in fault_surface at A
+group "fault" add node in fault_surface at D
+group "fault" add node in fault_surface at E
+# Leave out the portion of the fault surface in the viscoelastic layer
+# group "fault" add node in fault_surface at B
+# group "fault" add node in fault_surface at C
+nodeset 10 fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 133 109 139 111 89 87
+nodeset 11 face_xpos
+nodeset 11 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 121 147 127 149 97 104
+nodeset 12 face_xneg
+nodeset 12 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 106 120 113 128 35 47
+nodeset 13 face_ypos
+nodeset 13 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 144 132 151 140 49 33
+nodeset 14 face_yneg
+nodeset 14 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 103 92 96 84
+nodeset 15 face_zneg
+nodeset 15 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face without fault
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 16 face_zneg_nofault
+nodeset 16 name "face_zneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 78 62 70 54
+nodeset 17 face_zpos
+nodeset 17 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "spbm_hex8_unif1_20km.exo" dimension 3 overwrite

Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif2_20km.exo
===================================================================
(Binary files differ)


Property changes on: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif2_20km.exo
___________________________________________________________________
Name: svn:mime-type
   + application/octet-stream

Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif2_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif2_20km.jou	                        (rev 0)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif2_20km.jou	2009-03-03 02:19:00 UTC (rev 14196)
@@ -0,0 +1,163 @@
+## /Applications/Plotting/Cubit_10_2.app/Contents/MacOS/Cubit_10_2
+## Cubit Version 10.2
+## Cubit Build 8
+## Revised 10/17/2006 15:13:41 MST
+## Running 06/16/2008 02:47:47 PM
+## Command Options:
+## -warning = On
+## -information = On
+
+# Description:
+# This journal creates a mesh with dimensions of 2000 km (x), 1000 km (y),
+# and 200 km (z).  The mesh is divided into 3 horizontal layers:
+# a surface layer from 0 to 20 km, marking the bottom of the locking depth;
+# a middle layer from 20 to 40 km, marking the creeping section;
+# a bottom layer from 40 to 200 km, which approximates the viscoelastic
+# half-space.
+# The mesh is cut vertically by the fault surface, which lies along the
+# x-axis, and by a symmetry plane, which lies along the y-axis.
+# In this version of the mesh, the fault cuts through the entire mesh.
+# The entire volume is meshed at a resolution of 20 km.
+# The final mesh contains 56661 vertices and 50000 cells.
+# ----------------------------------------------------------------------
+# Create outer block
+# ----------------------------------------------------------------------
+
+# Block is 2000 km x 1000 km x 200 km
+# -1000 km <= x <= 1000 km
+# -500 km <= y <= 500 km
+# -200 km <= z <= 0 km
+reset
+brick x 2000000 y 1000000 z 200000
+
+# Translate block so the top is at z=0
+volume 1 move x 0 y 0 z -100000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create planar surface with plane xplane offset 0
+surface 7 name "fault_surface"
+create planar surface with plane zplane offset -40000
+surface 8 name "material_interface"
+create planar surface with plane yplane offset 0
+surface 9 name "symmetry_plane"
+create planar surface with plane zplane offset -20000
+surface 10 name "fault_bottom"
+
+# ----------------------------------------------------------------------
+# Cut volume with interfaces.
+# ----------------------------------------------------------------------
+compress all
+webcut volume all with plane fault_surface
+webcut volume all with plane material_interface
+webcut volume all with plane symmetry_plane
+webcut volume all with plane fault_bottom
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces.
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+
+# ----------------------------------------------------------------------
+# Remove surface bodies.
+# ----------------------------------------------------------------------
+delete body 2 3 4 5 7 8 9 11 12 14 16 17 18 20 21 23 26 27 29 31
+imprint all with volume all
+merge all
+
+# ----------------------------------------------------------------------
+# Set interval size for outer surfaces to 20 km
+# ----------------------------------------------------------------------
+volume all size 20000.0
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Smooth mesh.
+# ----------------------------------------------------------------------
+volume all smooth scheme smart laplacian
+smooth volume all
+volume all smooth scheme condition number
+set debug 91 on
+smooth volume all
+volume all smooth scheme mean ratio
+smooth volume all
+volume all smooth scheme condition number
+smooth volume all
+
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1 6 15 19 25 28 30 32
+block 1 name "elastic"
+block 2 volume 10 13 22 24
+block 2 name "viscoelastic"
+
+# ----------------------------------------------------------------------
+# Create fault nodeset
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+group "fault" add node in fault_surface at A
+group "fault" add node in fault_surface at D
+group "fault" add node in fault_surface at E
+# Include portion of fault mesh in the viscoelastic layer.
+group "fault" add node in fault_surface at B
+group "fault" add node in fault_surface at C
+nodeset 10 fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 133 109 139 111 89 87
+nodeset 11 face_xpos
+nodeset 11 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 121 147 127 149 97 104
+nodeset 12 face_xneg
+nodeset 12 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 106 120 113 128 35 47
+nodeset 13 face_ypos
+nodeset 13 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 144 132 151 140 49 33
+nodeset 14 face_yneg
+nodeset 14 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 103 92 96 84
+nodeset 15 face_zneg
+nodeset 15 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face without fault
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 16 face_zneg_nofault
+nodeset 16 name "face_zneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 78 62 70 54
+nodeset 17 face_zpos
+nodeset 17 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "spbm_hex8_unif2_20km.exo" dimension 3 overwrite

More information about the CIG-COMMITS mailing list