[cig-commits] r15433 - short/2.5D/benchmarks/savageprescott/meshes
willic3 at geodynamics.org
willic3 at geodynamics.org
Mon Jul 6 21:47:41 PDT 2009
Author: willic3
Date: 2009-07-06 21:47:40 -0700 (Mon, 06 Jul 2009)
New Revision: 15433
Added:
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_20km.exo.gz
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_20km.jou
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_20km.exo.gz
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_20km.jou
short/2.5D/benchmarks/savageprescott/meshes/spbm_tet4_graded_20km.exo.gz
short/2.5D/benchmarks/savageprescott/meshes/spbm_tet4_graded_20km.jou
short/2.5D/benchmarks/savageprescott/meshes/spbm_tet4_unif_20km.exo.gz
short/2.5D/benchmarks/savageprescott/meshes/spbm_tet4_unif_20km.jou
Removed:
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded1_20km.exo.gz
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded1_20km.jou
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded2_20km.exo.gz
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded2_20km.jou
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded3_20km.exo.gz
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded3_20km.jou
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded4_20km.exo.gz
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded4_20km.jou
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded5_20km.exo.gz
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded5_20km.jou
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif1_20km.exo.gz
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif1_20km.jou
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif2_20km.exo.gz
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif2_20km.jou
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif3_20km.exo.gz
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif3_20km.jou
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif4_20km.exo.gz
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif4_20km.jou
Log:
Deleted unused meshes and cleaned up existing hex mesh.
Created graded and uniform meshes for both hexes and tets.
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded1_20km.exo.gz
===================================================================
(Binary files differ)
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded1_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded1_20km.jou 2009-07-06 22:13:05 UTC (rev 15432)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded1_20km.jou 2009-07-07 04:47:40 UTC (rev 15433)
@@ -1,189 +0,0 @@
-## /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 500 km (x), 250 km (y), and 50 km (z)
-# is then put through 1 level of refinement.
-# The final mesh contains 78441 vertices and 70480 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 > -250000.0 and with x_coord < 250000.0 and with y_coord > -125000.0 and with y_coord < 125000.0 and with z_coord > -50000.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
-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"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_yneg_nofault" add node in face_yneg
-group "face_yneg_nofault" remove node in fault
-group "face_yneg_nofault" remove node in face_xpos
-group "face_yneg_nofault" remove node in face_xneg
-nodeset 18 face_yneg_nofault
-nodeset 18 name "face_yneg_nofault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_ypos_nofault" add node in face_ypos
-group "face_ypos_nofault" remove node in fault
-group "face_ypos_nofault" remove node in face_xpos
-group "face_ypos_nofault" remove node in face_xneg
-nodeset 19 face_ypos_nofault
-nodeset 19 name "face_ypos_nofault"
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "spbm_hex8_graded1_20km.exo" dimension 3 overwrite
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded2_20km.exo.gz
===================================================================
(Binary files differ)
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded2_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded2_20km.jou 2009-07-06 22:13:05 UTC (rev 15432)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded2_20km.jou 2009-07-07 04:47:40 UTC (rev 15433)
@@ -1,191 +0,0 @@
-## /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 500 km (x), 250 km (y), and 50 km (z)
-# is then put through 1 level of refinement.
-# The final mesh contains 78441 vertices and 70480 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 > -250000.0 and with x_coord < 250000.0 and with y_coord > -125000.0 and with y_coord < 125000.0 and with z_coord > -50000.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"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_yneg_nofault" add node in face_yneg
-group "face_yneg_nofault" remove node in fault
-group "face_yneg_nofault" remove node in face_xpos
-group "face_yneg_nofault" remove node in face_xneg
-nodeset 18 face_yneg_nofault
-nodeset 18 name "face_yneg_nofault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_ypos_nofault" add node in face_ypos
-group "face_ypos_nofault" remove node in fault
-group "face_ypos_nofault" remove node in face_xpos
-group "face_ypos_nofault" remove node in face_xneg
-nodeset 19 face_ypos_nofault
-nodeset 19 name "face_ypos_nofault"
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "spbm_hex8_graded2_20km.exo" dimension 3 overwrite
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded3_20km.exo.gz
===================================================================
(Binary files differ)
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded3_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded3_20km.jou 2009-07-06 22:13:05 UTC (rev 15432)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded3_20km.jou 2009-07-07 04:47:40 UTC (rev 15433)
@@ -1,189 +0,0 @@
-## /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 400 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 400 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 500 km (x), 250 km (y), and 100 km (z)
-# is then put through 1 level of refinement.
-# The final mesh contains 156633 vertices and 147110 cells.
-# ----------------------------------------------------------------------
-# Create outer block
-# ----------------------------------------------------------------------
-
-# Block is 2000 km x 1000 km x 400 km
-# -1000 km <= x <= 1000 km
-# -500 km <= y <= 500 km
-# -400 km <= z <= 0 km
-reset
-brick x 2000000 y 1000000 z 400000
-
-# Translate block so the top is at z=0
-volume 1 move x 0 y 0 z -200000
-
-# ----------------------------------------------------------------------
-# 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 > -240000.0 and with x_coord < 240000.0 and with y_coord > -120000.0 and with y_coord < 120000.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
-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"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_yneg_nofault" add node in face_yneg
-group "face_yneg_nofault" remove node in fault
-group "face_yneg_nofault" remove node in face_xpos
-group "face_yneg_nofault" remove node in face_xneg
-nodeset 18 face_yneg_nofault
-nodeset 18 name "face_yneg_nofault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_ypos_nofault" add node in face_ypos
-group "face_ypos_nofault" remove node in fault
-group "face_ypos_nofault" remove node in face_xpos
-group "face_ypos_nofault" remove node in face_xneg
-nodeset 19 face_ypos_nofault
-nodeset 19 name "face_ypos_nofault"
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "spbm_hex8_graded3_20km.exo" dimension 3 overwrite
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded4_20km.exo.gz
===================================================================
(Binary files differ)
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded4_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded4_20km.jou 2009-07-06 22:13:05 UTC (rev 15432)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded4_20km.jou 2009-07-07 04:47:40 UTC (rev 15433)
@@ -1,189 +0,0 @@
-## /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), 500 km (y),
-# and 800 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 800 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 500 km (x), 125 km (y), and 200 km (z)
-# is then put through 1 level of refinement.
-# The final mesh contains 159291 vertices and 150768 cells.
-# ----------------------------------------------------------------------
-# Create outer block
-# ----------------------------------------------------------------------
-
-# Block is 2000 km x 500 km x 800 km
-# -1000 km <= x <= 1000 km
-# -250 km <= y <= 250 km
-# -800 km <= z <= 0 km
-reset
-brick x 2000000 y 500000 z 800000
-
-# Translate block so the top is at z=0
-volume 1 move x 0 y 0 z -400000
-
-# ----------------------------------------------------------------------
-# 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 > -250000.0 and with x_coord < 250000.0 and with y_coord > -62500.0 and with y_coord < 62500.0 and with z_coord > -200000.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
-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"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_yneg_nofault" add node in face_yneg
-group "face_yneg_nofault" remove node in fault
-group "face_yneg_nofault" remove node in face_xpos
-group "face_yneg_nofault" remove node in face_xneg
-nodeset 18 face_yneg_nofault
-nodeset 18 name "face_yneg_nofault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_ypos_nofault" add node in face_ypos
-group "face_ypos_nofault" remove node in fault
-group "face_ypos_nofault" remove node in face_xpos
-group "face_ypos_nofault" remove node in face_xneg
-nodeset 19 face_ypos_nofault
-nodeset 19 name "face_ypos_nofault"
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "spbm_hex8_graded4_20km.exo" dimension 3 overwrite
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded5_20km.exo.gz
===================================================================
(Binary files differ)
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded5_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded5_20km.jou 2009-07-06 22:13:05 UTC (rev 15432)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded5_20km.jou 2009-07-07 04:47:40 UTC (rev 15433)
@@ -1,190 +0,0 @@
-## /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), 500 km (y),
-# and 800 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 800 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 500 km (x), 500 km (y), and 200 km (z)
-# is then put through 1 level of refinement.
-# Note that the entire region along the fault (y-axis) is refined.
-# The final mesh contains 307051 vertices and 294194 cells.
-# ----------------------------------------------------------------------
-# Create outer block
-# ----------------------------------------------------------------------
-
-# Block is 2000 km x 500 km x 800 km
-# -1000 km <= x <= 1000 km
-# -250 km <= y <= 250 km
-# -800 km <= z <= 0 km
-reset
-brick x 2000000 y 500000 z 800000
-
-# Translate block so the top is at z=0
-volume 1 move x 0 y 0 z -400000
-
-# ----------------------------------------------------------------------
-# 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 > -250000.0 and with x_coord < 250000.0 and with y_coord > -300000.0 and with y_coord < 300000.0 and with z_coord > -200000.0
-refine hex in inner_box depth 0
-
-# ----------------------------------------------------------------------
-# Smooth mesh.
-# ----------------------------------------------------------------------
-set debug 91 on
-volume all smooth scheme smart laplacian
-smooth volume all
-volume all smooth scheme condition number
-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
-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"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_yneg_nofault" add node in face_yneg
-group "face_yneg_nofault" remove node in fault
-group "face_yneg_nofault" remove node in face_xpos
-group "face_yneg_nofault" remove node in face_xneg
-nodeset 18 face_yneg_nofault
-nodeset 18 name "face_yneg_nofault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_ypos_nofault" add node in face_ypos
-group "face_ypos_nofault" remove node in fault
-group "face_ypos_nofault" remove node in face_xpos
-group "face_ypos_nofault" remove node in face_xneg
-nodeset 19 face_ypos_nofault
-nodeset 19 name "face_ypos_nofault"
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "spbm_hex8_graded5_20km.exo" dimension 3 overwrite
Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_20km.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_20km.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/x-gzip
Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_20km.jou (rev 0)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_20km.jou 2009-07-07 04:47:40 UTC (rev 15433)
@@ -0,0 +1,141 @@
+## /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 uses the geometry from file geometry.jou to create a mesh
+# with dimensions of 2000 km (x), 1000 km (y), and 400 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 400 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.
+# 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 480 km (x), 240 km (y), and 100 km (z)
+# is then put through 1 level of refinement.
+# The final mesh contains 155061 vertices and 145584 cells.
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set interval size for outer surfaces to 20 km and generate initial mesh
+# ----------------------------------------------------------------------
+volume all size 20000.0
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Refine region inside inner box.
+# ----------------------------------------------------------------------
+group "inner_box" add hex with x_coord > -240000.0 and with x_coord < 240000.0 and with y_coord > -120000.0 and with y_coord < 120000.0 and with z_coord > -100000.0
+refine hex in inner_box depth 0
+
+# ----------------------------------------------------------------------
+# Smooth mesh.
+# ----------------------------------------------------------------------
+set debug 91 on
+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
+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"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face without fault or x-boundaries
+# ----------------------------------------------------------------------
+group "face_yneg_nofault" add node in face_yneg
+group "face_yneg_nofault" remove node in fault
+group "face_yneg_nofault" remove node in face_xpos
+group "face_yneg_nofault" remove node in face_xneg
+nodeset 18 face_yneg_nofault
+nodeset 18 name "face_yneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face without fault or x-boundaries
+# ----------------------------------------------------------------------
+group "face_ypos_nofault" add node in face_ypos
+group "face_ypos_nofault" remove node in fault
+group "face_ypos_nofault" remove node in face_xpos
+group "face_ypos_nofault" remove node in face_xneg
+nodeset 19 face_ypos_nofault
+nodeset 19 name "face_ypos_nofault"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "spbm_hex8_graded_20km.exo" dimension 3 overwrite
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif1_20km.exo.gz
===================================================================
(Binary files differ)
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif1_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif1_20km.jou 2009-07-06 22:13:05 UTC (rev 15432)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif1_20km.jou 2009-07-07 04:47:40 UTC (rev 15433)
@@ -1,181 +0,0 @@
-## /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
-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"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_yneg_nofault" add node in face_yneg
-group "face_yneg_nofault" remove node in fault
-group "face_yneg_nofault" remove node in face_xpos
-group "face_yneg_nofault" remove node in face_xneg
-nodeset 18 face_yneg_nofault
-nodeset 18 name "face_yneg_nofault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_ypos_nofault" add node in face_ypos
-group "face_ypos_nofault" remove node in fault
-group "face_ypos_nofault" remove node in face_xpos
-group "face_ypos_nofault" remove node in face_xneg
-nodeset 19 face_ypos_nofault
-nodeset 19 name "face_ypos_nofault"
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "spbm_hex8_unif1_20km.exo" dimension 3 overwrite
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif2_20km.exo.gz
===================================================================
(Binary files differ)
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif2_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif2_20km.jou 2009-07-06 22:13:05 UTC (rev 15432)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif2_20km.jou 2009-07-07 04:47:40 UTC (rev 15433)
@@ -1,183 +0,0 @@
-## /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"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_yneg_nofault" add node in face_yneg
-group "face_yneg_nofault" remove node in fault
-group "face_yneg_nofault" remove node in face_xpos
-group "face_yneg_nofault" remove node in face_xneg
-nodeset 18 face_yneg_nofault
-nodeset 18 name "face_yneg_nofault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_ypos_nofault" add node in face_ypos
-group "face_ypos_nofault" remove node in fault
-group "face_ypos_nofault" remove node in face_xpos
-group "face_ypos_nofault" remove node in face_xneg
-nodeset 19 face_ypos_nofault
-nodeset 19 name "face_ypos_nofault"
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "spbm_hex8_unif2_20km.exo" dimension 3 overwrite
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif3_20km.exo.gz
===================================================================
(Binary files differ)
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif3_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif3_20km.jou 2009-07-06 22:13:05 UTC (rev 15432)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif3_20km.jou 2009-07-07 04:47:40 UTC (rev 15433)
@@ -1,181 +0,0 @@
-## /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 400 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 400 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 108171 vertices and 100000 cells.
-# ----------------------------------------------------------------------
-# Create outer block
-# ----------------------------------------------------------------------
-
-# Block is 2000 km x 1000 km x 400 km
-# -1000 km <= x <= 1000 km
-# -500 km <= y <= 500 km
-# -400 km <= z <= 0 km
-reset
-brick x 2000000 y 1000000 z 400000
-
-# Translate block so the top is at z=0
-volume 1 move x 0 y 0 z -200000
-
-# ----------------------------------------------------------------------
-# 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
-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"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_yneg_nofault" add node in face_yneg
-group "face_yneg_nofault" remove node in fault
-group "face_yneg_nofault" remove node in face_xpos
-group "face_yneg_nofault" remove node in face_xneg
-nodeset 18 face_yneg_nofault
-nodeset 18 name "face_yneg_nofault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_ypos_nofault" add node in face_ypos
-group "face_ypos_nofault" remove node in fault
-group "face_ypos_nofault" remove node in face_xpos
-group "face_ypos_nofault" remove node in face_xneg
-nodeset 19 face_ypos_nofault
-nodeset 19 name "face_ypos_nofault"
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "spbm_hex8_unif3_20km.exo" dimension 3 overwrite
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif4_20km.exo.gz
===================================================================
(Binary files differ)
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif4_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif4_20km.jou 2009-07-06 22:13:05 UTC (rev 15432)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif4_20km.jou 2009-07-07 04:47:40 UTC (rev 15433)
@@ -1,181 +0,0 @@
-## /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), 500 km (y),
-# and 800 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 800 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 111807 vertices and 104000 cells.
-# ----------------------------------------------------------------------
-# Create outer block
-# ----------------------------------------------------------------------
-
-# Block is 2000 km x 500 km x 800 km
-# -1000 km <= x <= 1000 km
-# -250 km <= y <= 250 km
-# -800 km <= z <= 0 km
-reset
-brick x 2000000 y 500000 z 800000
-
-# Translate block so the top is at z=0
-volume 1 move x 0 y 0 z -400000
-
-# ----------------------------------------------------------------------
-# 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
-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"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_yneg_nofault" add node in face_yneg
-group "face_yneg_nofault" remove node in fault
-group "face_yneg_nofault" remove node in face_xpos
-group "face_yneg_nofault" remove node in face_xneg
-nodeset 18 face_yneg_nofault
-nodeset 18 name "face_yneg_nofault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_ypos_nofault" add node in face_ypos
-group "face_ypos_nofault" remove node in fault
-group "face_ypos_nofault" remove node in face_xpos
-group "face_ypos_nofault" remove node in face_xneg
-nodeset 19 face_ypos_nofault
-nodeset 19 name "face_ypos_nofault"
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "spbm_hex8_unif4_20km.exo" dimension 3 overwrite
Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_20km.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_20km.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/x-gzip
Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_20km.jou (rev 0)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_20km.jou 2009-07-07 04:47:40 UTC (rev 15433)
@@ -0,0 +1,126 @@
+## /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 uses the geometry from file geometry.jou to create a mesh
+# with dimensions of 2000 km (x), 1000 km (y), and 400 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 400 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.
+# The fault is restricted to the upper two (elastic) layers.
+#
+# The volume is meshed at a uniform resolution of 20 km.
+# The final mesh contains 108171 vertices and 100000 cells.
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set interval size for entire volume to 20 km and generate mesh
+# ----------------------------------------------------------------------
+volume all size 20000.0
+mesh 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
+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"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face without fault or x-boundaries
+# ----------------------------------------------------------------------
+group "face_yneg_nofault" add node in face_yneg
+group "face_yneg_nofault" remove node in fault
+group "face_yneg_nofault" remove node in face_xpos
+group "face_yneg_nofault" remove node in face_xneg
+nodeset 18 face_yneg_nofault
+nodeset 18 name "face_yneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face without fault or x-boundaries
+# ----------------------------------------------------------------------
+group "face_ypos_nofault" add node in face_ypos
+group "face_ypos_nofault" remove node in fault
+group "face_ypos_nofault" remove node in face_xpos
+group "face_ypos_nofault" remove node in face_xneg
+nodeset 19 face_ypos_nofault
+nodeset 19 name "face_ypos_nofault"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "spbm_hex8_unif_20km.exo" dimension 3 overwrite
Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_tet4_graded_20km.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/2.5D/benchmarks/savageprescott/meshes/spbm_tet4_graded_20km.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/x-gzip
Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_tet4_graded_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_tet4_graded_20km.jou (rev 0)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_tet4_graded_20km.jou 2009-07-07 04:47:40 UTC (rev 15433)
@@ -0,0 +1,155 @@
+## /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 uses the geometry from file geometry.jou to create a mesh
+# with dimensions of 2000 km (x), 1000 km (y), and 400 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 400 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.
+# 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 480 km (x), 240 km (y), and 100 km (z)
+# is then put through 1 level of refinement.
+# The final mesh contains 125760 vertices and 705181 cells.
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set initial interval size to 20 km and set volume meshing scheme to
+# tetmesh
+# ----------------------------------------------------------------------
+volume all size 20000.0
+volume all scheme tetmesh
+
+# ----------------------------------------------------------------------
+# Set surface meshing scheme to qtri with a map base scheme.
+# This gives us regularly-spaced vertices that should have the same
+# locations as those for the hex meshes.
+# ----------------------------------------------------------------------
+surface all scheme qtri base scheme map
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Mesh the volume.
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Refine region inside inner box.
+# ----------------------------------------------------------------------
+group "inner_box" add tet with x_coord > -240000.0 and with x_coord < 240000.0 and with y_coord > -120000.0 and with y_coord < 120000.0 and with z_coord > -100000.0
+refine tet in inner_box depth 0
+
+# ----------------------------------------------------------------------
+# Smooth the mesh.
+# ----------------------------------------------------------------------
+set debug 91 on
+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
+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"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face without fault or x-boundaries
+# ----------------------------------------------------------------------
+group "face_yneg_nofault" add node in face_yneg
+group "face_yneg_nofault" remove node in fault
+group "face_yneg_nofault" remove node in face_xpos
+group "face_yneg_nofault" remove node in face_xneg
+nodeset 18 face_yneg_nofault
+nodeset 18 name "face_yneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face without fault or x-boundaries
+# ----------------------------------------------------------------------
+group "face_ypos_nofault" add node in face_ypos
+group "face_ypos_nofault" remove node in fault
+group "face_ypos_nofault" remove node in face_xpos
+group "face_ypos_nofault" remove node in face_xneg
+nodeset 19 face_ypos_nofault
+nodeset 19 name "face_ypos_nofault"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "spbm_tet4_graded_20km.exo" dimension 3 overwrite
Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_tet4_unif_20km.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/2.5D/benchmarks/savageprescott/meshes/spbm_tet4_unif_20km.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/x-gzip
Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_tet4_unif_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_tet4_unif_20km.jou (rev 0)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_tet4_unif_20km.jou 2009-07-07 04:47:40 UTC (rev 15433)
@@ -0,0 +1,146 @@
+## /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 uses the geometry from file geometry.jou to create a mesh
+# with dimensions of 2000 km (x), 1000 km (y), and 400 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 400 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.
+# The fault is restricted to the upper two (elastic) layers.
+#
+# The volume is meshed at a uniform resolution of 20 km.
+# The final mesh contains 105334 vertices and 586720 cells.
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set interval size to 20 km and set volume meshing scheme to tetmesh
+# ----------------------------------------------------------------------
+volume all size 20000.0
+volume all scheme tetmesh
+
+# ----------------------------------------------------------------------
+# Set surface meshing scheme to qtri with a map base scheme.
+# This gives us regularly-spaced vertices that should have the same
+# locations as those for the hex meshes.
+# ----------------------------------------------------------------------
+surface all scheme qtri base scheme map
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Mesh the volume.
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Smooth the mesh.
+# ----------------------------------------------------------------------
+set debug 91 on
+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
+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"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face without fault or x-boundaries
+# ----------------------------------------------------------------------
+group "face_yneg_nofault" add node in face_yneg
+group "face_yneg_nofault" remove node in fault
+group "face_yneg_nofault" remove node in face_xpos
+group "face_yneg_nofault" remove node in face_xneg
+nodeset 18 face_yneg_nofault
+nodeset 18 name "face_yneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face without fault or x-boundaries
+# ----------------------------------------------------------------------
+group "face_ypos_nofault" add node in face_ypos
+group "face_ypos_nofault" remove node in fault
+group "face_ypos_nofault" remove node in face_xpos
+group "face_ypos_nofault" remove node in face_xneg
+nodeset 19 face_ypos_nofault
+nodeset 19 name "face_ypos_nofault"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "spbm_tet4_unif_20km.exo" dimension 3 overwrite
More information about the CIG-COMMITS
mailing list