[cig-commits] r18928 - short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes
willic3 at geodynamics.org
willic3 at geodynamics.org
Sat Sep 17 11:58:29 PDT 2011
Author: willic3
Date: 2011-09-17 11:58:29 -0700 (Sat, 17 Sep 2011)
New Revision: 18928
Added:
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.inp.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.jou
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_graded_4km.exo.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_graded_4km.inp.gz
Removed:
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_graded_20km.exo.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_graded_20km.inp.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_refined_20km.jou
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_unif_20km.exo.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_unif_20km.inp.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_unif_20km.jou
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_graded_20km.exo.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_graded_20km.inp.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_refine1_20km.jou
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_unif_20km.exo.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_unif_20km.inp.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_unif_20km.jou
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_4km.exo.gz
Modified:
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.exo.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_graded_4km.jou
Log:
Removed all the uniform meshes and previous meshes.
Now all we have is a graded tet mesh and a hex mesh with one level of
refinement in the center.
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.exo.gz
===================================================================
(Binary files differ)
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.inp.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.inp.gz
___________________________________________________________________
Name: svn:mime-type
+ application/x-gzip
Copied: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.jou (from rev 18926, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_refined_20km.jou)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.jou 2011-09-17 18:58:29 UTC (rev 18928)
@@ -0,0 +1,132 @@
+# 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 2 horizontal layers:
+# a surface layer from 0 to 40 km, marking elastic layer;
+# a bottom layer from 40 to 400 km, which approximates the viscoelastic
+# half-space.
+# The fault plane is also cut horizontally at 20 km depth, representing
+# the division between the upper locked section and the lower creeping
+# section.
+#
+# 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 (elastic) layer.
+#
+# 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 beta 2.2
+smooth volume all
+
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1 6 15 19
+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
+group "fault" add node in surface 85
+group "fault" add node in surface 91
+group "fault" remove node with z_coord < -82000.0
+nodeset 10 fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 61 55 89 87
+nodeset 11 face_xpos
+nodeset 11 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 73 75 97 104
+nodeset 12 face_xneg
+nodeset 12 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 29 45 35 47
+nodeset 13 face_ypos
+nodeset 13 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 43 31 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
+# ----------------------------------------------------------------------
+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 "hex8_refined_20km.exo" dimension 3 overwrite
+export Abaqus "hex8_refined_20km.inp"
+
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_graded_20km.exo.gz
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_graded_20km.inp.gz
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_refined_20km.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_refined_20km.jou 2011-09-17 18:48:29 UTC (rev 18927)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_refined_20km.jou 2011-09-17 18:58:29 UTC (rev 18928)
@@ -1,127 +0,0 @@
-# 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 2 horizontal layers:
-# a surface layer from 0 to 40 km, marking elastic layer;
-# a bottom layer from 40 to 400 km, which approximates the viscoelastic
-# half-space.
-# The fault plane is also cut horizontally at 20 km depth, representing
-# the division between the upper locked section and the lower creeping
-# section.
-#
-# 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 (elastic) layer.
-#
-# 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
-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 61 55 89 87
-nodeset 11 face_xpos
-nodeset 11 name "face_xpos"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -x face
-# ----------------------------------------------------------------------
-group "face_xneg" add node in surface 73 75 97 104
-nodeset 12 face_xneg
-nodeset 12 name "face_xneg"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face
-# ----------------------------------------------------------------------
-group "face_ypos" add node in surface 29 45 35 47
-nodeset 13 face_ypos
-nodeset 13 name "face_ypos"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face
-# ----------------------------------------------------------------------
-group "face_yneg" add node in surface 43 31 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
-# ----------------------------------------------------------------------
-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 "hex8_refined_20km.exo" dimension 3 overwrite
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_unif_20km.exo.gz
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_unif_20km.inp.gz
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_unif_20km.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_unif_20km.jou 2011-09-17 18:48:29 UTC (rev 18927)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_hex8_unif_20km.jou 2011-09-17 18:58:29 UTC (rev 18928)
@@ -1,127 +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 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
-export Abaqus "spbm_hex8_unif_20km.inp"
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_graded_20km.exo.gz
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_graded_20km.inp.gz
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_refine1_20km.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_refine1_20km.jou 2011-09-17 18:48:29 UTC (rev 18927)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_refine1_20km.jou 2011-09-17 18:58:29 UTC (rev 18928)
@@ -1,156 +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 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
-export Abaqus "spbm_tet4_graded_20km.inp"
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_unif_20km.exo.gz
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_unif_20km.inp.gz
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_unif_20km.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_unif_20km.jou 2011-09-17 18:48:29 UTC (rev 18927)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_unif_20km.jou 2011-09-17 18:58:29 UTC (rev 18928)
@@ -1,147 +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 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
-export Abaqus "spbm_tet4_unif_20km.inp"
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_4km.exo.gz
===================================================================
(Binary files differ)
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_graded_4km.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_graded_4km.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/x-gzip
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_graded_4km.inp.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_graded_4km.inp.gz
___________________________________________________________________
Name: svn:mime-type
+ application/x-gzip
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_graded_4km.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_graded_4km.jou 2011-09-17 18:48:29 UTC (rev 18927)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_graded_4km.jou 2011-09-17 18:58:29 UTC (rev 18928)
@@ -137,6 +137,9 @@
group "fault" add node in fault_surface at D
group "fault" add node in fault_surface at A
group "fault" add node in fault_surface at E
+group "fault" add node in surface 85
+group "fault" add node in surface 91
+group "fault" remove node with z_coord < -82000.0
nodeset 10 fault
nodeset 10 name "fault"
@@ -205,5 +208,7 @@
# ----------------------------------------------------------------------
# Export exodus file
# ----------------------------------------------------------------------
-export mesh "tet4_4km.exo" dimension 3 overwrite
+export mesh "tet4_graded_4km.exo" dimension 3 overwrite
+export Abaqus "tet4_graded_4km.inp"
+
More information about the CIG-COMMITS
mailing list